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ABSTRACT 


This Is the third report In a series of three workshops, sponsored by the 
National Aeronautics and Space \dmi?ilscrat Ion , to investigate the state of the 
art In global sea surface temperature measurements from space. Three workshops 
were necessary to process and analyze sufficient dat3 from which to draw con- 
clusions on the accuracy and reliability of the satellite measurements. In 
this workshop (Workshop 111), the final two (out of a total of four) months of 
satellite and in situ data chosen for study were processed and evaluated. Re- 
sults from the AVHRR, HIRS, SMMR, and VAS sensors, in comparison with In situ 
data from ships, XBTs, and buoys, confirmed satellite rms accuracies in the 
0.5 to 1.0°C range, but with variable biases. These accuracies may degrade 
under adverse conditions for specific sensors. A variety of color maps, plots, 
and statistical tables are provided for detailed study of the individual sensor 
SST measurements. 



FOREWORD 


This report documents the proceedings of a Satellite-Derived Sea Surface 
Temperature Workshop, held at the Jet Propulsion Laboratory, Pasadena, 
California, on February 22 to 24, 1984. The workshop was the third and final 
of a series designed to compare directly SST measurements from existing 
satellite sensors, thereby evaluating their global accuracies and enabling 
Informed decisions to be made concerning future sensor development. Motivation 
for the comparisons arose from reports of approximately 1°C accuracy in SST 
measurement by four different satellite sensors, and from the need for ocean- 
ographers and climate scientists to understand conditions under which better or 
worse performance could be expected from each sensor type. The workshops were 
convened under the sponsorship of the Oceanic Processes Branch, Office of Space 
Science and Applications of the National Aeronautics and Space Administration. 

A first planning meeting for the workshop series was held at the 
NASA/Coddard Space Flight Center on April 16, 1982, at which the principal 
workshop participants met to discuss approaches to the proposed satellite and 
in situ data comparisons. Subsequent workshop organization took place at JPL 
to facilitate archiving and processing of the satellite and in situ data sets 
on the JPL Pilot Ocean Data System. 

A substantial effort was required to assimilate, process, and analyze 
several months of satellite and in situ data of divergent sampling densities, 
spatial resolutions, and formats. The first workshop, held January 27 to 28, 
1983, was thus limited to a single satellite sensor (SMMR) and one month of 
data (November 1979), and to evaluation of approaches developed to compare the 
satellite data sets with each other and with climatology, ship, XBT, and buoy 
observations. The second workshop, held June 22 to 24, 1983 included data from 
the SMMR, HIRS, and AVHRR, for November 1979 and December 1981. Finally, the 
third workshop was held February 22 to 24, 1984, and included data from the 
SMMR, HIRS, AVHRR, and VAS , for the additional months of March and July 1982. 

An improved set of display products for analysis was developed for this 
workshop. 

- Since the completion of these workshops, increasing attention has been 
paid by the scientific research community to the question of utilizing satel- 
lite SST data in large-scale ocean and atmosphere observing programs. A Tropi- 
cal Oceans and Global Atmosphere (TOGA) Workshop on Sea Surface Temperature and 
Net Surface Radiation was held in La Jolla, California, March 28 to 30, 1984, 
and a Committee for Space Research (COSPAR) International Workshop on 
Satellite-Derived Sea Surface Temperatures for Global Climate Applications was 
held in Washington, D.C. , May 29 to 31 , 1985. Ensuring the long-term avail- 
ability and accuracy of satellite SST measurements is now a high priority for 
ocean and climate research applications. 



Participation In the three NASA/JPL Workshops was varied. Those wno par- 
ticipated In Workshop III and contributed to this report are listed in Appen- 
dix l). i.te ettorts ot all those who contributed to the organization ot the 
worksnops are gratetully acknowledged. Special tnanks are due to Jell Hi l land 
and the statl oi the JPL Pilot Ocean Data System, through whose tireless 
et forts the wine array oi processed SST comparison produ *s was generated and 
uaue available tor the worksnops. 


Lnl G. NJoku 
Workshop Chairman 
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SECTION I 


SUMMARY 
E. Njoku 

Jet Propulsion Laboratory, California Institute of Technology 


In the series of three Satellite-Derived Sea Surface Temperature (SST) 
Workshops, held at the Jet Propulsion Laboratory (JPL) between January 1983 
and February 1934, satellite measurements of SST were reviewed and evaluated 
in detail. The emphasis was on global-scale evaluations, to complement the 
many investigations previously carried out using high-resolution or regional 
data, and to address the SST measurement objectives of large-scale ocean and 
climate programs. 

Four satellite sensors and their associated retrieval techniques were 
reviewed: the Advanced Very High Resolution Radiometer (AVHRR), the High- 

Resolution Infrared Sounder/Microwave Sounding Unit (HIRS/MSU), the Visible- 
Infrared Spin-Scan Radiometer Atmospheric Sounder (VAS), and the Scanning 
Multichannel Microwave Radiometer ( SMMR) . Data from these sensors were com- 
pared with each other and with in situ data from ships, expendable bathythermo- 
graphs ( XBTs) , and drifting buoys. Four months of data were studied: 

November 1979, December 1981, March 1982, and July 1982. 

Principal investigators for each sensor provided JPL with "raw" SST data 
to be processed on the JPL Pilot Ocean Data System (PODS). A variety of dis- 
play products was generated from these data so that results from all sensors 
could be compared in a common format. The display products included color 
maps, histograms, scatterplots, and tables of comparison statistics. Data 
were compared either by point-to-point match-ups (raw data comparisons) or as 
monthly averaged fields on a 2° x 2° latitude-longitude grid (binned data 
comparisons). Workshop Investigators then examined the display products and 
sought to draw conclusions as to the accuracy and error characteristics of the 
sensor SST measurements. 

Conclusions and recommendations arising from Workshops I and II have been 
documented in the reports of those workshops (JPL 1983 and 1984) and in the 
report of a subsequent TOGA Workshop on Sea Surface Temperature and Net Surface 
Radiation (WCP 1984). Conclusions and recommendations from Workshop III have 
been summarized by Njoku (1985) and detailed results from the workshop series 
will be published as a special collection of papers in the November 1985 issue 
of the Journal of Geophysical Research ( JGR Oceans ). This report, therefore, 
serves mainly to provide details of the new data and analysis methods used in 
Workshop III and to supplement discussions appearing in the JGR issue. 

As a very broad and general conclusion, which will be expanded upon in the 
following sections, the workshops have shown that present satellite sensors can 
measure global SST with rms accuracies in the range of 0.5 to 1.0°C. Future 
emphasis must, however, be placed on improved validation and monitoring tech- 
niques to understand the nature of residual spatial and temporal bias varia- 
tions. The present accuracies and geographical distributions of in situ 
sensors are inadequate for this purpose. Fortunately, large-scale ocean/ 
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DATA SETS AND PRODUCTS 
J. Hlllnnd And E. NJoku 

Jet Propulsion Laboratory, California Institute of Technology 


A. BACKGROUND 

Since 1981, the Advanced Very High Resolution Radiometer (AVHRR) on the 
KOAA satellites has made SST measurements in the Infrared portion of the spec- 
trum using a multi-channel technique developed by McClain et al. (1983). In 
order to determine SST more accurately under partially cloudy conditions, 
infrared soundings from the High-Resolution Infrared Sounder (HIRS) and micro- 
wave soundings from the Microwave Sounding Unit (MSU), also on NOAA satellites, 
have been combined in a scheme described by Susskind et al. (1984). Another 
instrument, the Visible-Infrared Spin Scan Atmospheric Sounder (VAS), has pro- 
vided SST retrievals from geostationary orbit (Smith and Woolf 1982). The 
ocean surface has also been viewed in the microwave portion of the spectrum by 
the Scanning Multichannel Microwave Radiometer aboard Nimbus-7. Ullhelt et 
al. (1983) have described the SMMR SST retrieval techniques in detail. 

In contrast to satellite methods, in situ data collected from ships, ex- 
pendable bathythermographs (XBTs) , and moored or drifting buoys have provided 
direct bulk measurements of SST. These platforms have long Berved as ocean- 
ographers' primary cools. Hence, a large body of knowledge has been compiled 
on in situ accuracies, with estimates in the range U.2 to 1.0°C (Saur 1963; 
Tabata 1982). 

Against this background of spacecraft and in situ measurements, partici- 
pants at the JPL workshops sought to review the sensor performances (including 
calibration problems), understand the different SST retrieval algorithms, 
evaluate the sensor SST accuracies, and discuss directions for future sensor 
development. The issue of utilization of satellite SSTs in climate, air-sea 
interaction, and mesoscaic oceanography studies was not the main focus of the 
workshops, but did have a bearing on the recommendations that arose from the 
discussions. Workshop planning and initial results from SMMR were discussed 
in JPL (1983). Workshop II results were more comprehensive as a result of 
refinements in the analysis procedures. In addition, more data were available 
for analysis with the acquisition of AVHRR (MCSST), HIRS/MSU, and VAS data 
sets as described in JPL (1984), In Workshop III, the subject of this report, 
the acquisition and analysis of data were completed, recommendations for future 
research and sensor development were discussed, and plans were made for event- 
ual publication of results from the workshop series in the open literature. 


B. DATA SET CHARAC1KRISTICS 

Each sensor collected data in a unique manner, due to resolution and scan- 
ning methods, as well as a result of sensor processing and duty cycles (the 
percentage of time the sensor provided data satisfactorily relative to the 
total time). Therefore, the data distribution varied greatly. Table 2-1 sum- 
marizes pertinent sampling parameters for each sensor. The comparison shows 
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Table 2-1. Sensors Evaluated During the SST Workshops and 
Resolution .if Derived SSTs 


Platform 

Sensor 

Duty Cycle 
(X) 

Spatial Resolution 
of Derived SST (km) 

Coverage 

TIROS-N) 

(AVHRR 

100 

25 

Global 

NOAA-7 ) 

(HIRS/MSU 

100 

125 

Global 

Nimbus-7 

SMMR 

30 

150 

Global 

GOES 

VAS 

75 

50 

N.W. Atlantic 
S.E. Pacific 

Ships 

Thermometer 



Clobal (mostly 
northern hemisphere) 

XBT 

Thermistor 



N. Pacific 

FCCE 

Thermistor 



Southern hemisphere 


that the high resolution, nearly continuous duty cycle of the visible-infrared 
sensors yields an enormous number of discrete radiances. Spacecraft-treasured 
radiances were averaged as part of the conversion process to geophysically 
meaningful temperatures before delivery to JPL, thus somewhat reducing the data 
volume. Details of the sensor modes of operation were reported in JPL (1983 
and 1984), but will be summarized below for reference. 

1 . AVHRR 

Global, duy/night coverage across a 2,500 km wide swath at 4 km 
resolution characterizes the fundamental global sampling of AVHRR instruments 
aboard the NOAA satellites. Prior to mid-November 1981 the Clobal Operational 
Sea Surface Temperature Computation (GOSSTCOHP) provided 50 km resolution SST 
retrievals using a single window (centered at 11 urn) algorithm. After this 
date che improved five-channel ins': rumen t was used to derive SST from the 3.7, 
11, and 12 Pm windows, utilizing the triple-window technique known as the 
multi-channel sea surface temperature (HCSST) algorithm. A spatial resolution 
of 25 km was retained. The standard NOAA GOSSTCOHP and HCSST products 
provided geolocatod SSTs and supporting parameters such as platform source, 
data quality, and day/night status. 

2 . H IRS/HSU 

The HIRS and HSU instruments flown on TIROS-N and NOAA-7 served as 
sources for derived SSTs. The large number of Infrared and microwave channels 
are combined in a physical algorithm to produce surface temperatures under 
clear or cloudy conditions. Arrays of HIRS soundings (instantaneous iield of 
view 17.4 km at nadir) are averaged across the 2,300 km wide swath to form 
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SSTs at a spatial resolution of 125 ka. Because the retrievals are spatially 
averaged, points within 60 ka of land tend to be contaminated by varner land 
temperatures. For workshop purposes, opace/tlee location and quality param- 
eters were provided so that Information could be segregated for various study 
oonths, regions, and sampling conditions such as day/nlght and clear/cloudy. 

3. SfWR 

Dual polarized microwave radiance measurements at 6,6, 10.7, 18, 
and 21 CHz were the fundamental Input to SKMR SST algorithms. The SMMR samples 
along a nadir-centered 780 km wide swath with a spatial resolution of 150 km at 
6.6 GHz. Various quality control criteria were applied to the data by the 
Nimbus-7 algorithm development team, but the most obvious and influential, with 
regard to workshop processing, was a land proximity mask. All data within 600 
km of land. Including large islands, were eliminated due to possible antenna 
sldelobe contamination. Furthermore, In order to distinguish the highest 
quality values, data sent to JPL were flagged for day/ twilight/night status 
and cell (1-5) location In the swath. The instrument was turned on and off 
every other day due to spacecraft/power limitations. Furthermore, the end 
cells In the swath were deemed unrcliabile due to polarization correction 
errors. Thus, the overall duty cycle was reduced to 302. 

4. VAS 

The Geostationary Operational Earth Satellite (GOES), carrying VAS, 
provided a stable platform for scanning the full disc. Daytime-only IR and 
VIS data were collected as part of the normal operations. These data were 
screened for cloud-free areas as part of the SST derivation scheme. Three of 
the twelve thermal bands sensed by VAS were used to derive SSTs at a spatial 
resolution of ~50 km. Finally, retrievals from the eastern tropical Pacific 
and northwestern Atlantic were provided for evaluation. 

5. In Situ 

The primary surface data set consisted of ship intake temperature 
measurements collected from radio reports by the Fleet Numerical Oceanography 
Center (FNOC). Typically, intake temperatures are accurate to the nearest 
1°C. Additionally, biases on the order of tenths of a degree Celsius have 
been reported (JPL 1983). However, these data are the sole source of global 
in-situ measurements and at best provide spotty spatial coverage in the 
southern hemisphere. The temporal resolution of most reports is six hours. 
Marine reports were closely scrutinized for pathological errors related to 
erroneous 6hip locations and extreme temperatures. 

Complementary data sets consisting of XBT drops across the northern and 
tropical Pacific and measurements from drifting buoys in the southern Pacific 
provided additional validation information. XBT drops between North America, 
Japan, Hawaii, Tahiti, and the Panama Canal are made about every 200 km, 
yielding 400 to 1 ,000 observations during any month. Reported accuracies are 
0.1°C with biases of about +0.2°C relative to salini ty-tenperature-depth 
(STD) instruments. 
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Drifting buoyo launched during the First CARP Clobal Experiment (FGCE) 
supplemented ship observations in the southern hemisphere from 20 S to 65 S. 
Hore than 100 buoys reported SSTs at roughly 6-hr Intervals throughout each 
study month. In an examination of the buoy program (Carrett 1981), comparisons 
made with ship measurements within 1 hour and 100 lua of the buoy observation 
yielded a worst-case standard deviation of 1 .48°C and an average bias of 
+0.28°C relative to the ships. Allowing for buoy temperature sensor stabili- 
zation reduced the standard deviation and bias to 1.15°C and -r0.75°C when 
compared to Intake temperatures and 0.56°C and +0.18°C, respectively, when 
compared to bucket measurements. It should be noted that these statistics 
were determined for buoy SSTs taken 24 hours after the ship recording. 

C. ANALYSIS TECHNIQUES 

Very different spatial and temporal sampling characteristics, as well as 
a low slgnal-co-noise ratio, characterized the data. Hence, analysis 
procedures were designed to reduce the noise by forming monthly, 2° latitude 
by 2° longitude average SST anomaly fields. Noise levels were determined 
from point-to-point or "spot" comparisons of SST anomalies. An anomaly is 
defined as the departure of absolute SST from climatology. Anomalies were 
computed by linearly interpolating, in space and time, a 1° by 1° clima- 
tology generated by Reynolds (1982) to the irregularly spaced satellite or 
surface point and subtracting the climatology from the measured temperature, T. 
The resultant value, hereafter referred to as a "raw" anomaly, AT, was used 
as the fundamental signal rather than the absolute SST. Thus, a picture of 
ocean variability was depicted by each sensor. 

A variety of statistical and display routines, summarized in Table 2-2, 
was used to portray raw-anomaly quantitative results and spatial features. 

First and second moments were computed in the usual Gaussian sense: 






( 2 - 1 ) 



and 


AT ^ 
rmsd 


[1/n ^2 < AT i ~ AT) 2 J 1/2 

1-1 


( 2 - 2 ) 


where n is the total number of points, ATj is the raw anomaly located at 
latitude y and longitude x at time t, and AT rnsc j is the root-mean-square 
deviation about the mean. 
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Monthly average fields were formed from raw 
points that fell within a 2° by 2° cell centered 
longitudes* The average temperature or "binned" 
longitude, 1, is simply expressed as 


ra 

k-1 


where m is the total number of points in cell j , 
anomaly located at latitude y and longitude x within the cell centered at J , 1* 

Statistics of the binning process were retained for the purpose of compar- 
ing sampling charateristics , temperature extremes, and data dispersion within 
each cell. Raw anomalies exceeding +5.75°C were eliminated before binned 
anomaly fields were formed, because the natural variability of the ocean is 
typically much less than this magnitude. It follows that any signal of this 
intensity is the result of poor sensor performance or algorithm deficiencies, 
except perhaps in the case of a strong El Nirfo. No further quantitative edit- 
ing was performed on SST anomalies. However, the data were stratified into 
latitude/longitude bands and segregated based on a qualitative interpretation 
of the status flag associated with each observation. Table 2-2 summarizes the 
screening procedures applied to each data set. Field-data statistics were 
calculated as in Eqs. (2-1) and (2-2) by substituting the mean anomalies for 
the raw values. 

Sensor and algorithm performance were measured relative to climatology and 
to each sensor. Statistics of the relationships quantified the bias, standard 
deviation, and correlation. The correlation between any two sensors (clima- 
tology was treated as a sensor) is given by: 


anomalies by averaging all 
on odd latitudes and 
anomaly for latitude, j, and 


(2-3) 


i , and T yx is the raw 


N _ _ . N _ N 

2 T i„ t 2„ - » 2 £ * 2 „ 

n-1 n»l m-1 



where subscripts l and 2 refer to sensor pairs and subscripts n and m 
the cells common to both sensors. 


(2-4) 


denote 
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Table 2-2 • Analysis and Graphical Display Techniques Applied Co 
Satellite-Derived SSTs 


Techniques 


Data Type Analyzed 

Raw Raw Binned 

Absolute Anomaly Anomaly 


Compute raw anomalies by interpolating 
climatology to satellite point and subtracting 

Form 2° latitude x 2° longitude monthly 
averages (binning) 

Calculate absolute SST by adding binned 
ancmalies to 2° binned climatological SST 

Prepare a histogram of SST and summary 
statistics; mean, standard deviation, and 
number of observations 

Prepare a contour map of binned absolute SST 

Prepare a thematic map of binned anomalies 

Prepare a thematic map of anomaly differences 

Draw a thematic map of number density within 
a 2° cell 

Make a scatter diagram of binned anomaly SST 
versus Reynolds climatological SST: summary 

statistics, bias, standard deviation about 
bias, correlation and number of observations 

Prepare cross correlation tables: statistics, 

correlation, bias, standard deviation about 
bias, and number of 2° cells 


Calculate error partitioning tables: statistics 

and overall ras error contributed by each sensor 
and rms error for each sensor combined with two 
other sensors 

Make a scatter diagram of sensor versus 
sensor: statistics, correlation, bias, 

standard deviation about bias, and number of 
2° cells 

Hake a scatter diagram of binned anomaly 
differences for a sensor pair versus number 
of observations in 2° ceils for either sensor 
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Another information extraction technique, error partitioning, was used to 
ascribe a measurement error to each sensor* This method employs sensor trip- 
lets in a set of three simultaneous equations which can be solved for each 
sensor's contribution to the total error in an rms sense* The mean-square 
difference between binned anomalies from two sensors is expressed as: 


°12 - «T, - V 2 > - I/n]T <7 1» - 


n»l 


(2-5) 


where < > indicates the sample mean computed over all binned points, n, 
that are common to the sensor triplet, and Tj and T 2 are the binned SST 
anomalies* The error in the measured anomaly is 



(2-6) 


where T^ is the measurement from sensor k and T t is the true anomaly. 
Then, it is possible to express Eq. (2-5) as: 


°12 " <(e l ~ e 2 )2> “ <E 1 2> + <G 2 2> 


(2-7) 


It has been assumed that sensor errors are uncorrelated; hence, the mean cross- 
product of the errors is zero. Similar expressions may be derived for D 13 
and D 23 , thus forming a set of three equations that can be solved for the 
sensor errors, <e£>, k ** 1,2,3 [JPL (1983), Appendix Gj. 

In this manner, error estimates for sensors forming the triplet can be de- 
termined. For M sensors there will be (M^ - 3M + 2)12 possible triplet com- 
binations containing a given sensor. An overall error estimate can be obtained 
for each sensor by averaging the partitioned error for that sensor in each 
triplet combination. 

These analytical techniques were applied to the binned anomaly fields and 
raw anomalies. Results were displayed from field data on monthly global and 
regional scales and from raw data in the form of "spot" comparisons on spatial 
and temporal scales commensurate with sensor sampling and geophysical variabil- 
ity. The space-time scales used for preparing analysis products are presented 
in Table 2-3. 



Table 2-3. Data Analysis Products Spatial and Temporal Length Scales^ 


Analysis Product Global^'^ 

North 

Pacific 3 

South 

Pacific^ 

North 

Atlantic 

Monthly 
3 Average 

Anomaly histogram 

X 

X 

X 

X 

X 

Absolute SST contour map 

X 

X 

X 

X 

X 

Thematic map of binned 
anomalies 

X 

X 

X 

X 

X 

Thematic map of anomaly 
differences 

X 




X 

Thematic map of number 
density 

X 




X 

Scatter diagram of anomaly 
SST versus climatological 
SST 

X 

X 

X 

X 

X 

Cross correlation table 

X 

X 

X 

X 

X 

Error partitioning table 

X 

X 

X 

X 

X 

Scatter diagram sensor 
versus sensor 


X 

X 

X 

Spot merge 
+6, +12 hr; 
20, TOO km 

Scatter diagram anomaly 
difference versus numerical 
observ. m 2° cell 


X 

X 

X 



1. For some products the Pacific study area was separated into 3 regions of 

latitude to separate the tropics and extratropics. These regions were lati- 
tudes: (1) below 20°S, (2) between 20°S and 20°N, and (3) above 20°N. 


2. Global study area; 60° S to 60° N, 0° to 360° E. 

3. North Pacific study area: 0° to 60° N, 100° to 290° E. 

4. South Pacific study area: 60° S to 0°, 100° to 290° E. 

5. North Atlantic study area: 0° to 60° N, 290° to 360° E. 

6. Global and regional thematic maps and scatter diagrams within 20° 
latitude bands extended to 60° latitude. All other products terminated 
at 55° latitude to eliminate possibly spurious points due to sea ice. 
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Spot comparisons of the raw data were made using a number of different 
space-time "windows" in which sensor data points were matched within a given 
time tolerance (+ hr) and radial distance (km). Table 2-4 provides the window 
specifications for matching any given sensor to the sensors listed in the 
Table. 

D. WORKSHOP III RESULTS 

Detailed analysis results from Workshop III have been compiled into a 
special collection of papers to be published in the November 1985 issue of the 
Journal of Ceophyslcal Research ( JGR Oceans) . In this report, some additional 
sets of data products are provided. Tables of correlations and other statis- 
tics can be found in ?pendix A. Error partitioning results are provided in 
Appendix B, and selected color images of SST anomaly fields are shown in 
Appendix C. A lengthy set of histograms and scatterplots were also generated 
for the Workshop and can be individually copied and supplied by request to the 
authors. Individual summaries of findings from the Workshop data are provided 
in Sections III to VIII. 

1. Data Sets 

Table 2-5 shows the satellite and in situ data sets processed dur- 
ing the three JPL workshops. In Workshop III, data from all sensors for the 
months of March and July 1982, and also HIRS data for December 1981, were 
evaluated. Details of the procedures used to derive the SST data can be found 
in Sections III to VI, or in previous workshop reports (JPL 1983 and 1984). 
Prior to Workshop III some changes were made in the HIRS SST algorithm to 
improve its performance. These changes included a higher-resolution climatol- 
ogy to define land points, tightened criteria for internal consistency check- 
ing, and retention of more data samples in each cloud analysis area (see Sec- 
tion IV). This new version of HIRS data was named "Version 2." In addition, 
both HIRS versions were provided with data quality weights which could be 
applied in the binning process. Those SST data with weights applied were 
referred to as HIRS "weighted." 

Another data set unique to Workshop III wa3 provided by T. Wilheit. This 
"SMMR/ship" data set consisted of a blended monthly averaged SST field on a 
2° x 2° grid, for the Pacific Ocean region, comprising original raw SST 
data from the SMMR and from FNOC ships. The mechanism for blending the SMMR 
and sitip data is outlined in Section VI and, in essence, uses the ship data, 
where available, in an objective analysis scheme to remove spatial biases in 
the SMMR SST field. 

2. Analysis Products 

A great improvement in Workshop III was the availability of color 
images, or "thematic maps," to display global SST anoraa.’ • fields and data 
distributions. Table 2-6 lists the complete set of global color maps that 
were produced, Including retrospective processing of maps for November 1979 
and December 1981. Those maps, which have beer, referenced by discussions in 
Sections III through VI, are shown in Appendix C. The maps are self- 
explanatory. 
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Table 2-4. Tine-Space Windows Used for Merging Sensor Data to Data from 
Other SensorB Listed in the Table 


Sensor 

Time 

<+hr) 

Radial Distance 
(km) 

AVHRR 

12 

100 

Ships 

6 

100 

XBT 

12 

20, 100 

FGGE Buoys 

12 

20 


Table 2-5. 

Data Sets 

Available for Workshop Processing 


November 

1979 

December 

1981 

March 

1982 

July 

1982 

SMMR 

1 

2 

3 

3 

AVHRR 

2 

2 

3 

3 

HIRS/MSU 

2 

3 

3 

3 

VAS 

- 

- 

3 

3 

FN0C Ship 

1 

2 

3 

3 

FGGE Buoy 

2 

- 

- 

- 

TRANS PAC XBT 

2 

2 

3 

3 


Key: 1 ** available for Workshop I 

2 - available for Workshop II 

3 - available for Workshop III 
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Table 2-6. SST Workshop Color Maps (Global) 



Nov 79 

Dec 81 

Mar 82 

Jul 82 

Data Distribution: 
AVHRR 

/ 

/ 

✓ 

✓ 

H1RS 

✓ 

✓ 

✓ 

/ 

SMMR 

✓ 

✓ 


✓ 

VAS 

0 

0 

/ 

/ 

Ship 

/ 

/ 

✓ 

/ 

XBT 

✓ 

✓ 

✓ 

/ 

FGGE Buoy 

✓ 

0 

0 

0 

SST Anomalies: 

AVHRR 

/ 

✓ 

✓ 

✓ 

H1RS 

/ 

✓ 

✓ 

/ 

SMMR 

✓ 

✓ 

/ 

/ 

SMMR/Ship 

0 

0 

✓ 

✓ 

VAS 

0 

0 

✓ 

✓ 

Ship 

/ 

✓ 

✓ 

/ 

XBT 

✓ 

✓ 

/ 

✓ 

FGGE Buoy 

/ 

0 

0 

0 

SST Anomaly Differences: 

AVHRR (day-night) 

0 

✓ 

/ 

/ 

HIRS-AVHRR 

/ 

J 

/ 

✓ 

SMMR- AVHRR 

/ 

✓ 

✓ 

✓ 

SMMR/Ship-AVHRR 

0 

0 

✓ 

✓ 

VAS-AVHRR 

0 

0 

/ 

✓ 

HIRS-SMMR 

✓ 

/ 

/ 

✓ 

AVHRR-Ship 

✓ 

✓ 

✓ 

/ 

HIRS-Ship 

✓ 

✓ 

/ 

/ 

SMMR-Ship 

✓ 

✓ 

/ 

✓ 

SMMR/ Ship-Ship 

0 

0 

✓ 

/ 

VAS-Ship 

0 

0 

✓ 

✓ 

XBT-Ship 

/ 

✓ 

✓ 

✓ 

XBT-AVHRR 

/ 

✓ 

/ 

/ 


Key: / Available 

0 Not available 
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The Cables in Appendix A provide staciscics (as discussed in Section II. C) 
fpr pair-wise sensor comparisons of binned data for all months. All data have 
been masked up to 600 ka from land to conform to the SMMR areal coverage. 
Another benefit of the masking is that high-gradient coastal current regions, 
which could cause sampling errors in the 2°-binned analyses, are largely 
eliminated. For each month, there is one global set of statistics and another 
set for each of four ocean regions. This allows the contrasting of statistics 
for tropics and extratropics, or for North Atlantic and North Pacific. A 
second table (A-2) was generated using the same data as for Table A-l, but 
with a 3 x 3 cell two-dimensional spatial filter applied. The weights for the 
filter were as follows (normalized by a factor of 16): 


1 2 1 

2 4 2 

1 2 1 


The motivation for applying the filter was to investigate any reduction in the 
standard deviation or increase in the correlation therefrom. (Climate-scale 
models can accommodate data on 2° x 5° or even 5° x 5° spatial scales.) 
However, care must be taken in interpreting the statistics of Table A-2 since 
the 2°-binned sensor data are no longer strictly independent after the 
spatial filter has been applied. Finally, Table A-3 has been included to show 
the effects of using the HIRS algorithm quality weights on the March 1982 data 
(see Section IV). 

The error-partitioning tables in Appendix B follow from the discussion in 
Section IIC. Only the four data sets with truly global coverage have been 
included. Results could be significantly skewed by the addition of ship data, 
unless all data were restricted to the North Atlantic or North Pacific. 
Separate tables for these regions were not produced, however, at the 
Workshop. As in Appendix A, separate tables were generated for unsmoothed and 
smoothed data sets. The HIRS data used was Version 2. For each sensor, the 
partitioned error from each of three possible triplets (or triads) was 
averaged to give an overall average rms error. These average errors have been 
collected together for convenience in Table B-3. 

E. DISCUSSION 

The data products (maps, plots, statistics) generated for the Workshops 
led to very detailed discussions concerning sensor calibration, algorithm 
performance, error characteristics, error sources, validation problems, and 
future research. Participants were encouraged to submit summaries of their 
evaluations for the Workshop report. These investigator summaries are 
provided in Sections III through VIII. (Some investigators chose to postpone 
publication of their analyses until the JGR Oceans special issue.) 
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SECTION III 

ADVANCED VERY HIGH RESOLUTION RADIOMETER (AVHRR) 
SENSOR COMPARISON CHARACTERISTICS 

E. Paul McClain 

National Oceanic and Atmospheric Administration 


A. INTRODUCTION j 

Multi-channel sea surface temperature (MCSST) techniques make use of the 
AVHRR on the NOAA satellites (Schwalb 1978), which delivers i.l-km and 4-km 
resolution measurements in four (or five) channels: 0.58-0.68, 0.725-1.10, 

3.55-3.93, and 10.3-11.3 um (also 11.5-12.5 pm for five-channel AVHRRs). 

Various combinations of the four or five AVHRR data channels are used in the 
important initial cloud-filtering stages. Thresholds of the bi-directional 
reflectance in the visual and reflected-IR channels, and of brightness 
temperature in the thermal-IR channels, have been established for cloud-free 
conditions. Homogeneity tests take advantage of the low radiometer noise 
levels and the high degree of spatial uniformity in the ocean surface reflec- 
tance and temperature fields in the absence of clouds. Details of the cloud 
tests are given in McClain et al. 1983, and in Part II, pp. 1-8, in the SST 
Workshop II Report (JPL 1984). 

The brightness temperatures measured in given window channel of the ; 

AVHRR are corrected for atmospheric attenuation by use of the brightness tern- j 

peratures from two or all three of the atmospheric windows, each of these I 

6pectral bands being characterized by a different atmospheric transmittance. K 

The relationships between the atmospheric correction and various combinations 
of brightness temperature differences is linear with exceptionally small scat- j 

ter, as was determined from simulation data bases (McClain 198 1 ) . Satellite 
and buoy measurements matched to within 25 km and 24 hours have been used to ■ i 

i $ 

derive a small but significant temperature-dependent bias correction term for 
the simulation equations. Further details of deriving multichannel sea 

surface temperatures from AVHRR measurements are given in McClain et al. 1983. : 

B. ADVANTAGES OF THE AVHRR (MCSST) TECHNIQUE 

The principal advantages of the AVHRR-derived sea surface temperatures 
are high resolution, broad geographic coverage, general consistency, and good 
accuracy. The MCSSTs have been produced on an operational basis since Novem- 
ber, 1981, and they are available globally from N0AA/NESDIS as a monthly mean 
(65N-65S, 2.5° grid) or as a weekly composite (70N-70S, 100-km grid), and as 
selected regional (50-km grid) or local charts (14-km grid). The basic re- 
trievals are obtained from 2x2 arrays of AVHRR data over an area nominally 
8 km on a side, then one or more of the retrievals per 25-km box are resolved 
onto the various grid intervals listed above. MCSSTs can be obtained within 
10-20 km of a land, ice, or cloud edge; and except in regions of extremely 
persistent and continuous cloud cover, most areas of the world are sampled at 
least once on a 100-km grid every 5-10 days. 
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In comparison with in situ sensors, particularly those on ships of 
opportunity, the AVHRR with its nearly continuous onboard calibration 
generates an internally uniform set of brightness temperature measurements 
orbit after orbit, day after day. Except for quite unusual circumstances, 
such as the El Chlchon eruptions and the electrical Interference problems that 
have plagued the 3.7 ut data during certain periods (discussed in next 
section), the operationally derived MCSSTs generally comprise a spatially and 
temporally consistent data base. The few changes that have been made in the 
operational algorithms have affected the root mean square (rms) differences 
with respect to drifting buoy temperatures only at the <0.25C level. Recent 
drifting buoy spot comparlsions over a wide range of temperatures, geographic 
area, and seasons consistently indicate biases of < 0 . 1 C and rms diifcrcnces 
(or scatters) of 0.5-0.6C (Strong and McClain 1984). Comparisons with 
screened ship observations, after removal of a common ship-based climatology 
to derive anomalies, are summarized for the various periods and regions in 
Table 3-1 along with the statistics for the other sensors studied during the 
several JPL SST Workshops. It Is evident that the AVHRR, compared with the 
other sensors, was almost always and everywhere oer the globe characterized by 
the largest number of matchups, the lowest bl3S and scatter, and the highest 
correlation relative to the Pazan set of screened ship observations (the 
prnclpal exception to this occurred In July 1982 in connection with the El 
Chlchon eruption effects). 


Table 3-1 gives the various statistics for each sensor with respect to 
ship matchups more than 600 km from any land or ice surface. The use of a 600 
km mask is necessary in order to equalize the coverage for the AVHRR, 

H1RS/MSU, and VAS , which can obtain observations near coastlines and ice 
edges, with that available for the SMMR, which is constrained to operate at 
least 600 km away. Furthermore, it should be noted that there are no 
measurements from the SMMR over the North Atlantic in November 1979, and the 
VAS coverage is limited to two areas in March and July of 1982, one about 
25° latitude by 30° longitude on a side in the southwestern North Atlantic 
and the other about 40° latitude by 40° longitude in the extreme eastern 
equatorial and northeastern South Pacific. 


Table 3-1 also enables comparison of sensor statistics for ship matchups 
>600 km from land or ice and those for the same matchups after a special 3x3 
weighted smoother^ is applied. This procedure was used on all available 
2° latitude-longitude bins >600 km from land or ice. A sharp drop in sample 
size resulted from loss of outer rows and columns of bins and from no computa- 
tion being made for arrays where inadequate in aitu or satellite data resulted 
in no bln average being computed for one or more bins of the 3x3 array. This 
was a particularly acute problem in the case of the FGGE buoy set; c.g., the 
unsmoothed global data set of N«400 for the AVHRR matchups was reduced to N**l 
for the smoothed set. 


*A smoothing set of weights is applied to the 
monthly mean SST anomalies (Tj, T 2 , ...» T 9 ) 
in overlapping 3x3 arrays of the two-degree bins. 


Ti T 2 T 3 
T 4 T 5 T 6 

t 7 t 8 t 9 


t ' 5 - 4*T 5 + 2*(T 2 + T a + T 6 + Tg) + T t + T 3 + T 7 + T 9 ) /16 
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Table 3-1. Sensor Sea Surface Temperatures Compared With Pazan Screened Ship 
Temperatures (>5 per bln). Values In Parentheses Are After a 3x3 
Center-Weighted Smoother Was Applied (A - AVHHR; S - SMMR; V - VAS, 
H - HIRS/MSU). 
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The effect of the smoother on the bias is mixed at best, but there is 
definite (often dramatic) improvement in the scatter (i.c. standard deviation) 
and cross correlation statistics in nearly every case for every sensor. This 
is noteworthy in particular, for it demonstrates Just how important a factor 
is the uneven quality of the in situ temperature observations when one is 
using them in an attempt to validate satellite-derived temperatures. 

Table 3-2 displays the statistics for AVHRR/ship matchups >600 km from 
land or ice with various types of in-sltu data (viz., Pazan ships, >5 and >20 
per cell; FGGE buoys; and Transpac XBTs); and Table 3-3 shows AVHRR/ship 
matchup statistics separated into day and night. Table 3-4 enables comparison 
of two types of in situ observations, viz., Pazan ships the Transpac XBTs. 

C. LIMITATIONS OF AVHRR (MCSST) TECHNIQUE 


Perhaps chief among the limitations of the MCSST, or any other infrared 
method, is lack of retrievals in areas of persistent and essentially continuous 
cloud cover. The relatively high resolution of the AVHRR does enable more 
retrievals to be made in areas of patchy cloud cover than can be done wth the 
other sensors. The cloud detection tests appear to work well in that retrievals 
are seldom made using cloud-contaminated data, including contamination by thin 
cirrus or sub-resolution cumulus fields. 


Severe volcanic eruptions or dust storms can produce extraordinary aerosol 
loadings in the atmosphere and thereby greatly increase attenuation of the in- 
frared signal reaching the satellite, as well as interfere with those cloud de- 
tection tests that depend on visible-band measurements. El Chlchon, because of 
the large mass of H 2 SO 4 droplets found at very high altitudes in the atmo- 
sphere, was particulary severe in its impact on the MCSSTs, especially in the 
Northern Hemisphere tropics and subtropics. Daytime retrievals were virtually 
eliminated between 5N-30N for up to six months, and nighttime retrieval!, were 
biased too low by up to 2-4C from April to October 1982 (Strong et al. 1983). 
Recent research Indicates the very real possibility of using a different 
formulation of the triple-window MCSST equation, one that appears to be nearly 
insensitive to the concentration of El Chichon type aerosols (Walton 1985). 
Furthermore, daytime visible band data from the AVHRR can probably be used to 
obtain a point-to-point measure of the aerosol loading, thus leading to other 
possibilities for correction of the retrieval temperatures (Stowe 1984). 


Use of any triple-window equation obviously needs noise-f ee measurements 
in all three IR-window channels. The noise level has been exceptionally low 
<0.1 K, in the 11 and 12um channels. The 3.7 pm window data, however, tend to 
have an acceptable noise figure (< 0.2 K) during the firBt 12 months or so after 
each satellite launch, but then become increasingly contaminated by electrical 
interference thereafter. Fortunately, so-called "outgasslng" procedures that 
were implemented recently 2 successfully reduced the complex but coherent noise 
in the 3.7 pm data to levels comparable to chose measured immediately after 
1/unch. More recent experience indicates that outgassing must probably be re- 
peated every 6-12 months to control this problem. Although the noise level in 
this channel was moderately high In July 1982, this was reflected primarily in 
a loss of nighttime MCSSTs from failure of the MCSST Intercomparison Test 
(see McClain In Part II, pp. 1-8, JPL 1984) rather than a decrease in accuracy. 
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2 0n 6/15/83 for NOAA- 6 , on 9/27/83 for N0AA-7 , and on 9/8/83 for N0AA-8. 
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Table 3-2. AVHRR Multichannel Sea Surface Temperatures Compared With In Situ 

Temperatures. Values In Parentheses Are After a 3x3 Center-Weighted 
Smoother Was Applied (P5 » Pazan Screened Ships. >5/Ein; P20 - Pazan 

Screened Ships, >2Q/Bin; F ■ FGGE Buoys; T ■ Transpac XBT's). 
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Table 3-3. AVHRR Multichannel Sea Surface Temperatures 

Compared With Pazcn Screened Ship Temperatures 
(»5/Bin) by Day and by Night. Values in 
Parentheses Are After a 3x3 Center-Weighted 
Smoother Was Applied. 
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Table 3-4. Pazan Screened Ship Temperatures (>5/Bin) Compared With 

Transpac XBTo and AVHRR Multichannel Sea Surface Temperatures. 
Values In Parentheses Are After a 3x3 Center-Weighted Smoother 
Was Applied (A - AVURR; T - Transpac XBTs). 
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Although the brightness temperature at 3.7 pm can be significantly ele- 
vated above the emission value by even small amounts of specularly reflected 
solar radiation (>1.7Z normalized bi-directional reflectance) and present 
operational MCSST procedures use only split-window, i.e., 11 and 12 pm 
brightness temperatures in the daytime for this reason, recent tests indicate 
that no detectable solar contamination exists at satellite zenith angler 
greater than 3-10° on the anti-solar side of the satellite subpoint track. 

Thus dual-window (3.7 and 11 pm) or triple-window MCSST equations can be 
used in the daytime over large areas of the globe. 

Concern has been raised over use of the standard MCSST equations, particu- 
larly the split-window one, when steep temperature inversions are present in 
the atmosphere just over a water surface (see Condal et al. in Part II, pp. 
29-31, JPL 1984). Such conditions are prevalent at times over the Great Lakes 
in late spring and early summer, and can occur in any mid-latitude coastal 
waters during this season if offshore winds are blowing. Recent investigations 
of this with buoy data from the Great Lakes in 1982 and 1983 found that the 
operational split-window MCSST equation actually performed rather well except 
under the most extreme inversion conditions. No reliable methods for satellite 
detection of the presence of these extreme conditions, and correcting for their 
effects on the temperature retrieval, have been found to date. 

The old problem remains of skin-versus-bulk temperature, as satellite IR 
techniques yield skin temperatures at depths of less than a millimeter and in 
8 itu methods of observation give bulk temperatures at depths ranging from a 
few centimeters (towed thermistors), to a meter or two (buoys), to three to ten 
meters (ship Intakes). The use in the MCSST metho l of a temperature-dependent 
bias correction derived from satellite/buoy matchup data presumably incorpor- 
ates some sort of average skin-vs-l m depth temperature adjustment, but this 
effect is almost always negative (i.e., the skin is cool) and generally amounts 
to 0.1-0.5C in magnitude (Robinson et al. 1984). It is not uncommon for the 
top few tens of cm of the water surface to become heated under low amounts of 
cloudiness and when very light winds result in little mechanical stirring; 
this has been termed the "diurnal thermocline" (RobinBon et al. 1984), 

D. ABILITY TO SATISFY USER NEEDS 

Obviously the higher the temporal and spatial resolution, and accuracy, 
and the more complete the coverage under all meteorological conditions, the 
more user needs will be satisfied by the satellite-derived SST measurements. 
Lesser resolution data in time and space presumably can always be derived from 
the original observations. With its capability to produce relatively high 
resolution daily to weekly observations very close to coastlines or ice edges, 
the AVHRR can satisfy many users, both oceanographic (including fisheries) and 
meteorological, although extensive clcud cover may be constraining in some 
areas and times. Good global MCSST coverage for the relatively low-resolution 
monthly mean charts is virtually always available for climatological users. 
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E. ERROR CHARACTERISTICS OF MCSST'S FROM STATISTICAL TAFiLES 

In addition to the various statistics summarized in Table 3-1 , which all 
refer to the Pazan set of screened ship data, and where all anomalies were de- 
rived using the Reynolds' climatology (Reynolds 1982), numerous ship and satel- 
lite anomaly fields, and associated difference fields, were also produced in 
connection with the workshops. The ship anomaly fields suffer from the tradi- 
tional lack of observational coverage in the Southern Hemisphere and some other 
areas. The Transpac XBT fields are severely constrained by limited geograph- 
ical coverage. With this in mind, a discussion of each data period follows. 

1. November 1979 Data Period 

This first data period is different from the other three in that no 
MCSSTs were yet available, the NOAA operational product at that time being an 
improved AVHRR/HIRS version of the earlier GOSSTCOMP (^Global Operational SST 
Comp utation) , which was based on data from the SR/VTPR instruments on the 
pre-TIROS-N generation of NOAA operational polar satellites (Walton 1980). 
Compared with the SR/VTPR data, the AVHRR/HIRS measurements are of better 
quality and higher spatial resolution. Consequently, although still based on 
a single AVHRR window channel, the satellite retrievals based on the GOSSTCOMP 
procedures from 1979 through mid November of 1981 were generally of higher 
quality than those in 1978 and earlier. 

Table 3-1 indicates that AVHRR biases are comparable or a bit larger than 
those for HIRS/MSU, but significantly smaller than those for SMMR. AVHRR 
scatter and correlation values are much better than for HIRS/MSU and very much 
better than for SMMR. 

Table 3-1 also demonstrates the very substantial improvement in scatter 
and cross correlation figures for all three sensors and in nearly every area 
that results from the 3x3 weighted smoothing procedure (the mid-Pacific 
values are of dubious reliability because of the small sample): tu? AVHRR 

statistics are especially impressive (scatter of 0.33-0.38C and cross 
correlation of 0.41-0.85), whether on an absolute basis or compared with the 
other sensors. Biases generally worsen somewhat for all sensors when the 
3x3 smoother is used. 

In Table 3-2 (AVHRR only) comparisons by in situ data type are limited by 
sampling constraints, especially after the 3x3 smoother is applied, so not 
all regions are represented. There is very little difference in bias from one 
type of in-situ data to another, but scatter values are clearly worse relative 
to FGGE buoys 3 and a bit worse relative to Transpac XBTs. Cro3s correlation 
comparison are mixed. Scatter and correlation figures generally improve sub- 
stantially when the 3x3 smoother is used. 


3 This was unexpected as previous studies have found lower bias and scatter 
relative to buoys than to ships (Strong and McClain 1984). An independent 
study of the FGGE buoy data set for Nov. 1979 for another purpose by a NOAA/ 
NESDIS contractor found that some of the FGGE buoy observations were seriously 
in error. This is significant because neither the buoy3 nor the XBTs were 
given the same kind of extensive screening as were the Pazan ships. 
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2. December 1981 Data Period 

This was the first full month of operational MCSST processing; and 
as in the Nov. 79 period • Table 3-1 indicates that AVHRR biases were comparable 
to HIRS/MSU ones and significantly better than those associated with SMMR, and 
the AVHRR biases tend to be better with the 600 km mask in place. Likewise, 
scatter magnitudes are much better for AVHRR than for HIRS/MSU, tne latter 
being significantly better in turn than for SMMR. Cross correlations are gen- 
erally comparable for HIRS/MSU and SMMR (0.13-0.33), being much lower than for 
AVHRR (0.74-0.77). Scatter and cross correlation figures almost always improve 
substantially when the 3x3 smoother is used, but the relative rankings given 
above still hold. The 3x3 weighted smoother statistics for the AVHRR are 
quite respectable indeed, being 0.18-0.29C for scatter and 0.91-0.93 for cross 
correlation. 

Table 3-2 again haB adequate AVHRR bln matchup sampling in only some of 
the regions especially after the smoother is applied, and there are no FGCE 
buoys. Bias and cross correlation values are comparable for all the in situ 
data t>ets , but scatter is somewhat worse for the Transpac XBT than the Pazan 
ships, although even this reverses when the 3x3 smoother is used. As before, 
the 3x3 smoother greatly improves the AVHRR/ship comparisons of scatter and 
cross correlation (see remarks on Table 3-1 just above), but generally worsened 
the biases somewhat. 

3. March 1982 Data Period 

This is the last full month of operational MCSSTs prior to the erup- 
tion of El Chichon and about the beginning of the period when the 3.7 pm 
noise level began to climb significantly. The statistics in Table 3-1 indicate 
generally lower biases for the AVHRR than for the HIRS/MSU, SMMR, or VAS. The 
AVHRR scatter values are superior to those of the other sensors , and they are 
particuarly noteworthy when the 3x3 smootht r is used, viz., 0.21-0.29C (see 
Table 3-1). The AVHRR cross correlation values are comparable with those of 
VAS and substantially better than those for HIRS/MSU and SMHR, and as with 
previous periods, those correlations associated with the 3 x 3. smoothing are 
higher (reaching 0.58-0.84 for the AVHRR). 

Table 3-2, which has the same coverage restrictions as in the previous 
period, shows somewhat lower bias values for Transpac XBTs than for Pazan 
ships, and again has about comparable cross correlations and higher scatter 
values relative to the Transpac XBTs. 

4. July 1982 Data Period 

By this month the El Chichon volcanic aerosol cloud had girdled the 
Earth several times, but generally had remained just north of the equator, and 
this eevereiy reduced the numbers of daytime AVHRR retrievals in the region of 
the aerosol cloud. Nighttime MCSST observations in that same zone exhibit 
large positive biases (AVHRR lower than in situ data because of aerosol attenu- 
ation in the high stratosphere). Chief effects of the increased 3.7 pm noise 
level are reouced observational densities and somewhat higher scatter, the 
latter from occasional erroneous passing of the uniform low stratus test, both 
at night. Table 3-1 reflects these factors in the large positive AVHRR bias 
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values. Unfortunately, the 3x3 smoothing procedures results in a zero sample 
in the Mid-Pacific zone. AVHRR scatter is adversely affected also, making it 
higher than the VAS, comparable to the HIRS/MSU, but still rather better than 
the SMMR. AVHRR cross correlations are poorer than in the previous periods, 
but still comparable to or even better than for the other sensors. 

Table 3-3 shown AVHRR statistics are comparable relative to all the in 
situ data, and that all are generally worse than for the previous months 
studied. As with the other periods, scatter and cross correlation statistics 
are significantly improved by the 3x3 smoothing scheme, whereas the effects 
of the smoother on the biasing is mixed. 

5. Diurnal Variations 

Table 3-3 summarizes the AVHRR/ship comparisons for the three data 
periods when both day and night statistics are available. Excluding for the 
moment the month of July 1982, in which El Chichon had a strong influence that 
differed by day and night, there appears a systematic diurnal difference in the 
bias in the March period of 1982, larger positive (ship SST higher than MCSST) 
at night than in the daytime, that is not evident in December 1981. Standard 
deviations are somewhat lower, and cross correlations somewhat higher, at night 
for all data periods. Both these statistics generally improve substantially 
when the 3x3 weighted smoother is used. The rather drastic change from a 
moderate to a large negative bias in the daytime to an even larger positive 
bias at night during June 1982 reflects the impact of El Chichon. Monitoring 
of drifting buoy/MCSST matchups during the first half of 1982 also detected a 
nighttime positive bias of about 0.4C in the NOAA operational product. This 
bias was effectively removed after mid-September 1982 by re-derivation of the 
temperature-dependent bias correction (see discussion in Sec. F.5). 

6. Summary 

With the partial exception of the El Chichon impacted July 1982 period, 
the error characteristics of the operational AVHRR-based SSTs as exhibited by 
the statistical measures given in Tables 3-1 and 3-2 are predominantly as 
follows: biases are generally a few tenths of a degree and positive in sign 

(AVHRR lower than in situ temperatures) except for the earliest data period 
when the GOSSTCOMP instead of the MCSST method was in use. Scatter (i.e. 
standard deviation) magnitude tended to be in the range 0.5-0.6C, and when the 
latter are subjected to a 3 x 3 weighted smoothing scheme, the scatter further 
decreases to a remarkable 0.2-0. 4! Cross correlations, except in the tropics 
where the extremely small range in surface temperature generally mediates 
against high values for this type of statistic, generally fall in the range 
0.3-0. 7; use of the 3x3 smoothing results in a further increase to 0.5-0. 9. 

It helps to put the foregoing discussion of the AVHRR/MCSST matchup com- 
parisons with the Pazan ship and Transpac XBT data sets into better perspective 
if one looks at the statistics for the matchups between the Pazan screened ship 
observations and the Transpac XBT measurements (see Table 3-4). The biases 
associated with the Transpac XBTs tend to be one to three tenths of a degree 
smaller than those found with the AVHRR, but the standard deviations and cross 
correlations of the AVHRR are consistently and significantly better than those 
relative to the XBTs. As has been noted in the previous remarks, the 3x3 
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weighted smoother has variable results on che bias values for both data sets, 
generally worsening them somewhat, while substantially improving both scatter 
and cross correlation figures. The XBT standard deviations tend to have larger 
reductions than the AVHRR when the smoother is applied, thus rendering them 
quite comparable in magnitude. This was particularly the case in the July 1982 
data period when the scatter for the XBTs became smaller than for the AVHRR, 
presumably because of the adverse effects of El Chichon on the MCSSTs. 

The foregoing once again emphasizes the uncertainties in assessing the 
"true" accuracy of the satellite-derived SSTs when the correlative in situ 
measurements evidently contain significant but unknown errors of their own. It 
would be desirable to cross compare the Pazan ship, Transpac XBT, and FGGE buoy 
data secs in this connection, but the number of ship/buoy, and XBT/buoy bin 
matchups is N**l and N-14, respectively. The statistics for the latter are not 
particularly impressive: B - -0.49, S - 1.11, and C « -0.03, but confidence in 

them is low with such a small sample. 

Unlike the favorable effects on the other statistical measures for all 
data periods and regions, use of the 3x3 smoothing had mixed results on the 
bias values, generally tending to worsen them slightly. The notable exception 
is the July 1982 data period when the El Chichon eruption cloud belt adversely 
affected all the AVHRR statistics generally. Even in this case the large 
positive biases were decreased somewhat by the smoother. Use of only those 
two-degree bins having more than 20 ship observations appeared to have only a 
slight, and mixed, effect on all the statistical parameters. 

F. REGIONAL AVHRR ERROR CHARCTERISTICS FROM GLOBAL ANOMALY CHARTS 

When using the globai charts of SST anomalies or anomaly differences be- 
tween AVHRR and in situ data sources, it is necessary to be reminded that valid 
comparisons can be made only in areas where there is an adequate distribution 
of both types of observations. Whereas the density of AVHRR-based SSTs was 
generally good to excellent on a world-wide basis (the exception being the Nov 
1979 GOSSTCOMP distribution, which is poor south of 45S and in parts of the 
tropical belt, particularly the zone from 90E to 180E) , the ship SST distribu- 
tion for any of the data periods is of adequate density only in parts of the 
North Atlantic and North Pacific and a few narrow tracks elsewhere (Figure C-l 
is an example, see Appendix C). The Transpac XBT data coverage is even poorer 
than that of the ships, being concentrated along a few heavily travelled ship- 
ping lanes, and only a few of the two-degree bins have >4-6 observations for a 
given month (Figure C-l is an example). The data distribution for the P GGE 
buoys (Nov. 1979 only) is all in the Southern Hemispere, but it provides some 
spotty coverage in regions rarely visited by ships (see Figure C-l). 

1. November 1979 Data Period 

There is general agreement between AVHRR- and ship-derived anomaly 
patterns , but some differences in amplitude are evident (Figure C-2). General 
agreement with the anomalies derived from Transpac XBTs and from the FGGE buoys 
(Figure C-2) is apparent also, though again not on a bin-by-bin basis. Some of 
the >3.5C positive anomalies seem suspect, particularly those on the edge of 
the AVHRR data void in the Southern Hemisphere. Parts of the South Atlantic 
area also appears too warm in the AVHRR. 


A*** 
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better than it was with the AVHRR/GOSSTCOMP in November 1979, particularly in j 

the central and western North Pacific, central and eastern South Pacific, and : 

the high latitudes of the Southern Hemisphere generally. There is still a i' 

paucity of observations in the deep tropics of the western Pacific (Indonesia/ 

Micronesia areas). As in the November 197J period, if one looks at the larger, ! 

more coherent anomaly features, there is good general agreement between the | 

satellite-derived pattern and the ship-derived pattern (Figure C-3), although 
the overall amplitude of the MCSST anomaly field tends to be greater. The 
strength of the large positive MCSST anomaly northeast of New Guinea, as well 
as a weaker one northeast of Madagascar, and a strong ship-based anomaly south- 
east of the tip of Africa, all appear suspect. 

3. March 1982 Data Period 

Although MCSST data densities during this period are poorer in the 
northern and eastern North Pacific than in December 1981 , they are higher in [ 

the Indian, South Atlantic, and Southern Oceana. As previously, the Indonesia/ 
Micronesia area is relatively poorly observe*. Once again, the larger scale 
patterns generally match, but the amplitudes of the MCSST maxima are generally 
somewhat greater than the ship ones (Figure C-4). The positive MCSST anomaly 
in the Indonesia/Micronesia area and westward appears to have no counterpart 
in the ship-based pattern. 

4. July 1982 Data Period j 

< 

i 

The drastic impact of the El Chichon aerosol cloud and the increased 
3.7 pm noise on MCSST observational densities is seen in Figure C-5. Excel- 
lent coverage remains, however, in the central North Atlantic and in the 
Southern Hemisphere down to at least 45S. If one disregards the "El Chichon 
negative anomaly" stretching around the globe between roughly 5N-30N, there is 
yet again reasonably good conformance betveen the MCSST-based and ship-based 
(Figure C-5) anomaly fields, although the amplitude of the MCSST negative 
anomaly in the mid North Pacific is somewhat too large. 

5. Diurnal Effects 

Global AVHRR day-minus-night charts (Figure C-9) were produced for 
the Dec. 1981, March 1982, and July 1982 data periods. Although some syste- j 

matic difference patterns are evident, their interpretation or explanation is ; 

hampered because no separate day and night anomaly charts were produced. ' 

Unfortunately, little recourse can be made to the separate day and night I 

tabulations in the statistical tables (e.g., see Table 3-4), as they are 
broken out only by rather large regions. 

A pronounced positive day-minus-night difference of several degrees ex- 
tends around the Earth between about 25-40S in the December 1981 data period. j 

This positive difference appears in a far weaker, more irregular or interrupted 
form and generally at somewha; lower latitudes of the Southern Hemisphere in 
the March 1982 period. It is absent there altogether in July 1982, but appears 
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In similar latitudes of the Northern Hemisphere that month. This seasonal be- 
havior, probably associated with northward shift of the moBt intense solar 
heating from the Southern Hemisphere at the time of the summer solstice there 
to the Northern Hemisphere at the time of its summer solstice, is consistent 
with a widespread "diurnal thermocline" developing in the weak wind regimes of 
the central subtropical anticyclone belts. 

The other major day-to-night difference that appears during all three per- 
iods is a positive one extending westward from the Indonesia area into the cen- 
tral Indian Ocean near the equator. In December 1981, this is clearly related 
to too large an atmosphere correction from a quadratic term in the operational 
daytime split-window equation (see discussion in Sec. G.2). In March and July 
of 1982, however, a linear form of the split-window equation was in daytime 
use, so an alternate explanation must be sought for what appears to be a per- 
sistent overcorrection for very high moisture. When the coefficients of the 
operational daytime and nighttime equations were adjusted slightly in October 
1982 on the basis of a larger and more representative sample of buoy/MCSST 
matchups (see McClain in Appendix C, pp. C-l to C-16, JPL, 1983), one of the 
results was to diminish this tendency toward a negative daytime and positive 
nighttime bias in those regions of the tropics that are extremely moist, and 
thus to lessen diurnal differences of the type noted here. 

6. Summary 

The larger-scale MCSST anomaly patterns are In fair to good agree- 
ment with the corresponding ship-derived fields almost everywhere and during 
every data period where there was adequate common coverage. Looking only at 
the small-scale, bin-to-bin values, it is difficult to see the relatively 
large degree of pattern similarity that really exists overall. In isolated 
areas there are suspect anomalies during the several data periods, and the 
amplitude of the AVHRR anomaly field seems somewhat larger overall than the 
ship field, but generally there is surprisingly good correspondence in view of 
the probable errors inherent in both fields. These errors combine (adding or 
subtracting) in a variable and unknown way in their difference field. The 
diurnal variations apparent in the global day-minus-night charts are sometimes 
difficult to account for in the absence of separate day and night anomaly 
fields, but the prevalence of positive day/night differences can often be 
attributed to the "diurnal thermocline" effect in areas where lack of wind has 
Inhibited mixing of the uppermost layer. 


G. POSSIBLE CAUSES FOR ERRORS IN THE AVHRR-BASED ANOMALY FIELDS 


► * 


The general subject of sources of error in SSTs obtained from the AVHRR 
using the GOSSTCOMP or MCSST methods is addressed in Section C, but a few 
additional remarks directed specifically to possible errors noted in the above 
discussion of anomaly fields are in order. 

1. November 1979 

The tendency for a positive anomaly rimming the edge of the data at 
A5-50S (Figure C-2) is probably something peculiar to the GOSSTCOMP method; as 
it does not appear in any of the three MCSST data periods. On the other hand, 
the rather scattered two-degree bins with ship-based anomalies do tend to 
support the presence of warmer than normal water in that part of the southern 
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Indian Ocean (Figure C-2). The relatively cooler water extending from near the 
southwest tip of Africa northwestward into the central South Atlantic is sup- 
ported in the AVHRR/GOSSTCOMP anomaly field (Figure C-2) only by a tongue of 
relatively less warm water. This could be a "diurnal thermocline" effect, as 
this area of the subtropical Atlantic was characterised by anticyclonic con- 
ditions with light winds at the surface and lew amounts of cloudiness 5 in 
November 1979. There is no obvious explanation for the lack in the AVHRR of 
any real indication of the strong negative anomaly evident in the ship data 
some distance to the east of the southern part of South America. 

2. December 1981 

The large positive anomaly northeast of New Guinea, and to a lesser 
extent the weaker ones just northeast of Madagascar and along the northeast 
coast of Australia (Figure C-3) do not appear to be supported by the ship-based 
anomaly field (Figure C-3), although ship data were very scarce in the first 
area mentioned. The Australian case would seem to be a manifestation of the 
"diurnal thermocline" in the IR-based SSTs. This area was under weak mean 
monthly pressure gradients near sea level (light winds), and there was a mean 
monthly albedo of <20%, indicating very little cloudiness. The New Guinea 
and Madagascar cases may also have an element of this effect, but the situation 
there was aggravated by the use of a quadratic term in the daytime split-window 
MCSST equation (see Sec. F.5). This term was found to produce erroneously high 
MCSSTs, but only in areas where atmospheric water vapor was exceedingly large 
(i.e. precipitble water >5 cm), a characteristic of the region extending 
westward from Micronesia/Indianesia into the central Indian Ocean. Further 
confirmation of a moisture maximum there is afforded by a SMMR-derived precip- 
itable water chart for December 1981 provided by NASA/GSFC. The positive anom- 
aly southeast of South Africa seems rather too large in magnitude in the ship- 
derived field, but there is a large amplitude and equally dubious negative 
anomaly in the MCSST field just to the south of it in a large region devoid of 
ship data. There is no obvious explanation of either of these. 

3. March 1982 

Aside from the previously noted tendency for greater amplitude in 
the MCSST-derived anomalies than in the ship-based ones (see particularly the 
negative anomaly features in the North Pacific and North Atlantic (Figure C-4) , 
the only unsupported feature is the positive anomaly in the far western Pacific 
and Indian Oceans centered roughly on the equator. This corresponds 
cllmatologically to the moistest portion of the tropics, and this is 
corroborated by the SMMR-based precipitable water charts supplied by 
NASA/GSFC. A similar, but rather stronger, positive anomaly in December 1981 
was largely accounted for by use of a quadratic daytime MCSST equation, 
discontinued in February 1982. The fnc t that an apparently erroneous positive 
anomaly persists in March and July 1982 indicates that the limited set of buoy 
matchups used in deriving the original temperature dependent bias corrections 
to the MCSST simulation equations evidently still did not adequately account 


5 This information obtained or inferred irorn charts of monthly mean pressure 
at sea level from NCAR, monthly mean SMMR wind speeds (NASA/GSFC), and of 
monthly mean albedo and outgoing long wave radiation from NOAA/NESDIS. 
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for Che atoo sphere/ocean conditions In this region. This positive anomaly 
vanished when re-derlvcd bias corrections, based on a ouch larger and more 
representative buoy matchup data base, were incorporated in the operational 
equations in October 1982. 

4. July 1982 

The large belt of negative anomaly stretching around the Earth from 
roughly 5-30N is, as discussed previously, almost entirely a consequence of the 
volcanic aerosol from eruptions of El Chichon in Mexico in early April 1982. 

The cold aerosol cloud almost eliminated daytime MCSSTs In this belt and 
severely attenuated the emitted radiation from the Earth's surface. The dif- 
ference between monthly mean SSTs derived solely from satellite MCSSTs and 
solely from ship and other in situ measurements has been used by Strong et al. 
(1983) to track the month-to-month coverage of the volcano cloud during 1982. 
His charts indicate a southward transport of the aerosol near central West 
Africa, and a northward transport In the west central Horth Pacific, which is 
consistent with distortion in the negative anomaly belt In those regions (Fig- 
ure C-5). The ship-based anomaly field (Figure C-5) also shows negative 
anomaly features in those two areas, but the greater extent and amplitude in 
the AVHRR chart is attributed to the added influence of tne volcano cloud. The 
source of the large negative MCSST anomaly along a large part of the extreme 
southern edge of the chart, especially southwest of Australia, la suspected to 
be deficiencies in the climatology there (the ship-derived anomaly chart has a 
data void all through that region). The N0AA/NESD1S anomaly charts, which are 
based on the Robinson/Bnur climatology, show only two small (in area) negative 
anomalies (maximum of -1.50 anywhere In this zonal belt, one southwest of 
Australia at about 50-53S and another southeast of New Zealand near 55-58S. 

The general tendency, in all data periods, for the amplitude of the AVHRR- 
based anomalies to be somewhat greater than the ship-derived anomalies, is pos- 
sioiy a consequence of to factors; (1) the sparsei and more irregular sampling 
of ship data going into the bin averages; and (2) the skin temperature. 

H. POSITIVE FINDINGS 

Despite the difficulty of comparing satellite sensors with differing spa- 
tial resolution and geographical covcrge, and matching each of these in turn 
with common in situ data sets that are themselves highly irregular in density 
of coverage in many regions and are of variable quality (furthermore, they are 
"spot" measurements at depths of one to several meters, whereas the AVHRR 
senses a kind of "skin" temperature over an area about ten kilometers on a 
side), there is surprisingly good correspondence in the sign and location of 
the major anomaly features; the correspondence in the amplitude of the anomaly 
maxima is only fair. 

As expressed by the bins, scatter, and cross-correlation statistics for 
the globe and for the North Pacific, Mid-.'*icif ic , South Pacific, and North At- 
lantic regions, the operational NOAA/NESDIS AVHRR product, the MCSST, generally 
compares better with the Pazan screened ship data set than do any of the other 
satellite sensors. This is especially true in the case of the standard devi- 
ation (scatter) and the cross correlation, with the AVHRR having been lower 
scatter and higher correlation with respect to the Pazan ships than do the 
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Transpac XBTs. This, and the fact that the AVHRR similarly has better statis- 
tics relative to the Pazan screened ships than it has relative to cither the 
Transpac XbTs or the FCGE buoys, suggests that the latter two data sets should 
have had an equally comprehensive screening to delete incorrect observations. 

When a modest amount of additional spatial smoothing (using the 3x3 
center-weighted smoother) is done, the AVHRR/MCSST figures for scatter and 
cross correlation become impressive even on an absolute basis, attaining 
values of 0.2-0.4C for the former statistic and 0.7-0. 9 for the latter. This 
high a correlation is noteworthy when one is reminded that it refers to . 
anomaly values, which have a much smaller range than the temperature values 
from which they were derived. 

1. RECOMMENDATIONS FOR IMPROVED COMPARISONS OF SENSOR AND IN SITU DATA 

In view of the substantial improvement in statistics that resulted fron 
application of the 3x3 center-weighted smoother, with its implications about 
the noise level remaining in the bin averages of ship and some of the sensor 
>iata, three recommendations are made: (1) Find a way to include 3x3 arrays 

having up to a total of say four side and/or corner bins missing from the 3x3 
array; this would cut down the erroneous reduction in sample size that results 
from the present way of doing the 3x3 smoothing (i.c., no bin can be miss- 
ing). (2) Try using 4 x 4 deg. lat./long. bins instead of 2 x 2 in the basic 
monthly mean anomaly chart. (3) If the effect on the statistics is salutory 
after doing either of the above, produce new global anomaly charts, either 
color-coded or contoured. 

If global day-minus-night charts arc to be properly interpreted, then it 
is necessary to produce separate daytime and nighttime global anomaly charts. 
Otherwise, one cannot tell if a feature in the difference field comes about 
because the daytime values are high, or the nighttime values are low, or both, 
or vice versa. 





SECTION IV 

RETRIEVAL OF SEA-SURFACE TEMPERATURES FROM HIRS2/MSU 

J. Sueskind and D. Reuter 
Goddard Laboratory for Atmospheric Sciences 

The methods used at the Goddard Laboratory for Atmospheres (GLA) to 
retrieve sea surface temperatures from HIRS2/MSU data for the four monthB of 
the NASA sea surface temperature intercomparison workshop are described. 

Results are shown comparing anomaly fields produced using data from ships, 
AVHRR, HIRS2/MSU, SMMR and VAS for the last three of these months. Fields 
from AVHRR and HIRS2 show the highest accuracy compared to ship fields. 

Errors in the HIRS2 fields appear more random while AVHRR data shows large 
area, spatially coherent errors. The random errors in the HIRS2 fields can be 
further reduced by performing the retrievals at a higher spatial resolution. 

A. OVERVIEW OF THF GLA RETRIEVAL SCHEME 

HIRS2 and MSU are the 20 channel infrared and 4 channel microwave passive 
sounders on the operational, low earth orbiting satellites. They monitor emis- 
sions, arising primarily from the earth's surface and the atmosphere up to the 
mid stratosphere. These, together with the SSU, a three channel pressure modu- 
lated infrared radiometer which monitors emission from the mid-upper strato- 
sphere, comprise the TOVS (TIROS Operational Vertical Sounder) system. 

The TOVS data are analyzed operationally by NOAA NESDIS to produce verti- 
cal temperature-humidity profiles using a method based primarily on statistical 
regression relationships between observed radiances and atmospheric parameters 
(Smith 1980). The approach used at GLA is fundamentally different from the 
current operational approach. Rather than rely on empirical relationships be- 
tween observations and meteorological conditions, we attempt to find surface 
and atmospheric conditions which, when substituted in the radiative transfer 
equations describing the dependence of the observations on the meteorological 
conditions, match the observations to a specified amount. A physically based 
retrieval scheme has a number of advantages over a statistically based scheme. 
The single most important advantage is the ability to correct for the effects 
of auxiliary factors such as surface temperature, surface emissivity, surface 
elevation, reflected solar radiation, satellite zenith angle, and most signifi- 
cant of all, clouds on the observations. All of these parameters are either 
solved for, or directly accounted for, together with the atmospheric tempera- 
ture profile, in an iterative scheme. As a result of this, the data analyzed 
produce not only global fields of atmospheric temperature profiles, which are 
necessary for the initialization of atmospheric models for numerical weather 
prediction, but also other auxiliary fields. These include the following 
monthly mean fields: sea/land surface temperature and their day-night dif- 

ference which, over land, is related to soil moisture; fractional cloud cover, 
cloud-top temperature, and cloud-top pressure and their day-night differences; 
and ice and snow cover which is derived from combined use of the surface emis- 
sivity at 50.3 GHz and the ground temperature. Another important advantage of 
the physical retrieval is the ability to identify those areas where a satis- 
factory solution to the radiative transfer equations cannot be found, in which 
case the retrieved parameters are flagged as questionable. 
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The decalls of an earlier version of Che GLA recrleval algorithm, as well 
as results of analysis of data from January 1979, are given in Susskind et al. 
(1984). The GLA processing syscem consists of the following main steps: 1) 

the forward calculation of radiances as a function of atmospheric and surface 
conditions; 2) accounting for cloud and surface radiative effects on the obser- 
vations; 3) determining atmospheric temperature profiles by the inverse solu- 
tion of the radiative transfer equation; and 4) determination of auxiliary 
meteorological parameters. Steps 1-3 are performed iteratively until conver- 
gence is reached. The effects of clouds on the infrared observations muBt be 
accounted for before either ground or atmospheric temperatures can be deter- 
mined. The cloud-filtering algorithm utilizes MSU channel 2, which has a 
transmittance of about 0.1 at the surface and is sensitive to the surface 
emissivity. Therefore, after the calculation of radiances, microwave surface 
emi8sivity is calculated next in the iterative scheme, using MSU channel 1, 
followed by the cloud correction, the retrieval of ground temperature, and 
finally, the update of atmospheric temperature profile. If sufficient agree- 
ment between observed and calculated radiances is found, the procedure is 
terminated and step 4 is performed. Otherwise, the iterative procedure is 
continued with recalculation of radiances, surface emissivity, etc. 

Table 4-1 shows the channels, centers, and peaks of the weighting func- 
tions dr/dinP, and radiance contribution function BdT/dinP, and other relevant 
information, for the channels on MSU and HIRS2. The current analysis docs not 
employ the SSU observations. Those channels utilized by GLA in analysis of 
the data and their primary function are indicated. 

The months of data for the sea-surface temperature workshop were November 
1979, December 1981, and March and July of 1982. HIRS2/MSU data for the first 
month were taken from TIROS-N. The remaining data were from NOAA-7 , on which 
the channels of HIRS are slightly but significantly different from those on 
TIROS-N. Consequently, one change had to be made to the processing system 
described in Susskind ct al. (1984) for the NOAA-7 data. Other changes and im- 
provements in the processing were also made during the course of the workshop. 

The data from November 1979, the first month of the workshop, would have 
been analyzed exactly as described in Susskind et al. (1984), but essentially 
one half of the month of data were missing including a big gap from November 10 
to November 17 and a number of smaller gaps. Therefore, in order to have more 
data points and cut down the random noise component of the monthly mean fields, 
retrievals were run in both the warmest (least cloudy) and second warmest 
(second least cloudy) 125 x 125 km quadrants of the 250 x 250 km grid (see Fig- 
ure 4-1), rather than in only the warmest quadrant as had been done earlier. 

Significant changes were made to the processing system used to analyze the 
data for December 1981 and March 1982. One change, affecting the clear column 
radiance algorithm, was made primarily because of the change in characteristics 
of channel 13 and channel 14 of KIRS2 on NOAA-7 from those of the same channels 
of HIRS2 on TIROS-N. The improved algorithm, shown in the next section, has 
been found to be superior even with TIROS-N data, and has now been incorpo- 
rated, together with other changes, for use in the re-analysis of TIROS-N data. 
The Becond change involved use. of the 1 1 pm window data, in addition to the 
3.7 and 4.0 um window data, in the retrieval of ground temperatures. In 
addition, individual soundings were assigned weights to be used in the genera- 
tion of monthly mean fields. After studying the December 1981 and March 1982 
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Table 4-1. HIRS2 and MSU Channels 


Channel 

A (pm) 

v(cm -1 ) 

Peak of 
dx/dlnp 
(mb) 

Peak of 
BdT/dlnp 
(mb) 

HI 

14.96 

668.40 

30 

20 

H2 a 

14.72 

679.20 

60 

50 

h3 

14.47 

691.10 

100 

100 

H4 a * c 

14.21 

703.60 

280 

360 

H5 C 

13.95 

716.10 

475 

575 

H 6 C 

13.65 

732.40 

725 

875 

H7 C 

13.36 

748.30 

Surface 

Surface 

H 8 b » c 

11.14 

897.70 

Window, sensitive to water vapor 


H9 

9,73 

1027.90 

Window, sensitive to O 3 


H10 

8.22 

1217.10 

Lower tropospheric water vapor 


HI 1 

7.33 

1363.70 

Middle tropospheric water vapor 


H12__ i 

6.74 

1484.40 

Upper tropospheric water vapor 


H13 a, d 

4.57 

2190.40 

Surface 

Surface 

Hl4 a * d 

4.52 

2212.60 

650 

Surface 

H15 a 

4.46 

2240.10 

340 

675 

H16 

4.39 

2276.30 

170 

425 

H17 

4.33 

2310.70 

15 

2 

H18 b 

3.98 

2512.00 

Window, sensitive to solar radiation 


H19 b 

3.74 

2671.80 

Window, sensitive to solar radiation 


Ml e 

0.596* 

50.30** 

Window, sensitive to surface emlssivity 

M2 d 

0.558* 

53.74** 

500 


M3 a 

0.546* 

54.96** 

300 


M4 a 

0.518* 

57.95 

70 



a used by GLA to compute temperature profiles 
bused by GLA to compute surface temperature 
c used by GLA to compute cloud fields 
d used by GLA in cloud correction 
e used by GLA to compute surface emlssivity 
*\ in cm 
**v in GHz 


sea surface temperature fields in Workshop III, a number of other changes were 
made to further reduce noise in the field. The newest system, called HIRS 
version 2 in the workshop, was used to reprocess March 1982 data and to 
process July 1982 data. For this reason, results for July 1982 and the 
reprocessed data for March 1982 are better indicators of the capabilities of 
HIRS2/MSU for retrieval of surface temperature than those of the earlier 
months. Modifications to Susskind et al. (1984) used in the analysis of 
workshop data for December 1981 and March 1982, and in the analysis of July 
1982 and re-analysis of March 1982, are described in the next two sections. 
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IMPROVEMENTS TO THE SYSTEM USED IN ANALYSIS OF DECEMBER 1981 AND MARCH 
1982 DATA 


HIRS2/MSU data for December 1981 and March 1982 were analyzed exactly as 
in Susskind et al. (1984), except for modifications which were made to improve 
the clear column radiance algorithm and the sea/land surface temperature 
algorithm. These are described below. 

1. Improved Estimation of Clear Column Infrared Radiances 

Infrared observations are highly sensitive to the presence of clouds in 
the field of view and, hence, accounting for their effects on' the observations 
becomes perhaps the most important step in the retrieval process. If one looks 
at an otherwise homogeneous but partially cloudy scene with cloud fraction o, 
then to a reasonable approximation, one can write the radiance observed in 
channel i to be 


R i - (1-a) Ri, C LR ♦ ® -i.CLD 


(4-1) 


where R^ is the radiance one would observe in a cloud-free area and R^ 
is the radiance one would observe in a completely cloud-covered area. Qyj 
depends not only on the atmospheric variables but also on the detailed proper- 
ties of the clouds. Rather than assume or attempt to determine the cloud pro- 
perties simultaneously with the determination of atmospheric and surface pro- 
pertie? , the method attempts to estimate, or "reconstruct" from the observed 
radiances, the clear column radiances which would have been observed if no 
clouds were present. These reconstructed clear column radiances, Rj[, are used 
in determiation of the atmospheric temperature profile as well as in the deter- 
mination of sea surface or ground temperature. The cloud field parameters are 
determined only after a complete atmospheric and surface solution is obtained. 
This method of treating cloud3 is fundamentally different from the approach 
used in analysis of AVHRR data, in which high spatial resolution is used to 
attempt to identify clear spots. Sea-surface temperatures are not determined 
in AVHRR spots thought to have any cloud contamination. 

A two f ield-of-view approach, similar to one originally introduced by 
Smith (1968), is used to^extrapolate observed radiances to obtain reconstructed 
clear column radiances, R^. We express Rj as 


Ri - Rij + n[R ltl - Rj. ,2 1 (4-2) 

where R^ ^ is the observed radiance for channel i in the field of view having 
the larger 11 pm radiance, and R^ 2 I s the observation of channel 1 in the 
second field of view. 


If R^ CLR is known for a given 
according to 


channel , 


then n can be solved for 
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R t,CLR ~ R l,l 
R l,l “ R i ,2 


Since n “ dj/(aj - 02 ), it should not be dependent on the channel used. In 
Susskind et al. (1984^, n was determined by combined use of HIRS2 channel 13 
and MSU channel 2. Ri3,CLR is computed in each iteration (N) using the Nth 
iterative temperature profile and ground temperature. If the guess is 
too warm (or cold), would be too large (or small). The effect on the 

computed brightness temperature of channel 13 of a bias in the temperature 
profile in the troposphere is accounted for, to first order, by modifying the 
brightness temperature according to 


0-0 + 
13 13.CLR 


9 M2 " 0 M2 


(4-4) 


where 0^3, CLR ls the equivalent brightness temperature to R^3 t CLR» ®M2 is 
the observed brightness temperature in MSU channel 2, which is sensitive to the 

average tropospheric temperature, and 6^2 the computed brightness 
temperature in MSU channel 2. This procedure works reasonably well even though 
channel id is sensitive primarily to radiation much closer to the surface than 
that of MSU channel 2. Then, nj 3 is computed according to 


n 13" 


B 1 3 1 °1 3 1 ~ R 13 , 1 
R 13,l “ R 13,2 


(4-5) 


where B(0'), the black body function of the corrected equivalent brightness 
temperature, is the corrected estimate of the clear column radiance. While 
results using equations (4-2) through (4-5) are quite good, it has been found 
that in cases where the initial guess has a lapse rate error, improved results 
are obtained by defining in an analogous manner to O 13 in equations 
(4-4) and (4-5) and setting n equal to the average of ^3 and weighted 
by the square of the difference in radiances for each channel in each field of 
view: 


n i3 (R 13,r R 13,2 ) + n i4 (R 14,l~ R 14,2 )2 

(R 13,1~ R 13,2 ) + (R 14,f R 14,2 )2 


(4-6) 


including channel 14 radiances in the determination of n has the effect of 
utilising a single infrared channel with a broader weighting function, more in 
line with that of the microwave channel, to correct for cloud effects. Com- 
bined use of channels 13 and 14 was especially necessary for N0AA-7 data be- 
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cause the weighting functions for both channels 13 and 14 were shifted lower in 
the atmosphere relative to TIROS-N and channel 13 provided even a poorer natch 
with HSU channel 2. The procedure gives improved results with TIROS-N data as 
well* 


2. Improved Determination of Sea-Surface Temperatures and Ground 
Temperatures 


Once n is obtained in any iteration, reconstructed clear column rad- 
iances, r£, are obtained for each channel using equation (4-2). These radi- 
ances are used to get the updated estimates of ground temperature and atmo- 
spheric temperature profile. One can solve for ground temperature at night 
from any window channel observation (8, 18, or 19) in a straightforward manner 
given the clear column radiance, the estimated temperature-humidity profile, 
and an estimate of the surface emissivity. During the day, solar radiation 
reflected off the ground contributes significantly to the observed radiances in 
the short-wave window channels 18 and 19. A sea surface or ground temperature 
can be obtained during the day by simultaneous use of the two shortwave channel 
observations, with the additional assumption that the surface bidirectional 
reflectance of radiation from the sun, in the direction of the satellite, 1 b 
the same at both 4.0 and 3.7 pm (Susskind et al. 1984). 

In Susskind et al. (1984) only the two short wave channels were used to 
obtain both day and night global surface temperatures. The 11 pm channel was 
not used because attenuation due to water vapor becomes very significant in 
humid atmospheres. The differences between the retrieved January 1979 monthly 
mean day and night sea surface temperature were almost all less than 1°C as 
expected. This gives evidence that the procedures used to correct the short 
wave observations for solar contamination are valid. 

While the results using this method were quite good, several new improve- 
ments have been made to the system. The first change involves including radi- 
ances in the 11 pm window channel 8 in the estimate of sea/land cjrface 
temperature. As alluded to earlier, under very humid conditions, significant 
attenuation of radiation leaving the surface by water vapor absorption at 
channel 8 frequencies produces a large potential source of error in the deter- 
mination of ground temperatures both because of uncertainties in water vapor 
distribution, and even more significantly, low response of the radiances to 
changes in surface temperatures. Nevertheless, inclusion of 11 pm radiances 
introduces additional surface temperature information, which becomes quite 
accurate at low and moderate water vapor conditions. In addition, during the 
day, under some conditions, the long wave observations are superior to short 
wave observations which have a potential error source due to the reflected 
solar radiation. In order to determine the proper mix of channels to be used 
to get ground temperatures in a given situation, ground temperatures were re- 
trieved for all channels, that is, three independent estimates, T s g, T s ig, 
and T 8 jg, were made at night, and two independent estimates, T g g and T g .g j,, 
were made during the day. Each estimate was given a weight, Q^, which de- 
creased according to the magnitude of the cloud correction and atmospheric and 
solar radiation corrections that had to be made to obtain the ground tempera- 
ture T s i from the observed brightness temperatures. The weight was defined 
as 
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where j i8 the difference between the reconstructed brightness tem- 
perature and that observed in field of view 1 for channel i, and 0^ - T s ^ 

is the difference between the reconstructed clear column brightness temperature 

used to obtain a surface temperature and the surface temperature obtained. The 
first term represents the cloud correction and the second factor represents the 
atmospheric and solar radiation correction. 

In the case of daytime retrievals, the brightness temperature of channel 
18, which is less affected by solar radiation than that of channel 19, is used 

in equation (4-7) together with Tjg jg, to get Qjg jg. For a given sound- 

ing, an estimate of surface temperature is given a’zero weight if its weight as 
obtained from equation (4-7) is not at least half the average of the rest of 
the weights. In addition, any weight less than 1/30 was set equal to zero. No 
surface temperature is retrieved if all weights are zero. Otherwise, the sur- 
face temperature T„ is taken as the weighted sum of T s ^ 

\ ■ <| «t \,J'\ «i <4 - 8) 


The entire sounding is also given a weight, W, which reflects the number of 
independent estimates of sea surface temperature as well as the weight of each 
estimate 


W - II Q. 2 ] 
i 1 


1/2 


(4-9) 


Tgir. weighted by W in generating monthly mean fields. Given a set of sea- 
surface temperatures Tj and weights W, for a grid box k, the (weighted) 

monthly mean sea surface temperature anomaly field should be produced 
according to 


-U 

\ 


1 “j !i j - “j 


(4-10) 


where T CL TM,j is the Reynolds (1982) sea surface climatology interpolated in 
time and space to the location of the sounding J. As a consequence of equa- 
tions (4-7) through (4-10), more optimal use is made of the three window chan- 
nels for a given sounding, and relatively clear and/or dry soundings are 
weighted more than cloudy or humid ones in generating monthly mean fields. In 
addition, night-time soundings are given more weight than day-time soundings 
(not necessarily a good result) because night soundings have 3 independent 
estimates instead of 2, which tends to increase W, and also because the night- 
time values of Q tend to be higher in shortwave channels than the day values. 
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W, the weight for each sounding, equation (4-9), was sent to the workshop 3 

along with the corresponding T s , equation (4-8), so that weighted anomaly i 

fields, equation (4-10), could be produced. Unfortunately, the workshop did 
not initially utilize the weights provided to them in generating the monthly 
mean anomaly fields and treated all soundings as having equal weight. 

Figure C-3 (see Appendix C) shows the anomaly fields derived by the work- 
shop from data from ships, AVHRR, HIRS and SMMR for December 1981. The small 
scale noisy nature of the HIRS field is evident. In addition, large spurious 
warm anomalies are observed in the HIRS field in the vicinity of land and 
especially in the Gulf of Mexico, the Sea of Japan, and the Yellow Sea. 

Similar characteristics are found in the March 1982 anomaly field generated 
using this method. The errors close to land were found to be the result of 
the inadvertent use of a 4° x 5° topography which was used to define land 
and water in the retrieval program. This land-water flag is used to specify 
the surface emissivity, which is taken as higher over ocean than over land. 

Using an emisrivity with too low a value in analysis of data off the coasts of < 

continents resulted in spuriously warm retrieved surface temperatures. It was 
also observed i-. Workshop III that the small scale noise in the HIRS fields is j 

greatly reduced by applying a 9 point smoothing to the anomaly fields I 


K ■ i \i s k i < 4 -“> 


where grid points 1 are adjacent to k and is a smoothing matrix. As a 
result of this, the workshop generated some statistics for smoothed anomaly 
fields derived from all sensors and ships. This produced some interesting 
findings which will be shown later. 

C. MODIFICATIONS MADE TO REPROCESS MARCH 1982 AND PROCESS JULY 1982 - HIRS 

VERSION 2 

As a result of the findings on the HIRS fields originally produced for 
December 1981 and March 1982, a number of further modifications were made to 
the program. The first change involved simply replacing the 4° x 5° top- 
ography by a 1° x 1° topography. This eliminated the large errors near 
the coasts. The second change involved the processing of more data to reduce 
the effects of random noise. It was found that retrievals performed in each 
of the three quadrants having the largest brightness temperatures (see Figure 
4-1) produced geophysical parameters of comparable accuracy to those obtained 
only from analysis of data from the warmest quadrant as seen by the 11 pm 
window. This tripling of the data density did little to change the sea-surface 
temperature anomaly patterns, but greatly reduced the random noise component. 

In addition, improved fields of ground temperature and ice and snow cover 
resulted from the increased data density. 

A consistency check was also added to the sea surface temperatures T s> i 
obtained for a given sounding with channel i. Ir the case of three estimates 
of sea surface temperature for a given sounding, Qj was set equal to zero if 
T 8 £ was different from T g , computed in equation (4-8) uning , by more 
than 1.5°C. In this case, T s was recomputed using equation (4-8) with 
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Qj_ * 0. In the case of two estimates of T a ^ for a given sounding, the 
less reliable sounding was eliminated if the difference between the two sound- 
ings was greater than 3°C. The less reliable sounding was defined as the 
one differing the most from climatology. This does not eliminate large dif- 
ferences from climatology as long as the result is supported by more than one 
channel for a given sounding. This system, called HIRS version 2, was used to 
reprocess March 1982 and process July 1982 data. As in the system described 
in Section A, equations (4-1) through (4-9) were used and temperatures and 
weights were sent to the workshop. Comparison of results for December 1981, 
March 1982, and July 1982 with those produced from ship data are shown In the 
next section for H1FJ5 fields and other fields in the workshop. 

D. COMPARISON OF RETRIEVED ANOMALY FIELDS WITH THOSE OBTAINED FROM SHIP DATA 

One way of judging the accuracy of sea surface temperatures retrieved by 
various sensors is to compare characteristics of the monthly mean anomaly 
fields derived from the sensors and from ship measurement. Climatology is 
also treated as a sensor in this comparison. Tables 4-2 to 4-4, taken from 
values generated by the JPL workshop, show: C, the correlation coefficient; 

B, the bias (sensor-ship); S, standard deviation; and N, the number of grid 
point samples (2° x 2° grid) comparing anomaly fields from a number of 
sensors with ship fields. Statistics are given for global colocations and 
also for those in the North Pacific and North Atlantic Oceans. Results are 
given for the unsmoothed field, and also for 9-point smoothed fields, in which 
case, both the ship field and retrieved field were smoothed and then compared 
to each other. All 8 points surrounding a grid point were needed to perform 
the smoothing. Therefore, only points for which unsmoothed anomalies were 
available for all surrounding points were included in the smoothed statistics. 
As a result of this, many ship and XBT observations were dropped from the 
smoothed fields because these fields have more data gaps, especially in the 
tropics and southern oceans. 

To interpret the significance of standard deviation from ships, an anomaly 
field produced by a sensor can be considered skillful if its standard deviation 
irom ships is at least as low as that of climatology, which is representative 
of the ocean signal. Statistical results for December 1981 are shown in Table 
4-2. It is interesting to see that climatology differs from ships with a 
standard deviation of about 0.6°C when no smoothing is applied and about 
0.4°C when smoothing is applied. This drop in "signal" is most likely due 
to noise in the unsmoothed ship field. It is also possible that further 
deviations from climatology exist in those grid points which are excluded from 
the statistics because ships did not report in all the surrounding grid points, 
but it is unlikely that this would explain all the reduction in standard 
deviation. 
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Table 4-2. Comparison of M^ithly Mean Anomaly Fields with Ships > 5/Cell, 
December 1S81 





Unsmoothed 



Smoothed 




Global 

N. Pac. 

N. Atl. 

Global 

N. Pac. 

N. Atl. 


C 

0.76 

0.77 

0.74 

0.91 

0.93 

0.93 

AVHRR 

8 

- 0.30 

- 0.44 

- 0.15 

- 0.33 

- 0.43 

- 0.19 


S 

0.50 

0.50 

0.41 

0.28 

0.29 

0.18 


N 

729 

376 

255 

235 

127 

102 


C 

0.21 

0.29 

0.13 

0.45 

0.56 

0.04 

H1RS 

B 

0.13 

0.31 

0.10 

0.21 

0.21 

0.24 


S 

0.88 

0.89 

0.77 

0.42 

0.45 

0.39 


N 

729 

376 

255 

235 

127 

102 


C 

0.21 

0.22 

0.33 

0.40 

0.34 

0.59 

SMMR 

B 

0.72 

1.08 

0.42 

0.71 

0.95 

0.47 


S 

1.17 

1.10 

1.15 

0.79 

0.72 

0.76 


N 

677 

361 

227 

226 

126 

96 


C 

0.63 

0.63 


c :.o 

0.50 


XBT 

B 

- 0.03 

-0.04 


0.24 

0.24 



S 

0.84 

0.84 


0.27 

0.27 



N 

158 

155 


8 

8 



C 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

CLIM 

B 

0.03 

0.18 

- 0.14 

0.08 

0.12 

0.01 


S 

0.61 

0.61 

0.55 

0.38 

0.41 

0.35 


N 

729 

376 

255 

235 

127 

102 


i 
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Table 4-3. Comparison of Monthly Mean Anomaly Fields with Ships > 5/Cell, 
March 1982 





Unsmoothed 



Smoothed 




Global 

N. Par.. 

N. Atl. 

Global 

N. Pac. 

N. Atl. 


C 

0.67 

0.52 

0.67 

0.77 

0.58 

0.80 

AVHRR 

B 

-0.36 

-0.50 

-0.29 

-0.44 

-0.54 

-C.30 


S 

0.51 

0.48 

0.42 

0.29 

0.29 

0.21 


N 

795 

434 

267 

368 

210 

153 


C 

0.10 

0.13 

-0.01 

0.40 

0.39 

0.43 

HIRS 

B 

0.30 

0.47 

0.16 

0.29 

0.42 

0.12 


S 

0.92 

0.95 

0.84 

0.4' 

0.41 

0.35 


N 

795 

434 

267 

363 

210 

153 



C 

0.55 

0.30 

0.66 

HIRS 

B 

0.29 

0.36 

0.21 

Version 2 

S 

0.31 

0.34 

0.22 


N 

368 

210 

153 



C 

0.24 

0.37 

-0.05 

0.15 

0.54 

-0.09 

SMMR 

B 

-0.21 

0.05 

-0.76 

-0.17 

0.13 

-0.77 

Night 

S 

1.11 

0.99 

1.19 

0.79 

0.67 

0.69 


N 

690 

392 

213 

300 

200 

95 



C 

0.58 

0.57 

0.75 

0.75 

SMMR/Ship 

B 

0.04 

0.03 

0.07 

0.07 

Composite 

S 

0.47 

0.46 

0.25 

0.25 


N 

438 

394 

207 

203 



C 

0.40 

0.50 

0.79 

0.79 

VAS 

B 

0.90 

0.89 

0.91 

0.91 


S 

0.56 

0.52 

0.26 

0.26 


N 

109 

106 

51 

51 



C 

0.39 

0.38 

0.70 

0.70 

XBT 

B 

-0.27 

-0.29 

-0.47 

-0.47 


S 

0.89 

0.91 

0.35 

0.35 


N 

242 

227 

18 

18 



C 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

CL1M 

B 

0.09 

0.27 

-0.05 

0.13 

0.29 

-0.10 


S 

0.52 

0.48 

0.42 

0.35 

0.32 

0.27 


N 

795 

434 

267 

368 

210 

153 
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Table 4-4. Comparison of Monthly Mean Anomaly Fields with Ships > 5/Cell, 
July 1982 





Unsmoothed 



Smoothed 




Global 

N. Pac. 

N. Atl. 

Global 

N. Pac. 

N. Atl. 


C 

0.62 

0.59 

0.66 

0.70 

0.63 

0.84 

AVHRR 

B 

-0.48 

-0.37 

-0.57 

-0.35 

-0.17 

-0.48 


S 

0.79 

0.93 

0.60 

0.52 

0.62 

0.37 


N 

644 

320 

258 

274 

117 

157 


C 

0.49 

0.43 

0.51 

0.78 

0.52 

0.85 

HIRS 

B 

-0.07 

0.01 

-0.08 

0.09 

0.14 

0.04 

Version 2 

S 

0.69 

0.72 

0.62 

0.38 

0.39 

0.36 


N 

662 

337 

259 

327 

170 

157 

HIRS 

C 

0.52 

0.45 

0.56 

0.79 

0.45 

0.88 

Version 2 

B 

-0.37 

-0.31 

-0.37 

-0.25 

-0.23 

-0.27 

Weighted 

S 

0.69 

0.74 

0.61 

0.38 

0.43 

0.30 


N 

662 

337 

259 

327 

170 

157 


C 

0.46 

0.54 

0.32 

0.55 

0.62 

0.66 

SMMR 

B 

-0.43 

-0.22 

-0.88 

-0.69 

-0.39 

-1.07 

Night 

S 

0.97 

0.87 

0.93 

0.60 

0.48 

0.51 


N 

522 

278 

193 

230 

127 

103 


C 

0.76 

0.74 


0.86 

0.86 


SMMR/ Ship 

B 

-0.04 

-0.06 


-0.10 

-0.10 



S 

0.53 

0.54 


0.27 

0.27 



N 

316 

282 


137 

137 



C 

0.49 


0.46 

0.42 


0.42 

VAS 

B 

0.48 


0.50 

0.55 


0.55 


S 

0.46 


0.45 

0.22 


0.22 


N 

92 


88 

38 


38 


C 

0.58 

0.55 





XBT 

B 

-0.22 

-0.23 






S 

0.94 

0.96 






N 

154 

146 






C 


0.00 

0.00 

0.00 

0.00 

0.00 

CLIM 

B 


0.67 

0.26 

0.70 

0.96 

0.40 


S 


0.73 

0.69 

0.63 

0.50 

0.64 


N 


338 

259 

336 

179 

157 
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One reason the global deviation from climatology Is small la that exten- 
sive areas have deviations very close to zero. On the other hand, large spa- 
tially coherent areas of anomalies > 0.5°C also exist, as evident in Figure 
C-3a showing the ship anomaly field for December 1981. A much better Indica- 
tor of skill is the correlation of the sensor anomaly field with that of ships. 
This correlation is high if the anomalies lie in the right place and are of 
the right sign. The correlations are reduced by random noise, especially if 
the noise is of comparcble magnitude to the signal. In addition, large areas 
of small anomaly will also reduce the correlations because the magnitude of 
the signals in these areas is comparable to noise. Climatology, by defini- 
tion, has zero anomaly correlation and zero skill. 

As seen in Table 4-2, smoothing has reduced the noise and increased the 
correlation coefficients for all cases (except XBT which contains only 8 points 
in the smoothed tield). Part of this is due to random noise reduction in the 
ship field and part due to noise reduction in the field being compared to 
ships. The AVHRR field shows skill both in the standard deviation sense and 
in the correlation sense in the unsmoothed fields, while both the HIRS and 
SMMR fields show poor skill in the unsmoothed fields. It is interesting to 
note the XBTs also appear noisy in the unsmoorhed fields with regard to stan- 
dard deviation, but not correlation coefficient. This may be in part due to 
noise, but it is also due to the fact that the portion of the North Pacific 
Ocean measured by XBTs had a larger anomaly than the North Pacific Ocean as a 
whole. In the smoothed fields, AVHRR shows remarkable skill in both catego- 
ries. In addition, the statistics for HIRS have improved greatly so as to 
show moderate skill in correlations and noise level. It should be noted that 
these statistics are for the unweighted HIRS2 field. The weighted field (not 
shown) provides better visual agreement with the ship field, but statistics 
were not initially computed using the weighted field. The SMMR field statis- 
tics also improved to show some skill in anomaly correlation, but the noise 
levels are high compared to the signal. 

While AVHRR appears extremely good statistically, some aspects of the 
anomaly field, shown in Figure C-3b, are disturbing. In particular, the region 
from about 75°E to I55°E, 5°N to 5°S , shows a coherent warm anomaly of between 
0.5° and 2.5°C, which is not supported by the few ship measurements in the 
area. In fact this area did not enter into the statistics for two reasons: 
first, the scarcity of ships in the area, and secondly because of the "SMMR 
Mask" (see Figure C-3d) which was applied to statistics of all fields so simi- 
lar areas could be compared statistically. This mask is applied to SMMR data 
because accurate SMMR retrievals cannot be obtained less than 600 km from land. 
A second problem is that the AVHRR data has regional biases. Table 4-2 shows 
a cold bias in both the Atlantic and Pacific Oceans, with the Pacific Ocean 
bias being about 0.2° larger. Figure C-3b shows the biases to be largest in 
the high latitudes, which represent the coldest, driest air. On the other 
hand, the western Pacific tropical area previously mentioned, which is the one 
with the highest water vapor content, was spuriously warm. It is possible 
that the algorithm used to correct the AVHRR soundings for water vapor under- 
corrccts the effect for low water vapor content and overcorrects in the case 
of high water vapor content. This is consistent with theoretical studies 
showing that the effects of water vapor content or> brightness temperature grow 
in less than a linear fashion. 
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Statistical results for March 1982 are shown in Table 4-3. This table 
includes results of the reprocessed March 1982 H1RS date, as described in Sec- 
tion B, called HIkS version 2 in the table. Only statistics for the smoothed 
version 2 data were generated. Also included in the table are statistics for 
the field derived using SMMK data and ship data simultaneously, and for the 
VAS retrievals, which were primarily in a small area in the North Atlantic 
Ocean. The SMMR/SH1P composite field was produced only in the Pacific Ocean 
because of calibration difficulties with the SMMR in the Atlantic Ocean. 

Comparison of ships with climatology shows smaller anomalies in March 1982 
than in December 1931. Indeed, the signal in the North Pacific and North 
Atlantic Oceans is on the order of O.J°C. Statistics for AVHRR in March 1982 
are similar to those of December 1981, but the correlations are reduced some- 
what, most likely because the signal is reduced. Statistics for the HIRS2 are 
again similar to those of December 1981. It should be noted that HIRS ver- 
sion 2 is significantly better than HIRS, both with regard to standard devia- 
tion and anomaly correlation. 

Smoothing has reduced the noise and increased the correlation coefficients 
for all cases. Part of t'lis is due to random noise reduction in the ship field 
and part due to noise reduction in the field being compared to ships. The un- 
smoothed AVHRR field in March 1982, shown in Figure C-4b, shows marginal skill 
in the standard deviation sense and good skill in the correlation sense. Sta- 
tistics for the unsmoothed March HIRS version 2 fields were not generated at 
the workshop. It is interesting to note the XBTs also appear noisy in the un- 
smoothed fields with regard to standard deviation and correlation coefficient. 

In the smoothed fields, AVHRR shows increased skill in both categories. 

The statistics for the smoothed HIRS2 fields are almost as good as those of 
AVHHR with regard to standard deviation, but the correlations are considerably 
lower. It should be noted that these statistics are for the unweighted HIRS2 
field shown in Figure C-6. The weighted field, shown in Figure C-6b, provides 
better visual agreement with the ship field, but statistics were not computed 
using the weighted field. The weighted SMMR field is much more noisy, but 
shows some correlation skill in the North Pacific Ocean. An additional field, 
produced using both SMMR data and ship data, shows good agreement in the 
vicinity of ships as expected. Surprisingly, the statistics are only margin- 
ally better than those of AVHRR or HIRS2, which did not have the benefit of 
including the ship data used to verify the fields. VAS shows good results for 
that portion of the North Atlantic Ocean where observations were made but has 
a disturbingly large bias of 0.9°C too warm. 

Figure C-4a indicates that March 1982 had only small areas showing more 
than 0.5° deviation from climatology, especially in the North Atlantic Ocean. 

In the North Pacific Ocean, there are small cold anomalies centered at about 
160°W, 35°N and 175°E, 28°N and a small warm anomaly about 115°W, 20°N. The 
weighted HIRS2/MSU field likewise shows small anomalies in the North Pacific 
Ocean, of the correct sign, centered at the appropriate locations. The anomaly 
centered at 35°N is weaker than in the ship field, and is less apparent in 
the HIUS2 field because It is between 0 and -0.5°, showing up as white in 
the picture, but on the cold side of the 0° contour. The small anomalies in 
the Atlantic Ocean are in almost perfect agreement with those in the ship 
field. The AVHRR field, on the other hand, indicates extensive areas of cold 
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anomaly in boch Che North AClancic and Pacific Ocean, sometimes exceeding 1.5° 
in magnitude. This feature is consistent with the cold bias of 0.3°C to 0.3°C 
indicated in Table 4-3 for AVHRR. The weighted HIRS2 field, like the ship 
field, shows a warm anomaly in the Bay of Bengal, which is absent in the AVHRR 
field, and an area of negligible anomaly along the equator from 75°E to 155°E. 
In contrast, the AVHRR field has a very large spurious warm anomaly in this 
region. Agreement of all fields in the southern hemisphere is reasonable, but 
the AVHRR warm anomalies are larger than those indicated by HIRS2/MSU or the 
ships. The unweighted HIRS2 field is somewhat warmer than the weighted field 
and has essentially no anomalies in the North Pacific Ocean. It is not sur- 
prising that small correlation coefficients were found in the North Pacific 
Ocean for the unweighted HIRS2 version 2 field. 

Statistics for July 1982 are shown in Table 4-4. For this month, statis- 
tics were generated for both the weighted and unweighted HIRS2 version 2 
fields. As observed from looking at the statistics comparing climatology with 
ships, July 1982 was much more anomalous than other months in the workshop 
with a standard deviation of about 0.6°C 3bout a cold bias (Table 4-4 shows 
climatology-ships) of 0.96°C in the North Pacific Ocean and 0.4°C in the North 
Atlantic Ocean. Therefore signals are large and one expects larger correla- 
tions between retrieved anomaly fields and ship anomaly fields. Correlations 
for AVHRR are simlar to those in March 1982, but the standard deviations are 
considerably larger than in the other months. HIRS2 statistics show improved 
correlation over March 1982, as expected, with standard deviations slightly 
degraded over those of March 1982. Statistics for the weighted HIRS2 field, 
relative to the unweighted field, show the weighted soundings have become 
colder on the average by about 0.35°C. Results have improved slightly in 
the Atlantic Ocean and degraded somewhat in the Pacific Ocean. SMMR results 
have improved over previous months, both in the standard deviation sense and 
in the sense of error correlation. The SMMR/ship field, as expected, shows 
good agreement with the ships used to produce the field and shows more of an 
improvement, with regard to standard deviation over the HIRS2 field, which did 
not have the benefit of including ships, than in March. VAS continues to have 
good agreement with regard to standard deviation in the small region of the 
North Atlantic Ocean where colocations with smoothed ships exist. VAS 
retrievals also continue to have a large warm bias, which is about 0.4°C less 
than in March 1982. 

Figures C-7a, b, and c and C-6c show the retrieved anomaly fields deter- 
mined from ships, AVHRR, HIRS, and weighted HIRS data. HIRS again shows very 
good agreement with the ship anomaly fields, which this time are quite large 
and extensive in both the North Atlantic and North Pacific Oceans. The HIRS 
field does have a few areas of spurious cold anomaly in July, located mainly 
within a thin latitude band running from about 15°N to 25°N over most of the 
ocean, and off the west coast of South America from 15°S to about the equator. 
These spurious anomalies, which are mostly of the order 0.5° - 1°C, may be due 
to the effects of aerosols put into the atmosphere by the eruption of El 
Chichon which occurred in Mexico in April 1982. The errors are amplified in 
the weighted field, which is colder in general. The effects of this eruption 
on the determination of sea surface temperatures are much more evident in the 
AVHRR anomaly pattern which shows a large block of spurious cold anomaly from 
10°N to 35°N running across the entire oceanic area, with magnitudes of the 
order of 2°C. 
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As shown in Figures C-6 and C-7, the weighted HIRS2 fields tend to be 
slightly cooler than the unweighted fields. The same result was found for 
December 1981. The weighted fields show better agreement with ship fields in 
December 1981 and March 1982. July 1982 is a special case where spurious cold 
anomalies already existed, possibly due to effects of the El Chichon eruption, 
and results were slightly worse in the weighted fields. 

Weighted fields count clear cases more than cloudy ones and night some- 
what more than day. A simple explanation for weighted fields being colder 
would be the heavier weight placed on nighttime soundings. To test this, 
weighted fields were produced using only nighttime data and using only daytime 
data and compared to the unweighted fields at the appropriate time of day. 

This comparison showed the weighted fields again to be cooler than the 
unweighted ones in both cases, with the effect being larger at night. It 
appears then that sea surface temperature soundings in an area are warmer 
under cloudy conditions than under clear conditions. This implies the effects 
of clouds are somewhat over-corrected in the system by the use of equations 
(4-2) through (4-6). 

E. SUMMARY AND FUTURE PROSPECTS 

Prospects for improved sea surface temperature fields obtained from 
sounders such as HIRS2/MSU appear bright. The current test shows HIRS2/MSU 
fields to be statistically comparable to those of AVHRR in March 1982 and 
superior in July 1982, where the El Chichon eruption had drastic consequences 
on the AVHRR data, but much smaller effects on HIRS2/MSU data. Errors in the 
AVHRR data tend to be systematic and spatially coherent, providing spurious 
warm anomalies in the most humid areas of the tropics and possibly spurious 
cold anomalies at high latitudes. HIRS2 errors tend to be more random in 
nature and resultant fields are improved considerably by averaging more data. 
Our plans in the near future involve reprocessing all data at a roughly 60 km 
spatial resolution, corresponding to 4 x 4 sets of HIRS2 spots broken up into 
2x2 quadrants. We have currently completely vectorized the retrieval code 
for use on the CYBER 205 computer and one-day global retrievals at 60 km reso- 
lution takes only 10 minutes CPU time. We feel this is the finest resolution 
with which retrievals can be done using HIRS2 data because of the need for at 
least four spots in a quadrant to allow for distinct separation of radiances 
into fields of view for the purpose of cloud filtering. These higher resolu- 
tion retrievals should produce much better anomaly fields both because of the 
reduction of random errors by a larger sample, and also because the individual 
retrievals will become more accurate due to increased homogeneity in the 
sounding area both with regard to clouds (assumption of a single cloud forma- 
tion is implicit in the cloud filtering algorithm) and the sea surface temper- 
ature itself. 

Hardware modifications in the near future will also improve results. The 
HIRS2 to be flown starting on NOAA-H (~1987) will have a modification in its 
channel 17, now centered at 2360 cm - *. The current channel monitors emission 
from the upper stratosphere and is highly affected by non-local thermodynamic 
equilibrium. Consequently little use has been made of observations in this 
channel either operationally or experimentally. The new frequency will be at 
2420 cm - * and can be used together with channel 18, at 2500 cm - *, to get 
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Improved sea surface temperatures during the day. The new combination is 
superior to the current channel 18,19 combination for several reasons. The 
new channel 17 is less sensitive to solar radiation than 19; therefore, the 
effect of solar radiation is smaller. Also, channel 17 is insensitive to 
water vapor (atmospheric N 2 is the biggest absorber), while channel 19 has 
moderate sensitivity to water vapor, which sometimes causes problems in tropi- 
cal areas during the day. The spurious cold anomaly in the central Pacific 
Ocean in March, derived from HIRS2*data, in fact appeared only in daytime 
soundings. This problem will not occur with the new channel. Another signif- 
icant improvement in the new channel 17 is that its bandpass is about 50 cm"** 
rather than 200 cm~* in channel 19. This makes the concept of "effective 
channel central frequency," which is used in converting back and forth from 
brightness temperature to radiance, more meaningful. The broader bandpass of 
channel 19 introduces an extra source of noise in the data. 


Perhaps the ultimate in measurement of sea surface temperatures will come 
from an advanced infrared sounder, AMTS, described in Chahine et al. (1984), 
which may fly on polar orbit in the 1990's. This instrument has a number of 
"super-window" channels with atmospheric transmittance *0.98 , even in very 
humid atmospheres, and 10 km spatial resolution with contiguous coverage. 
Simulation studies show soundings with about a 40 km resolution can be per- 
formed in up to 90% cloudiness with instantaneous accuracy of about 0.4°C. 

It is expected that monthly mean fields with accuracies of the order of 0.2°C 
can be obtained at a spatial resolution of about 50 km. This will go a long 
way toward meeting the needs of monitoring sea surface temperatures, at least 
in the climatological sense, and will make possible detailed studies relating 
the effects of sea surface temperatures anomalies in the tropics on atmo- 
spheric circulation. 
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A. SST INTERCOMPARISON WORKSHOP 
1. Introduction 

The results of the SST intercomparison workshop series (JPL 1983) 
are the first examination of monthly mean SSTs derived from MSI data provided 
by the VAS instrumer.t on the GOES series satellites. While VAS instruments 
are currently only on the U.S. geostationary satellites, limiting coverage to 
the western hemisphere, it is hoped that the success of VAS will encourage the 
European Space Agency, Japan, and India to consider installing a VAS instru- 
ment on their future geostationary satellites. Because the procedure to 
derive SSTs from VAS data is still in the developmental stage, several changes 
in the procedure were made between the processing of data for March 1982 and 
processing the data for July 1982. The most significant change was the use of 
the three window channel algorithm (3.9, 11.0, and 12.6 pm) in the 
processing of the July data as opposed to the use of only the two window 
channel (11.0/12.6 pm) algorithm for the March data. Initially only the 
two channel algorithm was used in order to extend the analysis of SST into 
areas of sunglint in the tropics. However, the analysis Of the March data 
showed that little additional data was gained by doing this. In addition, 
further satellite/buoy matches indicated that the triple window channel algo- 
rithm showed a smaller standard deviation than the two window channel algo- 
rithm and was less sensitive to the effects of volcanic aerosol contamination 
and low level inversion conditions. This is due to the smaller brightness 
temperature attenuation by aerosols and water vapor at 3.9 pm than at 11.0 
and 12.6 pm. Thus, the decision was made to use the best product (i.e., 
the three window channel algorithm) for processing the July data. 


2. March 1982 Results 

Two large regions were chosen for analysis of VAS data from GOES- 
East, one in the western North Atlantic and one in the eastern Tropical 
Pacific. Since ship observations of surface layer temperature provide the 
only long-term climatology of SST, Reynolds (1982) climatology has been used 
as a standard from which satellite SST monthly mean anomaly fields were pro- 
duced. Data from all sensors were binned on a two by two degree latitude/ 
longitude grid for each month. SKMR data were required to be more than 600 1cm 
from land in order to avoid contamination from land. Thematic contour charts 
of oensor anomaly fields from climatology for March are shown in Figure C-S 
(see Appendix C). VAS, AVHRR, and ship data all show a pattern of cold to 
warm to cold to warm proceeding southeast off the U.S. east coast; however the 



VAS data have a warm bias of 0.5 to l.0°C. In the South Pacific, the VAS data 
show only a slight warm bias and again are highly correlated with the AVHRR, 
ships, and XBTs. In particular, the VAS and AVHRR thematic contour anomaly 
charts show similar patterns with warm water along the coas* from 20 s to 30*S 
and extending to the west along 30®S, a pool of cold water along the coast 
from 0 to 1G*S, another cold anomaly offshore, and near normal conditions 
elsewhere. The HIRS data show generally weaker anomaly patterns and a warm 
bias near the coastlines due to problems in accurately specifying the land/ 
water boundaries (Susskind, personal communication). The HIRS data do show a 
warm anomaly along 30*S in the eastern South Pacific and a large warm anomaly 
in the western North Atlantic. Little correlation in patterns is found 
betwoen the VAS and the SMMR product. 

Table 5-1 summarizes the cross correlation statistics for each satellite 
verses ship-of-opportunity measurements for March 1982. VAS estimates of SST 
show a warm bias relative to ships for all regions ranging from 40.35 to 
1.73*C. The largest biases (1.73*C and 1.05*0 are found with the lowest 
numbers of matches (21 and 53) and also occur at the largest satellite zenith 
angles (North Pacific region 20-56*N and South Pacific regions 20-56 *S). This 
indicates that the magnitude of the warm bias for the two channel algorithm 
may increase with increasing satellite zenith angle but also suggests that 
noisy ship data may be partly responsible for some of the bias. 

The uniform warm bias in all regions, however, indicates a diurnal sam- 
pling bias and a possible bias in the matches used to tune the empirical algo- 
rithm. Satellite/buoy matches are continuing to be collected in order to 
ensure that a seasonally and geographically diverse set of matches is used to 
update the coefficients for the empirical algorithms. It does appear though 
that the diurnal sampling of VAS data is largely responsible for the warm 
bias. VAS data were generally processed at 1530 and 1830 GMT (1030 and 1330 
LST at the GOES-East subpoint) and only cloud-free observations were used. 
Thus, VAS SSTs might be expected to have a warm bias relative to estimates of 
SST that average day and night data. Diurnal heating of the ocean skin tem- 
perature as observed by satellite infrared data has also been reported by 
Strong (1984) and by Deschamps and Frouin (1984). Future intercomparisons 
must take into account possible diurnal sampling biases of each sensor. 

Additional cross correlation statistics for March show VAS with a scatter 
relative to ships of 0.79-1.24®C. The statistics show VAS well correlated 
with ships, and shows regional correlations very similar to those of the 
AVHRR. The one exception is the far South Pacific region (20-56*S). This 
again is the region of fewest matches and thus should be given little weight. 


3. July 1982 Results 

In the thematic anomaly charts for July (Figure C-8), the effects 
of the El Chichon volcanic aerosol are very evident in the AVHRR data as a 
zonal band of cold anomalies from 10-30®N. VAS data, however, do not show an 
analogous anomaly in those latitudes. This result is due to differences in 
the spectral channels of the VAS and AVHRR, differences in the processing 
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Table 5-1. Cross Correlations of Satellite SST Estimates Versus Ship SST Estimates for March 1982 



Number of 

Matches 


Bias 


Standard Deviation 

Cross Correlation 


AVHRR 

srota 

8 IRS 

VAS 

AVHRR 

SKHR 

sirs 

VAS 

AVHRR 

SNKR 

SIRS 

VAS 

AVHRR 

SMMR 

HIRS 

VAS 

Global 

4322 

1972 

— 

425 

-0.06 

-0.01 

— 

+0.63 

0.81 

1.20 

— 

0.96 

0.58 

0.25 

— 

0.59 

North Pacific 
(0-56*N) 

1563 

815 

— 

127 

-0.26 

-0.05 

— 

+0.52 

0.67 

0.99 

— 

0.92 

0.64 

0.29 

— 

0.61 

North Pacific 
(20-56*N) 

1033 

529 

1054 

53 

-0.39 

-0.01 

+0.54 

+1.05 

0.65 

1.03 

1.07 

0.89 

0.60 

0.36 

0.28 

0.63 

Tropical Pacific 
( 20*N-20*S ) 

837 

412 

858 

165 

40.10 

-0.22 

+0.23 

+0.20 

0.73 

0.87 

0.93 

0.90 

0.13 

0.08 

0.03 

0.10 

South Pacific 
(20-56*S) 

535 

202 

541 

21 

+0.24 

+0.55 

+0.05 

+1.73 

0.78 

1.19 

1.11 

1.24 

0.42 

0.08 

0.28 

-0.39 

South Pacific 
(0-56*S) 

984 

328 

— 

112 

+0.26 

+0.17 

— 

+0.52 

0.80 

1.14 

— 

1.20 

0.29 

0.11 

■ 

0.23 

VAS Pacific Region 
(14*N-30*S) 

178 

81 

— 

181 

+0.06 

-0.16 

— 

+0.35 

0.89 

0.91 

— 

1.06 

0.23 

-0.08 

— 

0.03 

Clobal AVHRR 
El Chichon Mask 

2214 

1088 

2229 

211 

♦0.03 

+0.12 

+0.20 

■rO.55 

0.86 

1.11 

1.04 

1.12 

0.48 

0.29 

0.20 

0.41 

North Atlantic 
(0-56*N) 

715 

315 

— 

186 

-0.33 

-0.92 

— 

♦0.76 

0.61 

1.18 

— 

0.79 

0.58 

-0.02 

— 

0.65 





algorithms, and differences in the average viewing geometry. The VAS, SMMR, 
and ship data all show a warm anomaly in the eastern tropical Pacific. 

In the North Atlantic VAS region, the VAS data appears to be slightly 
warmer than the ship data, but again shows similar patterns. The VAS ari 
AVHRR data show some correlation near the coast of the U.S., but meaningful 
comparisons between the two are hampered by the volcanic aerosol contamination 
in the AVHRR data. The anomaly patterns are much the same in a comparison of 
the VAS/SMMR data, however, the SMMR data is contaminated by "cold" instrument 
warmup noise in much of the North Atlantic (Milman, personal communication). 

In the VAS region of the Pacific, the VAS, SMMR, and ship data all show warm- 
ing. Here, the VAS and SMMR data show a high correlation with a pattern of 
warm anomalies along the coast and extending westward along the equator. In 
contrast, the HIRS data, while not showing any consistent bias in the El 
Chichon region, does show a large cold anomaly in this region. 

Cross correlation statistics for July 1982 are summarized in Table 5-2. 
VAS SSTs again show a slight warm bias in all regions. The very large warm 
bias at large local zenith angles evident in the March 1982 data, however, has 
been eliminated by the use of the three window channel algorithm. Little bias 
is evident in the region of the El Chichon volcanic aerosol (approximately 
10-30*H). In this region, the AVHRR data show a cold bias of 0.50-0.75*0 
relative to ships. The VAS standard deviations are also generally smaller in 
July than in March due to the use of the three window channel algorithm. The 
cross correlations of VAS data with ship data, however, are much weaker in 
July than March. 

After SST Intercomparison Workshop III, additional cross correlation 
tables were generated to try to answer some of the questions raised during the 
workshop. Most important to the interpretation of VAS data was the stratifi- 
cation of AVHRR data into day and night so that the daytime only VAS data 
could be directly compared to daytime only AVHRR data. Although the new cross 
correlation tables are masked to include only data greater than 600 km from 
land (to normalize the comparison between SMMR and the other sensors, but 
greatly reducing the number of VAS/ship matches), some trends are clearly evi- 
dent. In March 1982, AVHRR shows a global average day minus night difference 
relative to ships of +0.43*C. This reduces the VAS minus AVHRR day bias to 
+0.23*C. The VAS verses ship biases remain unchanged since ships measure SST 
at some depth beneath the surface and are relatively insensitive to diurnal 
heating of the ocean skin. In July, on a global basis, the AVHRR day product 
is 0.43*C wanner than ships while the AVHRR night product is 0.72*C colder 
than ships. There is no discemable bias between AVHRR day SSTs and VAS SSTo 
outside the El Chichon zone (i.e., in the South Pacific and North Atlantic), 
while within the El Chichon zone (the mid-Pacific) AVHRR day is 0.69*C colder 
than VAS and 0.50*C colder than ships. These data clearly show that the 
diurnal heating of the ocean skin is being detected by VAS and AVHRR, and 
demonstrates that most of the VAS warm bias relative to the other sensors is 
due to this diurnal variability. 
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Table 5-2. Cross Correlations of Satellite SST Estimates Versus Ship SST Estimates for July 1982 



Number of 

Katches 


Bias 



Standard Deviation 

Crof a Correlation 



AVURR 

SKMR 

SIRS 

VAS 

AVHRR 

SMMR 

HIRS 

VAS 

AVHRR 

SHKR 

HIRS 

VAS 

A7BR& 

smR 

HIRS 

VAS 

Global 

3962 

1826 

— 

437 

-0.54 

-0.18 

— 

+0.60 

0.90 

1.08 

— 

0.85 

0.44 

0.38 

— 

0.26 

North Pacific 
(0-56*N) 

136S 

708 

— 

116 

-0.69 

+0.26 

— 

+0.91 

0.95 

0.89 

— 

0.80 

0.41 

0.46 

— 

0.17 

North Pacific 
(20-56*N) 

514 

221 

— 

26 

-0.18 

+0.11 

— 

+0.40 

0.64 

1.13 

— 

1.25 

0.50 

0.24 

— 

-0.10 

Tropical Pacific 
(2ITN-20*S) 

779 

366 

— 

165 

-0.69 

+0.10 

— 

+0.77 

0.83 

0.83 

— 

1.00 

0.30 

0.27 

— 

0.07 

South Pacific 
C 20—56 *S ) 

958 

480 

— 

51 

-0.54 

-0.31 

— 

+0.61 

0.98 

0.92 

— 

0.50 

0.46 

0.44 

. — 

0.11 

South Pacific 
(0-56*S) 

883 

359 

— 

126 

-0.23 

+0.06 

— 

+0.50 

0.67 

1.03 

— 

1.05 

0.41 

0.26 

— 

0.00 

VAS Pacific Region 
a4-N-30‘S) 

162 

104 

— 

174 

-0.38 

+0.16 

— 

+0.68 

0.81 

1.10 

— 

1.08 

0.25 

0.05 

— 

0.01 

Global AVHRR 
El Chichon Mask 

2305 

1112 

— 

216 

-0.20 

+0.02 

— 

+0.49 

0.77 

1.08 

— 

0.90 

0.49 

0.41 

— 

0.26 

North Atlantic 
(0-56'N) 

695 

324 

— 

195 

-0.81 

-1.06 

—• 

+0.49 

0.87 

1.06 

— 

0.67 

0.38 

0.21 

— 

0.40 







B. EVALUATION OF OTHER PRODUCTS 

I. AVHRR 

The AVERR MCSST is the only operational satellite SST analysis cur- 
rently and is the most accurate and consistent product evaluated at the work- 
shop series. As with all SST data sources, care and understanding must be 
used when evaluating and applying this data. Studies such as that by Legeckis 
and Pichel (1984) are particularly useful in interpreting the weekly MCSST 
analyses. Users must also understand the nature and variability of ocean sur- 
face skin temperature measurements as opposed to ship bulk surface layer mea- 
surements. For example, the MCSST analysis for March 1982 has been criticized 
for showing a warm anomaly along the equator from the western Pacific into the 
western Indian Ocean; an area where ship climatology shows little monthly 
variability. The AVHRR day-night thematic contour analysis (Figure C-9), how- 
ever, shows that this warm anomaly may be the result of diurnal warming of the 
ocean surface. In fact, the AVHRR day-night analyses show a distinct diurnal 
pattern of solar heating from December 1981 to March 1982 to July 1982. In 
December 1981, a consistent zonal band of warm daytime SST anomalies is found 
from about 30-50*S, in the southern (summer) hemisphere. In March 1982, the 
warm anomaly has become more diffuse and shows the largest anomalies on the 
equator. By July 1982, the warm anomaly evident as a zonal band in the north- 
ern (summer) hemisphere. Diurnal variability of tue oceans skin is being mea- 
sured by satellite sensors, as is evident from the analysis of AVHRR day-night 
measurements . 


2. SMMR 

The problems with the SMMR SST product are largely due to instru- 
mental difficulties. The SMMR calibration biases are large and vary in time 
and space, and sidelobe interference requires observations to be greater than 
600 km from land. In spite of these difficulties, SMMR analyses of the 
Pacific and Indian Oceans appear reasonable. Unfortunately, the calibration 
problem makes it difficult to evaluate the problem of microwave emissivity 
changes of the ocean surface with wind speed, while the land mask restricts 
analysis of the important boundary currents. The SMMR/ship product is an 
improvement on SMMR alone, but it does not take full advantage of all the 
different sensors for measuring SST. 


3. HIRS/MSU 

Evaluation of the HIRS/MSU product is difficult because of changes 
in the product from one time to the next and because the data were presented 
late. The HIRS/MSU anomaly patterns generally look noisy and weaker than the 
anomaly patterns of the other sensors. In March 1982, the HIRS/MSU shows no 
correlation with any of the other products and a standard deviation from cli- 
matology of about 1°C. The July 1982 statistics are better, but the anomaly 
patterns are inconsistent, showing an overall cool bias. Particularly trou- 
blesome is a cool anomaly in the eastern Equatorial Pacific where all the 
other sensors show a warm anomaly. 
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C. RECOMMENDATIONS 
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l. Improvements in Infrared Sensors 

Recent theoretical and empirical studies of the infrared portion of 
the earth's spectrum have revealed that neither VAS nor AVHRR have the optimal 
channel selection for SST detection. Studies are now underway to determine 
which window regions using a filtered radiometer would yield the most accurate 
SSTs. In the long term, though, an infrared spectrometer interferometer 
instead of a filtered radiometer will be a much better instrument since, it 
would permit use of all portions of the infrared window regions to be utilized. 


2. A Combined Product 

Efforts should begin on a combined satellite SST product that takes 
advantage of the benefits of each sensing system discussed in the workshop 
series. Such an approach should use the raw data from each instrument, not 
just the finished products such as the SMMR/Ship composite. The McIDAS system 
has the capability of processing raw data from all sensors used in the work- 
shop series. It is time to begin a program to produce an operational SST 
analysis. 


3. Research Panel on SST Sensing 

A research panel to set research program goals, evaluate present 
systems, and recommend areas for further study should be set up under the 
direction of NSF or other appropriate agency. This panel should coordinate 
efforts between ongoing ocean research programs and the remote sensing commu- 
nity. This panel could also serve as the focus for the development of a com- 
bined SST product. 


t 

? 







: 

* 


i * 



/ 


v. 


N86-16860 

SECTION VI 

BLENDED SST FIELDS USING COMBINED Sf-MR/SHIP DATA 
(SUMMARY OF TECHNIQUE) 

T. Wilheit and D. Han 

NASA/Coddard Space Flight Center 
Grecnbelt, MD 


A. INTRODUCTION 

All remote sensors have bias problems which vary in space and time, albeit 
rather slowly. The SMMR can in this regard be considered a pathological 
example of a more general truth. Even ship measurements have a token bias of 
the order of a few tenths of a degree caused by heating of the thermistor by 
the hot engine room environment. This ship bias can be neglected if one dealt 
in SST anomaly (departure from climatology) since ships with a similar bias 
are the principal source of the data used to generate the climatology. The 
ships, however, are very poorly distributed globally, and individually 
represent very noisy measurements. Each remote sensor has its own distri- 
bution of data as well. The infrared measurements are very sparse in cloudy 
areas; the SMMR has coverage limitations due to the interference of land areas 
and only half-time operation. If one simply averaged the data from these 
various sources, the combination of the biases and varying sampling would 
introduce artifacts into the analyzed product which would reflect the sampling 
and which could be quite misleading. What is needed is an analysis scheme 
which uses the strength of each type of measurement to compensate for the 
weakness of the other. In particular, the remote sensors provide reasonable 
estimates of gradients of SST over length scales greater than the sensor 
resolution but smaller than hemispheric. Ships can be considered unbiased but 
are poorly distributed on small scales except in the densest shipping lanes. 

An analysis scheme based on a suggestion by Boll (1981) has been developed 
which exploits just this approximately complementarity to provide an improved 
SST analysis. The resulting accuracy appears to be in the required ~0.5*C 
range for the particular case of using SMMR and ship data to produce two 
degree latitude by two degrees longitude monthly analyses. Proper use of 
other data sources could only improve this product. Moreover, it could easily 
be argued that the mathematics used in this analysis technique are not particu- 
larly sophisticated or even rigorous and a more precise treatment may be of 
benefit. However the technique is very precisely targeted to the problems 
actually observed in the measurements, and extreme care must be taken, in any 
attempt to improve the mathematical foundation of the technique, not to lose 
sight of the measurement realities with which the technique is designed to 
deal. It could also be argued that biases of the sort observed here are 
inherent to all currently available remote sensing approaches to the measure- 
ment of any geophysical variable, and that the SST case studied here is only 
one of the simpler examples of a more general problem. 


B. DESCRIPTION 01 THE TECHNIQUE 

Consider an analysis algorithm as depicted in Figure 6-1. This analysis 
algorithm accepts measurement values with associated uncertainties at the grid 
points of the field. It also accepts values for the first and second 
derivatives of the field in finite difference fora, again with associated 
uncertainties. If there are enough values of the field and its derivative 
specified, then the field is overdeterained and the analysis algorithm derives 
a minimum weighted square error solution as the resultant field. An 
associated uncertainty field is a by-product of this computation. An analysis 
technique of this form is the basic building block of the Bias Removal 
Analysis Technique (BRAT). The specific finite differences used are: 

Differences of adjacent cells row-wise: 

R(i,j) « F(i+1, j) - f(i,j) 

Differences of adjacent cells column-wise: 

C(i, j) - F(i, j+1) - F(i, j) 

Differences of adjacent cells diagonally (both ways): 

Dl(i,j) - F(i+1, j+1) - F(i, j) 

D2(i,j) - F(i+1, j-1) - F(i, j) 

Two-space differences both column and row-wise: 

R2(i,j) = F(i+2, j) - F(i, j) 

C2(i, j) » F(i, j+2) - F(i, j) 

Finally, the Lambertian (a second difference operator): 

L(i,j) = 4F(i,j) - F(i+1, j) - F(i-1, j) - F(i,j+1) - F(i,j-1) 

The weight applied to any observation is the inverse square of its 
presumed uncertainty. For present purposes one degree celcius is U6ed 


GRADIENT 
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(ASSOCIATED 

UNCERTAINTY) 


BLACK 

BOX 


RESULTANT FIELD 

(ASSOCIATED 

UNCERTAINTY) 


ABSOLUTE 

MEASUREMENTS 

(ASSOCIATED 

UNCERTAINTIES) 


Figure 6-1. Generalized Analysis Algorithm uTiich is the Basic Building 
Block of the Bias Removal Analysis Technique. 
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as the uncertainty for individual SMMR and ship measurements. The weight of 
any of the differences is similarly the inverse square of the uncertainty 
which is derived directly from the uncertainties of the associated points in 
the field from which the finite differences are derived. The resulting 
problem is to solve for a set of values for F* (i • j ) which minimises the 
function: 


E - 


i.j 


[F(i,,Q - rHl.J )! 2 + [A(i,j) - A^i.j )) 2 + 


w F (i,j) 


w A U,j) 


+ [B(1 ♦ 

W B (i,J> 


w L U,j) 


where the unprimed quantities are given and the primed 

denotes summation over all grid points, and the W(i,j) 
weights. Traditionally, this equation is reduced to a 
results in a formal solution: 


are solved for. 



represent data 
matrix form which 


S - X F 


F o r 1 S 


If F is not particularly large the equation can actually be solved this way. 
However, in a practical case the field could be defined on a 50 x SO array; 
this would result in the F vector in the above equations having a length of 
2,500 and the matrix A would be 2,500 x 2,500. The direct inversion of such a 
matrix would hardly be a practical matter. One alternative is to solve for F 
by minimizing the error summation with respect to each element of F indivi- 
dually and to iterate through all the locations until satisfactory convergence 
is obtained. It can readily be shown that such a procedure is convergent (in 
a mathematical sense). It was however found that the convergence in the 
present application was unreasonably slow. Fortunately, a compromise between 
the two extremes of solving for all elements of F or a single element at a 
time is possible. It was found that solving for 5 elements at a time provided 
reasonable convergence and the inversion of a 5 x 5 matrix is a simple matter 
on any modern computer. 

Figure 6-2 illustrates the use of this basic building block in a typical 
application of BRAT. In this case we have treated the SMMR day and night data 
as though they cone from different sensors since the bias and drift problems are 
independent on the day and night portions of the orbit. In principle, other 
sensor types could be included at this level; each sensor should be self- 
consistent as far as bias and drift are concerned so this same day/night 
separation would seem appropriate for any current space-borne remote sensor. 

At this level, no derivative Information is available so climatology is used 
with an extremely small weight ( 10 ~ 6 ) simply to prevent, numerical problems. 
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Figure 6-2. Application of the Basic Building Block in the Bins 
Removal Analysis Technique 


Each application of the building block of this level results in a .sea 
surface temperature anomaly field with an associated uncertainty. Each of 
these resultant SST fields is differentiated and the veighta of the 
derivatives are determined from the weights of the elements of the SST field. 
The derivative fields are averaged with their associated weights. The 
weighted average of the resulting set of derivative fields is used as input to 
the final application of the building block along with the ship d/.ta. In this 
manner any constant bias in the £MMR or other remotely censed data is 
automatically removed, and any slowly changing bias is overpowered in the 
weighted least square error solution as long as there are enough ship data. 
Conversely, random noise in the ship data is smoothed out by the derivative 
constraints provided by the SMMR data. However, any bias in the ship data is 
passed through the analysis unmodified. 

In order to test this bias removal analysis technique (BRAT), we have 
chosen an area in the Pacific Ocean from 50°S to 50*N and from 13°E to 80 # W. 

We have produced monthly analyses of SST anomaly (with respect to the 
Robinson-Bauer climatology (Reynolds 1983) with a spatial grid size of 2 x 2 
degrees. We have used the day and night observations of the SMMR as different 
relative sensors ana the file of ship observations available from NOAA as the 
absolute sensors. We have used all beam positions in the SMMR night data but 
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a restricted set of day data to reduce the impact of ionospheric Faraday 
rotation on the resulting analysis. This will provide an extremely demanding 
test of the ability of the analysis to reduce large scale biases since the 
various SST workshops showed that in this respect the SMMR was the most 
severely afflicted of the remote censing techniques examined. 

Figure C-lOc (Appendix C) shows the results of application of BRAT to the 
data from July 1980. One can easily sec patterns of anomalies which look 
reasonable. They rarely exceed 2*C in magnitude and are spatially coherent. 

An exception one might be inclined to question is the cold anomaly extending 
across the South Pacific which exceeds 2.5*C. Because of the interference of 
large land areas with the S MMR retrievals, there are large data gaps around 
the continents. The data gap in the Micronesia area suggests excessive 
conservatism in protecting SMMR retrievals against this interference. 

Figures C-lOa and C-lOb show analyses for the same period based on 
ship-only and SKMR-only data respectively. Note that the ship data have large 
gaps particularly in the Southern Hemisphere. Note also that even where the 
ship data are dense the SSTs have a great deal of uncorrelated structure 
suggesting a noisy product. It is particularly worthy of note that the 
suspicious feature in the South Pacific discussed in the previous paragraph 
is, in fact, supported in the Ghip data although there is not nearly enough 
data to delineate its form. 

On the other hand, the SKMR-only product is quite smooth and delineates 
the anomaly pattern of figure C-lOa quite well. However there is a noticeable 
drift in the absolute value which is particularly noticeable in the North 
Pacific. This io precisely the sort of error that the bios removal analysis 
technique is designed to reduce. Quantitative estimates of the quality of 
these analysis products are difficult. However, estimates have been attempted 
and suggests an rms error of the order of 0.5*C. 
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INTERCOMPARISON OF GLOBAL SST FIELDS DERIVED FROM 
SATELLITE SENSORS AND SHI? OBSERVATIONS 

S. E. Pazan 

Scrippa Institution of Oceanography 
La Jolla, California 


A. INTRODUCTION 

Near global comparisons have been made between binned forms of SST 
obtained from different satellite instruments, ship meteorological observa- 
tions, and the output of a meteorological model. Binning hereinafter means 
gridding by simple summation of data values in a bin, which in this study was 
a 2* x 2* box or a 10* x 10* box surrounding a grid node. In the binned com- 
parisons, the question is asked: are the binned data from two instruments 

drawn from the fame population? The first two measures of statistical moments 
of the populations are used to test the hypothesis that they are. 

Comparisons were also attempted between raw forms of the various SST 
retrievals, using a structure function analysis. The results of this analysis 
were relatively uninteresting and are therefore omitted for the sake of 
brevity. 


B. DATA SOURCES 

Five different types of SST field estimates were available for intercom- 
parison. They were: 1. Multi-channel sea surface temperatures (MCSST) from 

the Advanced Very High Resolution Radiometer (AVHRR); 2. the Scanning Multi- 
channel Microwave Radiometer (SMMR); 3. the High-Resolution Infrared Sounder 
(HIRS); A. ship engine intake temperatures collected for marine weather ana- 
lyses; 5. Fleet Numerical Oceanography Center (FNOC) sea surface temperature 
monthly mean analyses. The instruments, data collection procedures, and error 
budgets of the respective SST analyses have been discussed elsewhere in the 
NASA/JPL SST workshop reports (JPL 1983 and 19 w<»j, and will not be discussed 
further here. 


C. ANALYSIS TECHNIQUES 

The "pool-permutation procedure" (PPP) described by Preisendorfer and 
Barnett (1983) has been used to determine whether the first and second moments 
of SST anomaly fields produced by the different instruments are sta- 
tistically identical. The PPP measure of the difference between first moments 
Is called "SITES", and the measure of the difference between second moments is 
called "SPRED." The PPP method allows a rigorous, nonparamctric test of the 
contention that the satellite-derived SST fields are similar to each other 
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and/or a SST field derived from conventional ship data. The power of the 
method (detection level) is estimable. Further, the technique is robust for 
relatively limited data sets (only 4 realizations in our case). 

The PPP test statistics values ("SITES" and "SPRED") fall into eight 
categories; the hypothesis of statistical identity is considered to be shown 
to be false only if the statistic falls in the last category, the category 
with the largest SITES or SPRED statistic value. If the PPP test were applied 
to a completely random set of data fields, any randomly selected permutation 
would fall into the last category 1 out of every S times; in this sense, the 
test has a confidence level of 7 out of 8, or 881. 

For further discussion of this method, the reader is referred to 
Preisendorfer and Barnett (1983). 


D. RESULTS 

The PPP tests were applied to the data sets discussed in Section B. As 
discussed above, the PPP test statistic values fall into 8 unique values, 
called "categories." 


1. Sensitivity Test 

Before drawing any conclusions from the results of PPP intercom- 
parisons it is useful to have some idea of the power of the test, in particu- 
lar testing the validity of applying this statistical test to a time series 
only four months long. Tests of the PPP were done using modified versions of 
the ship data sets, i.e., four month 'test data* sets were constructed using a 
transformation. T' « F(p,o) + T, where F represents a random population with a 
mean of p and a standard deviation of o. Several realizations of the PPP test 
were made for each pair of p and a values Table 7-1. 

For p not equal to zero and large o, the SITES test poorly discriminates 
between the T and T' fields; that is, the ability of this test to discriminate 
between fields actually degrades with increasing o. For instance, given a 
p •> 0.05, the SITES test was able to discriminate between T and T* for a o of 
0.06, but it was not able to discriminate for a a of 0.35. On the other hand, 
the SPRED test discriminates well between data sets differing by a random 
noise element of the mean bias is small. As the mean bias increases, this 
test is less capable of discriminating. Thus, when p =« 0 statistical identity 
iu rejected for o > 0.35 both in the northern and southern regions. However, 
if p is increased to 0.2, a must be 0.46 before statistical identity is uni- 
formly rejected. For a plausible bias, p, and random noise, o, the SPRED test 
also seems quite sensitive and complements the SITES test well. 


a. Comparisons of Global Satellite Data Fields . Results of the 
global SST comparisons are shown in Tables 7-2, 7-3, and 7-4. Since other 
results have indicated that AVHRR may have problems with water vapor contami- 



TABLE 7-1. Results of PPP Intercomparisons of Ship Injection Temperatures T 
With an Artifical Data Set T* , Where Each Temperature T*(*) * F(y,c) + T(*). 

F is a Random Function With a Mean cf y and a Standard Deviation of o. These 
Intercomparisons Were Run 20 Times for 20 Independent Random Realizations and 
the Number of Times a PPP Test Rejected the Statistical Identity of T* and T 
is Tabulated vs. the Number of Times a PPP Test Accepted the Statistical Iden- 
tity of T' and T. 





SITES Tests: 




Northern Region (30°N - 

50 # N) 





Number of 

Rejections/Number of 

Acceptances 




V* 

o = .06 

0.12 

0.17 

0.23 

0.35 

0.46 

0.00 

6/14 

1/19 

1/19 

6/14 

1/19 

4/16 

0.05 

20/ 0 

20/ 0 

20/ 0 

19/ 1 

13/ 7 

9/11 

0.10 

20/ 0 

20/ 0 

20/ 0 

20/ 0 

20/ 0 

19/ 1 

Southern Region (50*S - 

30 # S) 





Number of 

Rejections/Number of 

Acceptances 




W 

o = .06 

0.12 

0.17 

0.23 

0.35 

0.46 

0.00 

3/17 

5/15 

3/17 

3/17 

5/15 

3/17 

0.05 

20/ 0 

20/ 0 

17/ 3 

11/ 9 

3/17 

6/14 

0.10 

20/ 0 

20/ 0 

20/ 0 

20/ 0 

17/ 3 

14/ 6 




SPRED Tests: 



Northern Region (30*N - 

50'N) 





Number of 

Rejections/Number of 

Acceptances 





o = .06 

0.12 

0.17 

0.23 

0.35 

0.46 

0.00 

12/ 8 

19/ 1 

20/ 0 

20/ 0 

20/ 0 

20/ 0 

0.05 

0/20 

2/18 

15/ 5 

19/ 1 

20/ 0 

20/ 0 

0.10 

0/20 

0/20 

0/20 

13/ 7 

20/ 0 

20/ 0 

0.20 

0/20 

0/20 

0/20 

0/20 

17/ 3 

20/ 0 

0.40 

0/20 

0/20 

0/20 

0/20 

1/19 

18/ 2 

0.60 

0/20 

0/20 

0/20 

0/20 

0/20 

7/13 

Southern Region (50°S - 

30’S) 





Number of 

Rejections/Number of 

Acceptances 




y 

o = .06 

0.12 

0.17 

0.23 

0.35 

0.46 

0.00 

3/17 

7/13 

15/ 5 

18/ 2 

20/ 0 

20/ 0 

0.05 

0/20 

2/18 

12/ 8 

17/ 3 

20/ 0 

20/ 0 

0.10 

0/20 

0/20 

2/18 

13/ 7 

20/ 0 

20/ 0 

0.20 

0/20 

0/20 

0/20 

1/19 

14/ 6 

20/ 0 

0.40 

0/20 

0/20 

0/20 

0/20 

5/15 

19/ 1 

0.60 

0/20 

0/20 

0/20 

0/20 

0/20 

12/ 8 
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TABLE 7-2 


. Number of Space Grid Points. Below and Left of the Data Matrix 
Diagonal, the Number of Space Points in the Data Swarms Used in 
the PPP Tests are Printed. Row and Column Numbers Indicate These 
Instruments: 1) AVKRR, 2) SMMR, 3) HIRS, 4) Ship, and 5) FNOC. 


a) North Temperate 

(30*N-60*N, 

0*-360*E) 





1 

2 

3 

4 

AVHRR 

1 





SMMR 

2 

223 




HIRS 

3 

1140 

234 



Ship 

4 

984 

233 

1040 


FNOC 

5 

517 

235 

535 

525 


b) Global Tropics 

(30*S-30*N, 

0*-360*E) 





1 

2 

3 

4 

AVHRR 

1 





SMMR 

2 

1250 




HIRS 

3 

3210 

1430 



Ship 

4 

1110 

289 

1360 


FNOC 

5 

936 

824 

1150 

470 

c) South Temperate 

• (bO^S-SO'S, 

0 # -360*E) 





1 

2 

3 

4 

AVHRR 

1 





SMMR 

2 

470 




HIRS 

3 

2080 

470 



Ship 

4 

192 

5 

192 


FNOC 

5 






7-4 


TABLE 7-3. Pool Permutation Procedure, SITES Statistics on 2* Binned Data. 
Test the Hypothesis That Two Datasets are Drawn From the Sana Population. 

Below and Left of the First Data Matrix Diagonal, the Relative Category of the 
SITES Comparison Between the Observed Data Swarms is Printed; the PPP Test 
Statistics Fall Into Eight Unique Values, Called "Categories". Above the 
Right of the Data Matrix Diagonal, T(rue)/F(alse) is Printed, Indicating 
Whether or Not the Hypothesis is Supported or Rejected, Respectively. Below 
and Left of the Second Data Zlatrix Diagonal, the Percentage of the Artifi- 
cially Constructed Data Swarm Test Values Which Fall Below the Observed 
Swarms' Test Value is Printed; Above and to the Right of the Second Data 
Matrix Diagonal, the SITES Test Value for the Observed Data Swarms is Printed. 

**SITES Statistics** 


a) North Temperate 

(30*N-60*N, 

, 0 # -360*E) 



• 




Category /Hypothesis 


Percent 

/Test Value 





1 

2 

3 

4 

5 


1 

2 

3 

4 

5 

AVHRR 

1 


F 

F 

T 

F 

1 


1.5 

1.0 

0.32 

11 

SMMR 

2 

8 


F 

F 

F 

2 

88 


0.82 

1.2 

6.3 

HIRS 

3 

8 

8 


F 

F 

3 

88 

88 


1.0 

9.5 

Ship 

4 

6 

8 

8 


F 

.4 

59 

88 

88 


12 

FNOC 

5 

8 

8 

8 

8 


5 

88 

88 

88 

83 


b) Global Tropica 

(30*S-30*N, 

0*-360*E) 







Category/Hypothesis 


Percent/Test Value 





1 

2 

3 

4 

5 


1 

2 

3 

4 

5 

AVHRR 

1 


F 

F 

T 

F 

1 


0.86 

0.89 

0.63 

3.7 

SMMR 

2 

8 


F 

F 

F 

2 

82 


1.3 

0.98 

3.3 

HIRS 

3 

8 

8 


T 

F 

3 

82 

82 


1.0 

2.7 

Ship 

4 

6 

8 

7 


F 

4 

67 

82 

72 


3.4 

FNOC 

5 

8 

8 

8 

8 


5 

82 

82 

82 

82 


c) South Temperate 

(60°S-3C 

> # s. 

0*-360*E) 







Category/Hypothesis 


Percent/Test Value 





1 

2 

3 

4 



1 

2 

3 

4 


AVHRR 

1 


F 

F 

F 


1 


1.1 

1.7 

0.60 


SMMR 

2 

8 


F 

T 


2 

83 


1.3 

1.1 


HIRS 

3 

8 

8 


F 


3 

83 

83 


1.2 


Ship 

4 

8 

5 

8 



4 

83 

46 

83 
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Table 7-4. SITES Statistics on 10* Binned Data. (See Table 7-3 for 
Explanation of the Organization of This Table.) 

**SITES Statistics** 


a) North Temperate 

(30*N-60*N, O' 

*-360*E) 







Category/Hypothesis Test 

Percent/Test Value 





1 

2 

3 

4 

5 


1 

2 

3 

4 

5 

AVHRR 

1 


F 

F 

T 

F 

1 


2.4 

1.0 

0.17 

13 

SMMR 

2 

8 


F 

F 

F 

2 

88 


1.1 

1.8 

6.8 

HIRS 

3 

8 

8 


T 

F 

3 

88 

88 


1.2 

15 

Ship 

4 

2 

8 

7 


F 

4 

10 

88 

73 


17 

FNOC 

5 

8 

8 

8 

8 

1 

5 

88 

88 

88 

88 


b) Global 

Tropics (30*S-30*N, 0*- 

•360*E) 







Category/Hypothesis Test 

Percent/Test Value 





1 

2 

3 

4 

5 


1 

2 

3 

4 

5 

AVHRR 

1 


F 

F 

T 

F 

1 

82 

1.1 

1.2 

0.57 

4.0 

SMMR 

2 

8 


F 

F 

F 

2 

82 


1.8 

1.2 

5.1 

HIRS 

3 

8 

8 


T 

F 

3 

82 

82 


1.1 

4.0 

Ship 

4 

7 

7 

7 


F 

4 

67 

82 

72 


4.6 

FNOC 

5 

8 

8 

8 

8 

1 

5 

82 

82 

82 

82 


c) South Temperate 

(60*S-30*S, 0* 

-360*E) 







Category /Hypothesis 

Test 

Percent/Test Value 





1 

2 

3 

4 



1 

2 

3 

4 


AVHRR 

1 


F 

F 

T 


1 


1.3 

2.7 

0.59 


SMMR 

2 

8 


F 

T 


2 

83 


2.7 

1.2 


HIRS 

3 

8 

8 


T 


3 

83 

83 


1.5 


Ship 

4 

7 

5 

7 



4 

73 

46 

63 
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nation in the Tropics, the tests were taken on data swarms for three areas: 
the first north of 3C*N; the second between 30*S and 30*N; and the third south 
of 30*S. These teats are displayed in a compact matrix form, in which rows 
and columns each represent data swarms for particular instruments. The column 
and row order are always the some: AVHRR, SMMR (nighttime), HIRS, ship, and 

FNOC. For example, the entry for the first row and second column (and vice- 
versa) is an AVHRR-SMMR test result. Different results are shown above and 
below the matrix diagonal, in order to show results in a compact form. 

Figure 7-1 shows the spatial pattern of grid points used in each PPP com- 
parison. Since the time and space distribution for each data swarm differs, 
the spatial pattern of coincident grid points differs for each sensor combina- 
tion. The maximum number of space points is found in the AVHRR vs. RIRS com- 
parison. Because of the coastline mask on SKMR data, comparisons with SMMR 
have fewer space points. The numbor of observations in the tropics 
(Table 7-2(b)) tends to be somewhat larger than in the northern regions; SMMR 
comparisons no longer have the fewest spatial points, reflecting the fact that 
relatively less ocean lies adjacent to coastland in the tropics; therefore the 
SMMR land mask is not as important. The southern region (Table 7-2(c)) gen- 
erally has fewer data points in tho comparisons; its FNOC analysis is absent 
couth of the equator, while ship data is much sparser. SMMR comparisons are 
sparse in the southern hemisphere also. 


b. 2* SITES Tests . Results of the 2 * binned SITED test for the 
northern region (Table 7-3(a)) indicate that the data swarma are not drawn 
from the same population for any comparisons except the AVHRR - ship compari- 
son. The test category for the AVHRR - ship SITES comparison is 6, indicating 
an imperfect agreement. The SITES test values are all greater than 0.8 except 
for the AVHRR - ship comparison. 

Results of the 2* binned CITES test for the tropics (Table 7-3(b)) indi- 
cate that ship and AVHRR SST and ship and HIRS are statistically identical. 
Interestingly enough, the hypothesis of statistical identity was rejected for 
the AVHRR vs. HIRS. Since the El Chichon volcanic eruption contaminated AVHRR 
dat for July, 1982, this result demands a closer examination. Figure 7-l(c) 
shows the data distribution of ships vs. AVHRR and indicates very few spatial 
points were present in the tropics. This is probably due to the paucity of 
ship data. On the other hand, a color plate (Figure C— 7 ) of AVHRR SST - ship 
SST indicates that some of the areas of greatest difference in this month are 
not present in the joint distribution map, apparently because cne of the other 
months was deficient in either ship or AVHRR data in those regions. Also, 
'•uuth of the equator, the thematic map indicates that the SST differences are 
very small. The consequence of this is that over the four months the PPP 
tests are unable to discriminate effectively between tropical ship SST and 
AVHRR SST, in spite of widespread El Chichon contamination. 


C. 10* SITES Tests . SITES statistics for 10* binned data are 
shown in Table 7-4. In the northern and tropical regions the results are 
identical to 2* binned data SITES results. The southern region differs for 
JO* bins in that statistical identity is accepted between AVHRR and ships and 
HIRS and ships. 
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Figure 7-1. Spatial distributions of grid points used in PPP comparisons: 
a) AVHRR vs. SMMR; b) AVHRR vs. HIRS; c) AVHRR vs. ship; 
d) AVHRR vs. FNOC; e) SMMR vs. HIRS; f) SMMR vs. ship; 
g) SMMR vs. FNOC; h) HIRS vs. ship; i) HIRS vs. FNOC; j) ship 
vs. FNOC; k) ship vs. ship. 




7-8 



lATItuOf LATITVOt 



ORIGINAL PAGE IS 
OF POOR QUALITY 





Of 


Figure 7-1. (Continued) 


7-9 



d. I* SPRED Teats . Results of the SPRED test, i.e., the equal- 
ity of the second moments, for the northern region (Table 7-5(a)), indicate 
that the data sets are not drawn from the same statistical population for all 
SMMR comparisons except SMMR vs. AVHRR, and reject the statistical identity of 
ships and AVHRR. On the other hand, this test indicates statistical identity 
between any pair drawn from the AVHRR, ship, or HIRS data swarms. In the 
tropics, the SPREE test (Table 7-5(b)> rejects statistical identity only for 
HIRS vs. ship. The test otherwise accepts statistical Identity for every com- 
parison test. In the southern regions, the test (Table 7—5 ( c ) ) rejects sta- 
tistical identity in AVHRR vs. SMMR and HIRS vs. SMMR. Ship vs. SMMR indi- 
cates statistical identity, but with very few spatial points as mentioned 
before. 

However, as has been concluded previously, where SITES comparisons indi- 
cate a possible bias, SPRED test results indicating identity should be sus- 
pect. This would invalidate many of the positive SPRED comparisons. Con- 
versely, where SITES comparisons are good, very small SPRED values indicate no 
identity by the SPRED test. For this reason, ship vs. AVHRR comparisons, 
indicating statistical identity, are suspect. 


e. 10* SPRED Tests . Table 7-6 shows SPRED results for 10* 

binned data; the results are essentially the same as the 2* bin SPRED test. 


E. SUMMARY 

The SUES test accepts statistical Identity between AVHRR and ship data 
swarms. The SPRED test is consistent with these results, in the following 
sense. The presence or absence of a large SITES test value can bias the SPRED 
test towards acceptance of statistical identity, or rejection of statistical 
identity, respectively. Therefore, th : rejection of statistical identity in 
the SPRED tests between AVHRR and ships and the North temperate regions is 
probably due to the enhanced sensitivity of the SPRED test when SITES values 
are small. Contrariwise, acceptance of statistical identity in the SPRED test 
between SMMR and FNOC in the tropics, for instance, is made suspect by the 
clear rejection of statistical identity in the SITES tests between that par- 
ticular pair. 

HIRS and ships show good agreement in both SPRED and SITES at times, but 
this is more inconsistent than the PPP test results of AVHRR and ship 
comparisons. 

The general rejection of statistical identity in comparisons between 
satellite sensors themselves is noteworthy. Not even AVHRR and HIRS were 
statistically alike by the PPP tests. This seems to imply that the level of 
noise in AVHRR and HIRS was larger than that in the binned ship data. It 
implies that although at least two satellite sensors approach ship SST by the 
statistical measures used in this study, all of the binned satellite SSTs were 
noisier than the binned in situ ship SST. SMMR differs by particularly large 
PPP statistic values from other satellite sensors, and the FNOC SST fields 
show the least agreement with other sensors. 


Table 7-5. Pool Permutation Procedure, SPRED Statistics on 2* Binned Data. 

Test the Hypothesis That Two Datasets are Drawn From the Same 
Population. Below and Left of the First Data Matrix Diagonal, 
the Relative Category of the SPRED Comparison Between the Ob- 
served Data Swarms is Printed; Above and Right of the Data Matrix 
diagonal, T(rue)/F(alse) is Printed, Indicating Whether or Not 
the Hypothesis is Supported or Rejected, Respectively. Below 
and Left of the Second Data Matrix Diagonal, the Percentage of 
the ArtificiallyConstructed Data Swarm Test Values Which Fall 
Below the Observed Swarms' Test Value is Printed; Above and to 
the Right of the Second Data Matrix Diagonal, the SPRED Test 
Value for the Observed Data Swarms is Printed. 


**SPRED Statistics** 


(a) North Temperate (30*N-60 # N, 0*-360*E) 


Category/ 

Hypothesis Percent/Tet : Value 




1 

2 3 4 5 


1 

2 

3 A 

5 

AVHRP. 

1 


T T F T 

1 


0.1A 

0.66 x 10~5 0.066 

0.22 x 10-3 

SMMR 

2 

6 

F F F 

2 

63 


0.13 0.A5 

0.A7 

HIRS 

3 

1 

8 T T 

3 

1.0 

88 

0.055 

0.36 x 10-5 

Ship 

A 

8 

8 7 T 

A 

88 

88 

73 

0.037 

FNOC 

5 

1 

8 15 

5 

1.0 

88 

1.0 A9 


(b) Global 

Tropica 

(30*S- 

-30 *N, 

0*-360®E) 





Category/ 

Hypothesis 



Percent/Test Value 




1 2 3 A 5 


1 

2 3 

A 

5 

AVHRR 

1 T T T T 

l 


0.0087 0.57 x 10-** 

0.078 

0.12 

SMMR 

2 3 T T T 

2 

32 

0.0018 

0.017 

O.0A7 

HIRS 

3 ?. 1 FT 

3 

1.0 

1.0 

0.15 

0.048 

Ship 

A 4 6 8 T 

A 

26 

59 82 


0.0057 

FNOC 

5 A 7 5 2 

5 

37 

75 5A 

11 



(c) South Temperate (60 # S-30'S, 0*-360*E) 


Category/ 

Hypothesis 


Percent/Test Value 




1 

2 

3 

A 


1 

2 

AVHRR 

1 


F 

T 

T 

1 


0.035 

SMMR 

2 

8 


F 

T 

2 

83 


HIRS 

3 

A 

8 


T 

3 

Al 

83 

Ship 

A 

3 

2 

7 


A 

21 

21 


3 

0.0057 

0.11 


1 1 

0.00098 

0.0033 

0.018 


63 



Table 7-6. SPRED Statistics on 10* Binned Data. (S'ie Table 7-5 for 
Explanation of the Organization of This Table.) 

**SPRED Statistics** 


(a) North Temperate (30*N-60*N, 0*-360*E) 



Category/ 

Hypothesis 







Test 


Percent/Test Value 



1 2 3 4 5 


1 2 

3 

4 

5 

AVHRR 

1 T T F T 

1 

0.12 

0.055 

0.11 

0.033 

SMMR 

2 6 F F F 

2 

63 

0.47 

0.58 

0.27 

HIRS 

3 7 8 T F 

3 

71 88 


0.053 

0.19 

Ship 

4 8 8 7 F 

4 

88 88 

73 


0.29 

FNOC 

5 5 8 8 8 1 

5 

49 88 

88 

88 


(b) Global Tropics 

(30*S 

-30 *N, 0*-360*E) 





Category/ 

Hypothesis 







Test 


Percent/Test Value 




1 2 3 4 5 


1 2 

3 

4 

5 

AVHRR 

1 T T T T 

1 

0.051 

0.012 

0.19 

0.25 

SMMR 

2 5 T T T 

2 

57 

0.015 

0.0063 

0.030 

HIRS 

3 16 FT 

3 

1.0 64 


0.20 

0.0050 

Ship 

4 7 18 T 

4 

67 1.0 

82 


0.0096 

FNOC 

5 5 3 1 5 1 

5 

53 22 

1.0 44 


(c) South Temperate 

(60*! 

S-30*S, 0*-360*E) 





Category/ 

Hypothesis 







Test 


Percent/Test Value 




12 3 4 


1 2 

3 . 

4 


AVHRR 

1 T T F 

1 

0.012 

0.17 

0.076 


SMMR 

2 6 FT 

2 

66 

0.46 

0.032 


HIRS 

3 4 8 T 

3 

41 83 


0.016 


Ship 

4 8 15 

4 

83 1.0 

50 
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The similarity of 10* binned results and 2* binned results indicates that 
the sensitivity of the PPP test is not affected by changing bin size on these 
scales. This implies that the noise problems cannot be overcome easily by 
simply gathering more data. 
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SECTION VIII 


N86-1Q&62 


REQUIREMENTS FOR SEA SURFACE TEMPERATURE GROUND TRUTH 
IN THE INDONESIAN REGION 

J. Penrose 

Western Australia Institute of Technology 


The comparatively low density of ship and XBT observations in large areas 
of the Southern Hemisphere and the tropics limits the extent to which satel- 
lite SST estimates can be validated on a global basis. The additional diffi- 
culties associated with previous SST estimates in the tropics have been 
recognized for some time (Barnett et al., 1979) and are still residually 
apparent in the region centered on the Indonesian archipelago (JPL 1983, p. 
3-4) where both extensive cloud cover and high atmospheric water vapor content 
occur, notably in the October to April periods of each year. These factors 
affect AVHRR retrievals particularly and substantially degrade SST returns 
which, in drier areas appear to give useful agreement with ship derived 
values. SKMR returns and to a lesser extent HIRS returns are affected by the 
distributed nature of the land mass in the Indonesian region. AVHRR has the 
advantages in this region of high sampling density and small footprint. 

Present methods of SST retrieval applied in this region provide apparent high 
positive temperature anomalies where cloud cover allows retrievals to 
proceed. While this is presumably a consequence of over correction for 
atmospheric water vapor effects, adequate analysis is difficult to achieve 
because of limited ground truth availability in the region. 

Two related problems thus arise iu evaluating AVHRR SST retrievals in the 
Indonesian region. One concernc the regional nature of the processing appar- 
ently required for AVHRR data and the second concerns the requirements and 
availability of ground truth information. AVHRR regression coefficients 
derived in one set of regional environmental conditions can result in large 
SST errors when applied to data in other regions (JPL 1983, p. 5-3). Pending 
a fuller understanding of the global atmosphere it may bo necessary to develop 
regional processing algorithms for multichannel water vapor corrections, a 
process particularly called for in the Indonesian region. This in turn calls 
for improved ground truth data sets both of SST and atmospheric water vapor. 

The small AVHRR footprint is better suited to comparison with ship point 
measurements than the larger HIRS and SMMR footprints, SST retrievals from 
which are better compared with binned ship data. The nature of shipping in 
the Indonesian region suggests that point, rather than binned ship data are 
more likely to arise. This raises a difficulty already experienced in working 
with differently conditioned ground truth data sets. Reynolds (1S82) nns 
indicated that the widespread availability of ship near-surface temperature 
estimates, notwithstanding the substantial noise in such data, makes It more 
valuable than the more precise, but less dense XBT data. Further, single 
point observations of SST are usually assumed to be representative of an 
ensemble mean value over the spatial dimensions characterizing one or more 
pixels (JPL 1983, p. 5-7). This assumption may require examination in parts 
of the Indonesian region where comparatively large horizontal temperature 
gradients may be expected. 


8-1 


/ 


»• 


Reconngndattong : 

(1) Attention ohould be given to developing regional multichannel pro- 
cessing techniques applicable to the tropical region centered on 
Indonesia. This will call for improved ground truth data sets in 
the area. 

(2) Because the spatial density of ground truth data sets in the 
Indonesian region can be expected to remain low, attention should 
be given both to the development of high accuracy in ground truth 
wherever possible and to the statistical significance of compari- 
sons made with the data sets that do become available. 

(3) Wherever possible, S5T ground truth data in the Indonesian region 
should include measurements illustrative of the spatial character 
of the temperature field, notably over length scales comparable 
with the AVBRR footprint. 

(A) Cooperation in items (l)-(3) above should be sought from scientists 
and institutions active in the Indonesian region in order to maxi- 
mize the amount of ground truth data made available. 


8-2 



SECTION IX 


REFERENCES 


Barnett, T. P., W. C. Patzert, S. C. Webb, and B. R. Bean, 1979: Climatolog- 

ical usefulness of satellite determined sea surface temperatures in the 
tropical Pacific. Bull. Am. Hat. Soc. . 60 , 197-205. 

Chaline, M. T., N. L. Evans, V. Gilbert, and R. D. Haskins, 1984: Requirements 

for a passive IR advanced moisture and temperature sounder. Appl. Opt. . 
23, 979-989. 

COSPAR, (1985): International workshop on satellite derived sea surface tem- 

perature for global climate applications. Report of workshop held 
Hay 29-31, 1985, WashingtJn, D.C. (in preparation). 

Deschamps, P. Y. and R. Frouin, 1984: Large diurnal heating of the sea sur- 

face observed by the HCMR experiment. J. Phys. Oceanog. . 14 , 177-184. 

Garret, J., 1981: The performance of the FGGE drifting buoy system. Adv. 

Space Res. . 1, 87-94, Pergamon Press, Oxford. 

Holl, M. M., 1981: Alternating parallel analysis. In: T he Scientific/ 

Technical Report for the Expanded Ocean Thermal Structure Analysis System 
(Contract 4 i N0028-80-c-y751, U. S. Navy Fleet Numerical Oceanography 
Center), Meteorology International Inc., Monterey, CA. 

JPL, 1983: Satellite-derived sea surface temperature: Workshop-I. JPL Publi- 

cation 83-34 , Jet Propulsion Laboratory, California Inst, of Technology, 
Pasadena, CA, 147 pp. 

JPL, 1984: Satellite-derived sea surface temperature: Workshop-II. JPL Pub- 

lication 84-5 , Jet Propulsion Laboratory, California Inst, of Technology, 
Pasadena, CA, 145 pp. 

Legeckis, R. and W. Pichel, 1984: Monitoring of long waves in the eastern 

equatorial Pacific 1981-83 using satellite multi-channel sea surface tem- 
perature charts. NOAA Tech. Report, NESDIS 8 . NOAA/NESDIS, Washington, 

D.C. 

McClain, E. P., 1981: Multiple atmospheric window techniques for satellite- 

derived sea surface temperatures. Oceanography from Spac e, Plenum Press, 
New York, pp. 73-85. 

McClain, E. P., W. G. Pichel, C. C. Walton, Ahmad, and J. Sulton, 1983: 

Multi-channel improvements to satellite-derived global sea surface tem- 
peratures. Adv. Space Res. , 2, 43-47, Pergamon Press, Oxford. 

Njoku, E. G., 1985: Satellite derived sea surface temperature: Workshop com- 

parisons. Bull, Am. Met. Soc. , 66 , 274-281. 

Preisendorfer, R. and T. P. Barnett, 1983: Numerical model-reality intercora- 

parison tests using small-sample statistics. J. Atm. Sci. . 40 , 1884-1896. 


9-1 


Reynolds, R., 1982: A monthly averaged climatology of sea surface temperature. 

NOAA Tech. Report NWS 31 , KOAA/NWS, Silver Spring, MD, 35 pp. 

Reynolds, R. W. , 1983: A comparison of sea surface temperature climatologies. 

J. Clim. Appl. Meteoral. , 22 . 447-459. 

Robinson, I. S., N. C. Wells, and H. Chamock, 1984: The sea surface thermal 

boundary layer and its relevance to the measurement of sea surface tem- 
perature by airborne and spaceborne radiometers. Int. J. Rem. Sens. . 5, 
19-45. 

Saur, J. F. T. , 1963: A study of the quality of sea water temperatures 

reported in logs of ships' weather observations. J. Appl. Meteorol, . 2, 
417-425. ' " 

Schwalb, A., 1978: The TIROS N/NOAA A-G satellite series. NOAA Tech. Memo. 

NESS 95 , NOAA/NESDIS, Washington, D.C., 75 pp. ~ 

Smith, W. L., 1968: An improved method for calculating tropospheric tempera- 

ture and moisture from satellite radiometer measurements. Mon. Wea. 

Rev. . 96 . 387-396. 

Smith, W. L. , 1930: Operational sounding algorithms. VAS Demonstration Sound- 

ing Workshop. NASA Conf. Publ. 2157 , NASA/Goddard Space Flight Center, 
Greenbelt, MD, pp. 1-10. 

Smith, W. L. and H. M. Woolf, 1982: Algorithms used to retrieve surface skin 

temperature and vertical moisture and temperature profiles from VISSR 
Atmospheric Sounder (VAS) radiance observations. Proc. AMS 4th Conf. 
Atmos. Rad. . Torontc, Canada, American Meterological Soc., Bocton, MA. 

Stowe, L. L., 1984: Personal communication (address: Climate and Earth 

Sciences Laboratory, NOAA/NESDIS, E/RA11, Washington, D.C. 20233). 

Strong, A. E. , L. L. Stowe, and C. C. Walton, 1983: Using the N0AA-7 AVERR 

data to monitor El Chichon aerosol evaluation and subsequent sea surface 
temperature anomalies. Proc. 17th Int. Symp. Rem. Sens. Env. . Univ. of 
Michigan, Ann Arbor, MI, pp. 107-122. 

Strong, A. E. and E. P. McClain, 1984: Improved ocean surface temperatures 

from space - comparisons with drifting buoys. Bull. Am. Met. Soc. , 65 . 
138-142. 

Strong, A. E., 1984: Use of drifting buoys to improve accuracy of satellite 

sea surface temperature measurements. Trop. Ocean-Atm. Newsl. , 25 , 16-18. 

Susskind, J., J. Rosenfield, D. Reuter, and M. T. Chahine, 1984: Remote 

sensing of weather and climate parameters from HIRS2/MSU on TIROS-N. J. 
Geophys. Res. , 89 , 4677-4697. 

Tabata, S., 1982: An evaluation of the quality of sea surface temperatures and 

salinities measured at Station P and Line P in the northeast Pacific 
Ocean. J. Geophys. Rea .. 17 , pp. 374-385. 


9-2 



Walton, C. C., 1980: Deriving sea surface temperatures from TIROS-N data. 

Remote Sensing of Atmospheres and Oceans . Academic Press. New York 
pp. 547-578. 

Walton, C. C., 1985: Satellite measurements of sea surface temperatures in 

the presence of volcanic aerosol. J. Cllm. Appl. Meteoral .. 24, 501-507. 

WCP, 1984: Report of the TOGA workshop on sea surface temperature and net 

surface radiation. La Jolla, CA, March 1984. World Climate Program 
Report WCP-92 . WHO, C.P. 5, 1211 Geneva 20, Switzerland. 

Wilheit, T. T. , J. Gatlin, D. Han, B. M. Krupp, A. S. Milman, and E. Chang, 
1983: Retrieval of ocean surface parameters from the scanning multi- 

channel microwave radiometer (SKMR) on the Nimbus-7 satellite. IEEE 
Trans. Geosci. Rem. Sens. . 15, 225-224. 






9-3 


jr.. * 


i 


-• * ^ 




,y . 


APPENDIX A 
CORRELATION TABLES 




) 


(.1 


ft * * * '^*' ** (pi-- ■’■ »™ Y*M iCHl *■ J£wi .. uit,W*4fc- 'S»A<CJ7t. Mi'j ii "■« Ilui. ----- - 4 « — — t- t— -wA. * ". !a - - n] i»w|. f T^ i -^ i-tff SS iiVK- tfrfiS 




<M, kCL 


Table A-l. Statistical Conparlsons Between 2* Latitude-Longitude Binned 
Masked up to 600 ka fro® Land: (a) Kovecber 1979, Global. 


SST Ano&aly Fields. All Data 
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Statistical Comparisons Between 2* Latitude-Longitude Binned SST Anomaly Fields with 3x3 
Cell Spatial Smoothing. All Data Masked up to 600 1cm from Land: (a) November 1979, Global. 
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(g) December 1981, North Pacific 
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. 8.57 

C. 

. 8.88 

C. 

• 04 

C. 

0.88 

AVHK8 




B. 

• 51 

B. 

8 21 

B. 

-8.58 

B. 

1.29 

B. 

0.48 

D. 

8.87 





8. 

. • 39 

9. 

. 8 37 

8. 

. 8.49 

8. 

. 8 28 

6. 

• 37 

8. 

8.61 





N. 

226 

N. 

236 

M. 

1S8 

H. 

48 

H. 

1S7 

H. 

234 

8USSS1KD 






C. 

• 97 

C. 

. 8.44 

C. 

. 8.78 

C. 

8.85 

C. 

8 88 

HI X3 






fi. 

-8 32 

B. 

-1 97 

8. 

8.74 

B. 

-8.84 

B. 

8.48 







S. 

. 0 14 

S 

8 42 

8. 

a 8. IS 

8. 

• 36 

0. 

8.37 







N 

311 

H. 

231 

M 

48 

N. 

IS7 

N. 

311 

5USLKIND 








C. 

. 8.49 

C 

. 8 87 

C.. 

8.88 

C. 

8.88 

HITS 








B. 

-8.67 

8 

8.89 

B. 

• 27 

B. 

8 72 

WEIGHTED 








S. 

. 8.44 

S 

. 8 12 

8. 

8 31 

8. 

8 47 









N. 

231 

N 

48 

N. 

IS7 

N. 

311 

HU HAH 










C 

. 8.61 

C. 

8.66 

C. 

8.88 

SHH« NIGHT 










B 

1.71 

B. 

1 . 87 

B. 

1 38 

CELL 234 










6 

. 8.35 

8. 

8 St 

S 

• 42 











N 

38 

N. 

183 

M. 

231 

DATES 












C. 

8 42 

c. 

8 88 

VAS 












N. 

-8 5S 

B. 

-8 79 













S 

• 2? 

8 

8.24 













M 

30 

N. 

48 


PAZAN 

SHIPS 

>S/CELL 


C. . Ml 

B. 0 41 

ft . • 4.4 
N . 157 


clihatoi 
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Table A-3. Sane as Table A-2 Except HIRS and KIRS Version 2 Replaced by HIRS Version 2 Weighted and HIRS 
Version 2 Weighted Smoothed. (All Other Data Sets Smoothed.) (a) March 1982, Global. 


Twi Degree Average SSI Anomaly Crest Correlations 
Region Sb 0 deg N to S6.0 deg S, M deg K to 360 0 deg C 


At leewt 1 observe tlon(s> et common grid points 
<C--Cres* Correlation, B— Bias, S— Standard Deviation 
aboot bias, N — Number »f cannon valid ebser vat lens) 

Bias is average ever column temperatures nin«t raw tenptrottrot 



NOAA 

NOAA 

NOAA 

SUSSK1ND 

SUSSKIND 

MILMAN 

WlLHElV 

BATES 

PAZAN 

PAZAN 

TRANSPAC 

CLXMATOL 


AVHPR 

AVHHH 

AVHRR 


HIRS 

HIRS 

SMMR NIGHT HAN 

VAS 

SHIP3 

SHIPS 


XBT 




DAT 

NIGHT 





VERSION 

2CFLL 234 

SHMR/SHI? 


>8/CELl 

)?0/LLLL 






MAR 








COMPOSITE 








NOAA 


C 0 93 

c 

. . 8 90 

C. 

. 0 SB 

C . 0.4? 

C . . 1.62 

C. 0.60 C 

0.17 

C . . 1.74 

C 0.21 

C. 

. 0.91 

C. . 

0.00 

AVHRR 


D . -0.37 

B 

-0 18 

B. 

-0 S3 

B . -0 52 

B. -0 37 

B -0.21 B 

0 24 

B. 1 11 

B 1 13 

B. 

-0 2? 

B. 

-0 40 

DAT 


S t 31 

S 

» 17 

S 

8 71 

S . 0 66 

S . 0 77 

S . 0 54 8. 

0 65 

S. . 0 33 

8 . 0 21 

S 

. 0.27 

S. . 

0 64 

MAR 


N . 3383 

N 

. . 3404 

H . 

. 3404 

N . . 3404 

H . 2311 

N. 1666 N 

291 

N. 23? 

N. . y 

N. 

9 

N. . 

3404 


NOAA 

C. . 

o ?y 

c. 

. 0 45 

C 

. 0.4? 

C. , 

0.65 

C. . 

0.71 

C. . 

0 20 

C. . 

0 62 

C . 

0 94 

c . 

0.75 

C. . 

0 00 

AVMRR 

B. 

0 17 

B. 

-0 06 

b 

-0 06 

B 

9 97 

B. 

• 27 

b 

0.64 

B. 

1.57 

B 

0 67 

B. 

0.2S 

B. 

0 04 

NIGHT 

S . 

0 18 

S. 

. 0 72 

s. 

. 0 66 

3. 

0 72 

S . 

• 40 

S 

0 59 

S. . 

0 23 

S 

0 23 

8 . 

1 32 

S. . 

0 72 


N. , 

. 3921 

N. 

. 3921 

N. 

. 3?2t 

N. 

. 2807 

N. 

2103 

H. . 

311 

N. . 

368 

M. . 

23 

N . . 

IB 

M. . 

3921 


NOAA 

C. 

0 sS 

C. 

0 SO 

C. 

0 65 

C 

9.66 

C. 

9 21 

C. 

0.77 

C 

• 92 

c. 

• 70 

C. 

1 00 

AVHRR 

B. 

-0.24 

B. 

-0 24 

B 

-Q.I2 

B. 

0.08 

B 

• 51 

B. 

• 44 

B. 

• 60 

B. 

0.14 

1. 

-0.14 

. < 

S. 

0.76 

S. 

• 6? 

S. 

0.73 

8. 

• 48 

8 

0 59 

S. 

0 2? 

S 

• 28 

8. 

• 38 

8. 

• 79 


N. 

396? 

N. 

398? 

N. 

2G41 

it. 

2153 

H. 

201 

N. 

368 

N. 

23 

M. 

18 

N. 

398? 


SUSSKIND 

C . 

0 93 

C. 

0 17 

C 

0.15 

C. . 

0.48 

C. . 

0 46 

C . 

. 0 11 

C. 

. 0.40 

c. . 

• 00 

HIRS 

B. 

0 60 

B 

0 02 

B 

0.17 

1. 

• 41 

B. 

-0.29 

B. 

-• JB 

B. 

-0.76 

B. 

8 10 


S 

0 2S 

S . 

.1.00 

S . 

C.6H 

8. . 

0 .SI 

8. , 

. 0.41 

S . 

.0.47 

S 

. 0 52 

S . 

1 64 


H. . 

4021 

H 

2867 

N. , 

2190 

h. . 

Jli 

N. , 

368 

N. 

23 

N. 

18 

* . 

4021 

SJSSKINP 



C. . 

0.17 

C. . 

0.14 

C. . 

0.49 

C. . 

1 55 

C. 

. 0.30 

c. 

. 0 01 

C. . 

• 08 

HIRS i . . 



B. 

0.05 

B. 

5.21 

D. 

0.52 

B. 

-0 2? 

B. 

-0 27 

B. 

-I B2 

B. 

0.10 

VERSION 2 



S. . 

; •*> 

Z. . 

0.61 

S. . 

• 4V 

8. . 

O 31 

3.. 

. 0 32 

S 

. 0 51 

S . 

• 54 




N. 

' 2867 

N. 

. 2191 

N. . 

301 

N. 

368 

N. 

23 

N 

ta 

N. . 

4021 


MILMAN 


C. . 

0 02 

C. 

-c e4 

C . 

0 15 

c. 

-0 16 

C. 

. o.sv 

C 

0 00 

SNMR NIGHT 


B 

r 28 

8. 

0.8? 

B. 

• 17 

B. 

-0.8 

B 

-0 54 

a. 

-0.04 

CELL 234 

/ ’ * 

s. . 

0.39 

S 

. 0 83 

S . 

0 79 

S. 

. 0 4a 

S 

. 0 84 

s 

. 1 94 



N. . 

1618 

N. 

258 

N. 

300 

N. 

21 

N. 

18 

N 

<007 


UILHEIT/ 

C. 

•Mb 

C. 

. 0 75 

c. 

. 0 SI 

c. 

. 1 78 

C. 

. 1.00 

HAN 

b. 

6.3? 

». 

-0 07 

a 

-0.15 

B. 

-0 46 

I . 

-O 19 

SMMk/SHIP 

8. 

. 0.61 

s. 

. 0 25 

B. 

. 0 IS 

9. 

. 0 31 

I . 

. 0.46 

COMPOSITE 

N. 

234 

N. 

207 

N 

14 

M 

18 

h. 

. 2190 


Bms 

C. 

.2.7? 






. 0.00 

VAS . . . . . . , 

B. 

-O 91 






-0 60 


3 

. 0.26 






1 86 


M. 

51 





v 

301 

PAZAN 



C. . 

1 00 

C. 

. 0.70 


4.10 

SHIPS 



6. 

1 It 

B. 

-0 47 

’• 

13 

>S/CELL 



S. . 

• C0 

S 

. 1 3V 

V 

r 5 




N. . 

23 

H 

18 

N. 

. \ 

PAZAN .. .. 





C 

. 0.98 

C. 

. 0 0) 

SHIPS 





B 

-0 49 

B 

0 01 

>20/CEl.L 





S 

. 1 IS 

S. 

. 0 21 






N. 

4 

K. 

23 

TRANSPAC 







C. 

. I II 

XBT 







B. 

0.77 








3. 

. 1 48 








N. 

18 


CL IKATOl 





Table B-l. 


I 


w 

i 


I 


Partitioned rms Error Between Triplet Combinations of the 2 * - Binned SST Anomaly Data Sets 
AVHRR, SKMR, HIRS, and Climatology: (a) November 1979, AVHRR. 

RUS ERROR TABLE FOR DAT ft SET t 


FILE NAME: SSTSHOP: ISSTlftfiE6O004.BLV 

DESCRIPTION: NORA RUHRR 

LATITUDE BOUNDS* -55 TO 55 LONGITUDE BOUNDS* 0 TO 360 
OUERALL RUERRGE RMS ERROR- 0.53 


1 2 3 


1 



2 

3 

4 


3772 

3773 


0.34 


4538 


0.59 

0.67 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIRO 

LCUER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

ROW AIC COLUMN NUMBER OF TABLE DEFINE THE 2ND RNO 3RD DATA SETS IN TRIfd 


RON /COLUMN FILE NAME 


DESCRIPTION 


1 SSTSHOP : ISST IRAE60004 . ELV 

2 SSTSHOP. ISST] SDC60604. EL V 

3 SSTSHOP: ISST 3 he 966804 . BLV 

4 SSTSHOP: lSSTJtgz6O0.bin 


NORA AUHRR 
SNMR CELLS 234 
HIRS/NSU2 
CLIMATOLOGY 



(b) November 1979, SKHR 


I 


I 


V 


* 

•-1 


\ 


RfIS ERROR TABLE FOA DATA SET 2 




DESCRIPTION: SmR CELLS 234 ” 

LATITUDE BOUNDS- -33 TO 55 LONG I TUBE BOUNCE 


OUERRLL RUERRGE RMS ERROR- 0.67 

1 2 3 

4 

1 * ■ 0.01 

0.77 

2 —mm*. m , 

3 3772 

0.94 

4 3773 3649 

m 


0 TO 360 


"'™ R,CHr ^ cormr* Rns eb** 

“ LEFT TBImcLE v ,R8LE hoses com. mo po.nts m 

^ PHD COL™ SOSES of TRBLE CCE.HE W 2* M, 


ROU /COLUMN FILE NAME 


DESCRIPTION 


1 

2 

3 

4 


SSTSHOP: ISS71SDC60804. 81 

SSTSHOP: ISSTlhcg£ 60 O 4 gj 
SSTSHOP; ISSTJtgz6O0.bin 


NORA RUHRR 
STOW CaLS 234 
HIRS/MSU2 
aiMATOLOGY 


TRIAO 

TRIAD 






(d) November 1979, Climatology 


i 

! 


I 

03 

I 

U> 

•5 


t 

i 


i 


RMS ERROR TRBLE FOR ORrfl SET 4 


FILE NAME: SSTSHOP: tSSTltgz600.bin 

DESCRIPTION: CLIMATOLOGY 

LRTITUOE BOUND S= -55 TO 55 LONGITUDE BOLM)S= 0 TO 360 
OVERFILL ROERfiCE RMS ERROR= 0.63 


1 2 3 


t 


0.70 0.76 


4 


2 3773 0.45 

3 4533 3349 


4 


UPPER RIGHT TRIRNGLE OF TRBLE CONTAINS RNS ERROR IN TRIRO 

LOWER LEFT TRIRNGLE CF TRBLE CONTAINS NUMBER CONNON GRID POINTS IN TRIAD 

ROW RNO COLUMN NUfiBER OF TRBLE DEFINE THE 2ND RMO 3R0 DRTR SETS IN TRIRO 

ROW /COLUMN FILE NAME DESCRIPTION 


1 SSTSHOP : ISST lRfiEC-0304 . BLY 

2 SSTSHOP: ISST 1S0C6C904. BLY 

3 SSTSHOP : ISST )hcs603Q4 . BLY 

4 SSTSHOP : ISST J tgzCea .bin 


NORA AOHRR 
SMMR CELLS 234 
HIRS/MSU2 
CLIMATOLOGY 



(e) December 1981, AVHRR 


03 

» 

CT*> 


RflS ERROR TABLE FOR DATA SET I 


r I IU14 IC : 


onur : I oo i J Mttt 00003 BLV 
DESCRIPTION: NORA AVHRR 

LATITUDE BOUNDS* -55 TO 55 LONGITUDE BOUNDS* 
OUERALL AVERAGE RflS ERROR= 0.4Q 


0 TO 360 


1 2 



3840 

3840 


3 4 


0-38 0.45 

0.72 

4612 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIAD 
LGUER LEFT TRIANGLE OF TABLE CONTAINS NUMBER CGI mi GRID POINTS 
ROW AND COLUMN NUMBER OF TABLE DEFINE THE 2ND AND 3RD DATA SETS 


ROU/COLIJMN FILE NAME 


DESCRIPTION 


1 

2 

3 

4 


SSTSHOP : ISST 1PAE60O03 . BLV 
SSTSHOP: ISST 1SDC60603 . BLV 
SSTSHOP: ISST lhcg6O0G3. BLV 
SSTSHOP: ISSTJtgz683.bin 


NORA AUHRfi 
SI1MR CELLS 234 
HIRS/MSU2 
CLIMATOLOGV 


TRIAD 

TRIAD 



f 


r-B 




(f) December 1981, SHMR 


RMS ERROR TABLE FOR DATA SET 2 


FILE NAME : SSTSHOP : [ SST J SOCO00O3 . BLV 

DESCR i FT I ON : SMMR CELLS 234 

LATITUDE BOUNDS* -55 TO 55 LONGITUDE BOUNDS* 


0 TO 360 


OUERALL AVERAGE RMS ERROR* 0.95 


0.92 0.86 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIAO 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUMSER COMMON GRID POINTS IN TRIAD 

ROU AND COLUMN NUMBER OF TABLE DEFIfC THE 2f€ AND 3RO DATA SETS IN TRIAD 


ROW /COLUMN FILE NAME 


DESCRIPTION 


SSTSHOP: l SST 1 RRE6OO03 . BLV 
SSTSHOP : l SST 3S0C60803 . BLV 
SSTSHOP: IS3TlhcQ60e.03.BLV 
SSTSHOP: fSSTltqr609.bin 


NCAA RUHRR 
StUIR CELLS 234 
H1RS/MSU2 
CLIMATOLOGV 



B-8 


> - \ 

-'N 



/ V 


y ■ 


i prr77'T i 7 ' t r's r 1 1 1 * 

■ ' .•■ j . '» - 


~T??T 

t .* 


(g) December 1981, HIRS 

RMS ERROR TABLE FOR DATA SET 3 

F I LE NAME : S3TSHOP : ISST 1hcg60803 . BLV 

DESCRIPTION: HIRS/MSU2 

LATITUDE BOUNDS 3 -55 TO 55 LONGITUDE BOUNDS 3 0 TO 300 
OUERfiLL RUERRGE RI1S ERROR- 0.53 

I 2 3 

1 0.74 

2 3848 

4 4612 3852 

UPPER RIGHT TRIANGLE OF TABLE CONTAINS BUS ERROR IN TRIAD 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

ROW AND COLUMN NUMBER OF TABLE DEFINE THE 2ND AND 3RD DATA SETS IN TRIAD 


4 

0.40 

0.46 


ROW/COLUMN FILE NAME 


DESCRIPTION 


1 SSTSHOP : ISST 1AAEC-CO03 . BLV 

2 SSTSHOP : ISST 1SDC68803 6LY 

3 SSTSHOP : I SST ) hcg60003 . BLV 

4 SSTSHOP: ISSTltgzC00.bin 


NORA fiVHRA 
SMMR CELLS 234 
HIRS/MSU2 
CLIMATOLOGY 


cS.1 






n« >\ % 


*rr 

,V~L — 


to 

i 

\o 




< : . *\ ■ ■ 

{ ' v; ,w v r- ? .Tv y -vr -- * - 

’ ■ ■ ,‘. f • .j,:l'\r' - " . '» 


■.r't'Tm g V t i .Tr 


'•. <:l r A 


(h) December 1981 t Climatology 
RT1S ERROR TABLE FOR DRTR SET 4 


F I LE NRME : SSTSHOP : lSSTltgz600.bin 

DESCRIPT I OH: CLIMATOLOGY 

LRT I TUOE BOUNDS* -55 TO 55 LONGITUDE BOUNDS* 0 TO 360 


I* 

e 

e 

v 

* 

1 


H 

Y, 


OVERALL AVERAGE R*1S ERBOR= 6.50 


; y 

i 

1 2 3 4 / 

1 . 0.73 0.48 

2 3340 0.3O 


3 4612 3S52 


UPPER RIGHT TRIfiNGLE GF TftBLE CONTAINS RMS ERROR IN TRIRO 

LOWER LEFT TRIftNGLE OF TABLE CONTAINS HUMBER COMMON GRID POINTS IN TRIftO 

ROW AND COLUMN NUMBER OF TRBLE DEFINE THE 2ND AND 3RD DATA SETS IN TRIAD 


ROW/COLUMN FILE NAME 


DESCRIPTION 


1 SSTSHOP : 1 SST 1 RAE6O0C3 . BLY 

2 SSTSHOP : ISST 1SDC60003 . BLY 

3 SSTSHOP: ISST ]hcq60003.BLV 

4 SSTSHOP : ISST 1 tqzCGO .bin 


NORA RVHRfi 
SMMR CELLS 234 
HIRS/MSU2 
CLIMATOLOGY 



' \ 


H -■ -T^TCr 


rrr tt 
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(1) March 1982, AVHRR 


RMS ERROR TRBLE FOR DATA SET 1 


F I LE NAME : SSTSHOP : [SST I AAE600O2 . ELV 

DESCRIPT I OH: MORA A'JHRR 

LATITUDE BOUNDS 3 -55 TO 55 LONGITUDE BOUNDS 3 0 TO 360 
OUERALL AUERAGE RMS ERROR 3 0.60 


12 3 4 

2 0.60 0.62 

3 3780 0.59 

4 378i 4612 


UPPER RIGHT TRIRNGLE OF TABLE CONTAINS RMS ERROR IN TRIAD 

LOWER LEFT TRIRNGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

P.OU AND COLUMN hJMBER OF TABLE DEFINE THE 2ND AND 3R0 DATA SETS IN TRIAD 

ROW /COLUMN FILE NAME DESCRIPTION 


1 SSTSHOP : ISST 1AAE60002 . BLV 

2 SSTSHOP : ISSPISDC60002 . BLV 

3 SSTSHOP : ISST JhaQ66802 . BLV 

4 SSTSHOP: [SSTUgs600.bin 


NORA RUHRR 
SI1MR CELLS 234 
HIRS/MSU2 
CLIMATOLOGV 



(j) torch 1982, SMMR 


RI1S ERROR TABLE FOR DATA SET 2 


F I LE NAME : SSTSHOP : tSST JSDC6O002 . BLY 

OESCR I PT I OM ' SKUA CELLS 234 

LATITUDE BOUNDS- -55 TO 55 LONGITUDE BOUNDS* 0 TO 360 


OUERALL AUERAGE RNS ERROR* 8.98 

12 3 4 

: 0.98 0.97 

3 3780 0.98 


4 3781 3790 

UPPER RIGHT TRIANGLE OF TABLE CONTAINS RNS ERROR IN TRIAD 

LONER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

RON AND COLUMN NUMBER OF TABLE DEFINE THE 2ND AND 3RD DATA SETS IN TRIAD 

RON /COLUMN FILE NAME DESCRIPTION 


1 SSTSHOP : ISST1 RAE60O02 . BLV 

2 SSTSHOP: ISST1SOC60002. BLY 

3 SSTSHOP : ISST lhag60002 . BLV 

4 SSTSHOP; ISSTltgz600.bin 


NORA ADHRR 
Slim CELLS 234 
HIRS/MSU2 
CLIMATOLOGY 


1-12 


/ 


J 


* 




00 
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(k) March 1982, HIRS 


RMS ERROR TRELE FOR DRTR SET 3 


F I LE NAME : SSTSHOP : ISST ]hag6OO02 . BLV 

DESCRIPTION: HIRS/MSU2 

LRTITLIDE BOUNDS* -5? TO 55 LONG I RIDE BOUNDS 1 
OVERALL AVERAGE RMS ERROR- 0.68 


0 TO 360 


J 

12 3 4 '* 

1 0.65 0.72 

2 3780 0.66 | 


4 4612 3798 

UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIRO 


LOWER LEFT TRIRNGLE OF TABLE CONTRINS NUMBER COMMON GRID POINTS IN TRIAD 
ROW AND COLUMN NUMBER OF TRBLE DEFINE THE 2ND AND 3RD DATA SETS IN TRIAD 


ROU /COLUMN FILE NAME 


DESCRIPTION 


1 SSTSHOP : tSST IARE6CC-02 . BLV 

2 SSTSHOP: tSST JSDC6OO02 . ELY 

3 SSTSHOP : tSST ] hciq6O002 . BLV 

4 SSTSHOP : ISST 1 tqz600 .bin 


NORA RUHRR 
SMMR CELLS 234 
HIRS/MSU2 
CLIMATOLOGY 





B-13 
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(1) March 1982, Climatology 
RMS ERROR TfiBLE FOR DATA SET 4 


FILE NAME: SSTSHOP: rSSTHq26e0.bin 
DESCRIPTION: CLIMATOLOGY ' 

LATITUDE BOUNDS 3 -55 TO 55 LONGITUDE BOUNDS 3 0 TO 350 
OVERALL AVERAGE RMS ERROR 3 8.43 


1 2 

1 0.48 

2 3731 

3 4612 3790 


3 4 

0.50 

0.38 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIAD 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

ROW RND COLUMN NUMBER OF TABLE DEFINE THE 2ND fit ID 3RD DATA SETS IN TRIAD 


ROW /COLUMN FILE NAME 


DESCRIPTION 


1 SSTSHOP : t SST 1 AAE600O2 . ELY 

2 SSTSHOP : ISST 1SDC68002 . BLY 

3 SSTSHOP : (SST lhcq6OG02 . BLY 

4 SSTSHOP : l SST ] tqz603 . b i n 


NORA AVHRA 
SMMR CELLS 234 
H I RS/I1SU2 
CLIMATOLOGY 


\ 


•'i'VUTS,-' "Mi.VTai.TVH' 



«7l~ 


l 
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(m) July 1982, AVHRR 


RMS ERROR TABLE FOR OATR SET I 

F I LE HARE : SSTSHOF : ISST JRHE6000 1 . BLV 

DESCRIPTION: NORA RUHRR 

LATITUDE BOUNDS= -55 TO 55 LONGITUDE B0UNDS= 0 TO BOO 
OUERALL AUEP.AOE RMS ERROR = 0.S7 

2 3 4 


0.76 0.83 

2 13S 0.93 

2139 4694 

UPPER RIGHT TRIANGLE OF TABLE CONTAINS RI1S ERROR IN TRIAD 

LOUER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

ROU AND COLUMN NUMBER OF TABLE DEFINE THE 2ND AND 3RD DATA SETS IN TRIAD 

RGU/COLUMN FILE NAME DESCRIPTION 



1 SSTSHOP : ISST 1 AAE6O00 1 . BLV 

2 SSTSHOP : ISST 1SDC60OO t . BLV 

3 SSTSHOP : ISST Ihaa6800 1 . BLV 

4 SSTSHOP: lSSTltgz600.bin 


NORA RUHRR 
SMMR CELLS 234 
H I RS/MSU2 
CLIMATOLOGY 



) 



/ 


ca 

I 




(n) July 1982, SKMR 


RtlS ERROR TftBLE FOR DATA BET 2 


FILE nntlE: SSTSHOP. ISST JSDCG0001 .ELY 

DESCRIPT I OH: SI 111ft CELLS 234 

LATITUDE BOUNDS- -55 TO 55 LOHOITUDE BOUNOS- 0 TO 368 
OUERftLL RUERROE RMS ERROR- 1.01 

I 

2 

3 2138 

4 2139 


2 3 4 

1.85 0.95 

1.03 

2305 


UPPER RIGHT TRIANGLE OF TR8LE CONTAINS RMS ERROR 111 TRIRD 

LOUER LEFT TRIANGLE OF TABLE CONTAINS HUNT "A COMMON GRID POINTS IN TRIAD 

ROU ftNO COLUMN NUMBER OF TABLE OFF I IE THE 2N0 Al© 3f£» DATA SETS IN TRIAD 


R0W/C0LU1N FILE NAME 


DESCRIPTION 


1 SSTSHOP . ISST IRAE6800 1 . BLY 

2 SSTSHOP : ISST 1S0C6C88 1 . BLV 

3 SSTSHOP ISST IhogOCeO 1 . BLV 

4 SSTSHOP : ISST J tgz&OO .tin 


HOAR AUHAA 
SMMR CELLS 234 
HIRS/NSU2 
aiMATOLOGY 
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/ 


./ 



(o) July 1982, SIRS 


Kl$ ERROR TRBLE FOR DATA SET 3 


FILE KATIE: SSTSHOP- [SST Jhcg&OSQI ,BLV 
DESCRIPTION: HIRS/MSU2 

LATITUDE BOUNDS= -33 TO 33 LONGITUDE BOUNDS- 0 TO 360 
OUERftLL AUEREDE RMS ERROR- 0.75 


1 2 3 4 

1 0.74 0.76 

2 2138 0.75 

4 4094 2383 


UPPER RIGHT TRIRNGLE OF TRBLE CONTAINS RHS ERROR IN TRIfiO 

LOWER LEFT TRIfiMOLE OF TRBLE CONTAINS HUMBER COMMOM GRID POINTS IN TRIED 

ROW RNO COLUMN NUMBER OF TRBLE DEFINE THE 2ND AND 3R0 ORTR SETS IN TRIED 

ROU/COLUMN FILE NAME DESCRIPTION 

1 SSTSHOP : CSST lflftE6033 1 . BLY NORA EK'HFFl 

2 SSTSHOP: ISSTJSOC6000 I .BLY SMMR CELLS 224 

3 SSTSHOP: ISST Jhog60O9 I. BLY HIRS/NSU2 

4 SST SHOP : I SST 1 tgz&60 .bin CLIMATOLOGY 




(p) July 1982, Climatology 


RMS ERROR TRBLE FOR DRTR SET 4 


F i LE NAME : SSTSHOP : lSSTltgz600.bin 

DESCR I PT I ON : CLIMATOLOGY 

LATITUDE BOUNDS- -55 TO 55 LONGITUDE BOUt©S= 0 TO j68 
OVERFILL AVERAGE FITS ERROR 3 0.47 

1 2 3 4 

t 0.58 0.45 

2 2139 0.39 

3 4094 2305 

UPPER RIGHT TRIANGLE OF TABLE CONTAINS RNS ERROR IN TRIAD 

LOUER LEFT TRIRNGLE OF TRBLE CONTAINS NUflBER CONIOII GRIO POINTS IN TRIAD 

ROW AND COLUMN NUMBER OF TABLE DEFINE THE 21© AND 3PD DATA SETS IN TRIAD 


ROW /COLUMN FILE NAME 


DESCRIPTION 


1 SSTSHOP : ISST 1RAE6030 1 BLY 

2 SSTSHOP : ISST 1SOC6G00 1 . BLY 

3 SSTSHOP : ISST Jhog600O I . BLY 

4 SSTSHOP : ISST ) tgzC03 .bin 


rtcaa ruhrr 
SMIIR CELLS 234 
HIRS/MSU2 
CLIMATOLOGY 



Partitioned ras Error Between Triplet Combinations of the 2* - Binned SST Anomaly Data Sets 
With 3x3 Cell Spatial Filter Applied: AVHRR, SMMR, HIRS, and Climatology, (a) November 
1979, AVHRR. 


RMS ERROR TABLE FOA DATA SET 1 


FILE KATIE: SSTSHOP: (SST1RAES6000I.BLV 

DESCRIPTION: MQfifi AVHRR 

LATITUDE BOUTOS 3 -55 TO 55 LONGITUDE BOUNDS 3 0 TO 360 
OVERALL AVERAGE RMS ERROR= 0.82 


I 2 3 4 

2 0.70 0.88 

3 H 96 0.86 

4 1502 3089 


UPFER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIRO 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRIO POINTS IN TRIRD 

ROW AND COLUMN NUMBER OF TABLE DEFINE TIS 2ND RNO 3RD DATA SETS IN TRIAD 


FCU /COLUMN FILE NAME 


DESCRIPTION 


t SSTSHO? : l SST 1 AAEs6G00 1 . ELV 

2 SSTEHOP : ISST JSOCsoO80 J . BLV 

3 SSTSKOP : l SST 1 hcjqj6809 1 . BLV 

4 SSTSHOP: [SSTJtgis6C8.bin 


NORA AUHRA 
SMMR CELLS 234 
HIRS/MSU2 
CL I MAT uLOGV 
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(b) November 1979, SMMR 


RfIS ERROR TfiBLE FOR DfiTfl SET 2 


F I LE NAME : SSTSHOP : ISST 1SOCS6O0O 1 . BUY 

DESCRIPTION: SMMR CELLS 234 

LATITUDE BOUNDS 3 -33 TO 33 LONGITUDE BOUNDS 3 0 TO 369 
OUERALL RUERRGE RMS ERROR 3 0.8? 


1 


2 


3 


4 


1 0.98 


0.82 


3 1496 0.82 

4 1502 1846 


UPPER RIGHT TRIRNGLE OF TABLE CONTAINS RMS ERROR IN TRIAD 

LCtlER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRIO FDINTS IN TRIP® 

ROW AND COLUMN NUMBER OF TABLE DEFIffc THE 2ND AND 3RD DATA SETS IN TRIAD 


ROW /COLUMN FILE NAME 


DESCRIPTION 


1 SSTSHOP: ISST )AR£;&8001.BLV 

2 SSTSHOP : ISST JSDCs60O0 1 . BLY 

3 SSTSHOP: ISST )hagsG0001 .BLY 

4 SSTSHOP : ISST I tg2s600 .bin 


NOAA R’.Wfl 
SMMP. CELLS 234 
HIRS/MSU2 
CLIMATOLOGY 




/ . \ 
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(c) November 1979, HIRS 


fins ERROR TABLE FOR ORTfl SET 3 


FILE NAME: SSTSHOP : ISST 1 hogs6000 1 . BLY 

DESCRIPTION: HIRS/MSU2 

LATITUDE BOUT(OS= -55 TO 55 LONGITUDE BOUNO$= 0 TO 360 
OVERALL AVERAGE RMS ERROR* 0.37 

1 2 3 4 

1 0.26 0.37 

2 1496 -- — 0.48 

4 3089 1846 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIRO 

LOWER LEFT TR1RNGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

ROU AND COLUMN NUMBER OF TABLE DEFINE THE 2ND RND 3RD DATA SETS IN TRIAD 


ROW/COLUMN FILE NAME DESCRIPTION 


1 SSTSHOP : ISST ]AAEs6C08 1 . BLY NORA AUHRR 

2 SSTSHOP: ISST 1SDCS&0O01. ELY SMMR CELLS 234 

3 SSTSHOP: ISSTlhcgs6000l. BLV HIRS/MSU2 

4 SSTSHOP : ISST 1 tgzs600 .bin a I MHTOLOGY 


X 



V 

\ 
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(d) November 1979, Climatology 


BUS ERROR TR8LE FOR DRTR SET 4 


F I LE NAI1E : SSTSHOP : tSST 1 tgzs600 .bin 
DESCRIPTION: CLIMATOLOGY 

LATITUDE BOUNDS 3 -55 TO 55 LONGITUDE BGUNDS= 0 TO 360 
OVERALL AVERAGE RMS ERROR 3 0.42 


12 3 4 

1 0.44 0.4? 

2 1502 0.36 

3 3089 1846 


UPPER RIGHT TRIRNGLE OF TRBLt CONTRINS RMS ERROR IN TRIRO 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUI1BER COMMON GRID POINTS IN TRIAD 

ROW AND COLUHN NUMBER OF TABLE DEFINE THE 2>5D RHO 3RD DATA SETS IN TRIAD 


ROW/COLUNN 

FILE NAME 

DESCRIPTION 

1 

SSTSHOP : (SST 1 AAEs6080 1 . BLV 

NORA AUHFift 

2 

SSTSHOP : tSST lSCCs6e00 1 . BLV 

SMMR CELLS : 

3 

SSTSHOP : I SST ]hags6G00 1 . BLV 

HIRS/HSU2 

4 

SSTSHOP: 1SSTJtgzsO0O.bin 

CL 1 HRTOLOGY 


' "**** « 


A..-..,, 
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(e) December 1981, AVHRR 
RMS ERROR TABLE FOR ORTA SET 1 


F I LE NAME : SSTSHO? : [ SST )RREs50802 . BLV 

DESCRIPTION: NORA R'JHRR 

LATITUDE BOUNDS- -55 TO 55 LONGITUDE BOUNDS- 0 TO 360 
OUERRLL RUERAGE RMS ERROR- 0.4? 


1 2 3 4 

2 0 . 48 0.50 

3 2927 0.44 

4 2831 3970 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIAD 

LOUER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIRO 

ROW RNO COLUMN NUMBER OF TRSLE DEFINE THE 2ND RND 3RD DATA SETS. IN TRIAD 

ROW /COLUMN FILE NAME DESCRIPTION 


1 SST SHOP : t SST JAREsSO0e2. BLV 

2 SSTSHOP : ISST JSDCs63882 . BLV 

3 SST SHOP : [ SS'i ] MagsCO082 . BLV 

4 SSTSHCP : l SST ] tgr>603 .bin 


NORA RUHRR 
SMMR CELLS 234 
HIRS/MSU2 
CL I MATOLCGV 


f 

Q-Y 
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(f) December 1981 • SMMR 



/ 


I 

I 


i 


RMS ERROR TRBLE FOR DRTR SET 2 

LATITUDE BOUNDS” -55 TO 55 LONGITUDE BOUNDS” 0 TO 360 
OVERALL AVERAGE RMS ERROR” 0.77 


12 3 4 

t 0.73 0.76 

3 2827 0.77 

4 2831 2663 


UPPER RIGHT TRIANGLE OF TRBLE CONTAINS RMS ERROR IN TRIAD 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

ROW AND COLUMN NUMBER OF TRBLE DEFINE THE 2ND AND 3RD DRTA SETS IN TRIAD 

ROW /COLUMN FILE NAME DESCRIPTION 


1 SSTSHOP : (SST J AAEs680O2 . BLV 

2 SSTSHOP : I SST 1 SDCs68O02 . ELY 

3 SSTSHOP : I SST ] hags68002 . BL Y 

4 SSTSHOP . ISST 1 tgzs600 .bin 


NORA AVHRR 
SMMR CELLS 234 
HIRS/MSU2 
CLIMATOLOGY 



\ V 



n-z 



(g) December 1981, DIRS 

ans ERROR TABLE FOR ORTR SET 3 

FILE NAME: SSTSHOP: lSSTlhaqs68002.BLV 
DESCRIPTION: HIRS/MSU2 

LATITUDE BCL*NOS= -35 TO 53 LONGITUDE BOUNDS- 0 TO 300 
OUERALL AVERAGE RMS ERROR- 0.4? 


12 3 4 

t 0.42 0.55 

2 2827 0.44 

4 3970 2863 



UPPER RIGHT TRIfitICLE OF TABLE CONTAINS RMS ERROR IN TRIAD 

Lot-ten ucrT TmmoLE or thole ccrrrmrre nunoen comon onto roirrre m rmno 

ROW RNO COLUMN NUMBER OF TRSLE DEFINE THE 2ND AND 3RD DATA SETS IN TRIAD 


ROM /COLUMN FILE NAME DESCRIPTION 


1 SSTSHOP: lSST1AAEs60G02.BLY NORA AUHRR 

2 SSTSHOP: (SSTlSOCs60002.BLV SMMR CELLS 234 

3 SSTSHOP: ISSTlhogsf 0082. BLV HIRS/MSU2 

4 SSTSHOP: tSSTltgz36e0.bin CLIMATOLOGY 
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(h) December 1981, Climatology 


RMS ERROR TABLE FOR DRTR SET 4 


F I LE NAME : SSTSHOP : ISST 1 tgzs600 . 

OESCR I PT I OH : CLIMATOLOGY 

LATITUDE BOUNDS 3 -55 TO 55 


bin 

LONGITUDE BOUNDS= 


0 TO 360 


OUERALL AVERAGE RMS ERROR 3 0.38 


T 

2 3 

4 

1 

0.34 0.48 

— 

2 2831 

8.32 

— 

3 3970 

4 

2863 

__ 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIAD 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

ROW AND COLUMN NUMBER OF TABLE DEFINE THE 2ND AND 3RD DATA SETS IN TRIAD 


ROW/COLUNN FILE NAME 


DESCRIPTION 


1 SSTSHOP : tSST 1 AAEs68O02 . BLY 

2 SSTSHOP : ISST1 SDCs60002 . BLY 

3 SSTSHOP: ISST Jhags6OO02. ELY 

4 SSTSHOP : ISST 1 tgzsOOO .bin 


NORA AUHRR 
SMI1P. CELLS 234 
HIRS/MSU2 
CL I MRTOLOGY 
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(i) March 1982, AVHRR 


RMS ERROR TflBLE FOR DRTR SET 1 


FILENAME: SSTSHOP: ISSTlRREs60303.8LV 

DESCRIPTION: NORA AUHRR 

LATITUDE BOL1NOS- -55 TO 55 LONGITUDE BOUfOS- 
OUERALL AUERRGE RMS ERR0R= 0.42 


0 TO 360 


1 

2 

3 


2 3 4 

0,23 0.38 

2841 8.66 


4 


2841 3989 


UPPER RIGHT TRIRNGLE OF TRBLE CONTAINS RMS ERROR IN TRIAD 

LOUER LEFT TRIANGLE OF TR8LE CONTRINS NUMBER COMMON GRID POINTS IN TRIAD 

ROU RND COLUMN NUMBER OF TRBLE DEFINE THE 2ND RND 3RD DATA SETS IN TRIRD 


RON /COLUMN FILE NRME 


DESCRIPTION 


1 SSTSHOP: l£ST]fifl£s60G33.BLV NORR RUHRR 

2 SSTSHOP: [SSTiSOCs63003.BLV SMMR CELLS 234 

3 SSTSHOP: ISSTlhcgs68833.BLV HIRS/MSU2 

4 SSTSHOP : ( SST 1 tgzs603 .bin CL I MATOLOGV 
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(j) March 1982, SMMR 
RMS ERROR TABLE FOR DATA SET 2 


FILE NAME: SSTSHOP: tSSTJSDCs60003 BLV 

DESCRIPTION: SMMR CELLS 234 

LATITUDE 60UNDS= -55 TO 55 LONGITUDE BGUNDS= 


OUERALL RUERAGE RMS ERROR= 0.75 


0 TO 368 


1 

1 

2 

3 2841 

4 2841 


2 3 4 

0.71 0.64 


0.91 

2867 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIAD 
LOUER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS 
ROM AND COLUMN NUMBER OF TABLE DEFINE THE 2ND AND 3RD DATA SETS 


RuU /COLUMN FILE NAME 


DESCRIPTION 


1 SSTSHOP : ISST ]AAEs68003. BLV 

2 SSTSHOP: ISST isocs&eogs. BLV 

3 SSTSHOP: ISST Jhcgs60003, BLV 

4 SSTSHOP: ISST] tgzs600.bin 


NORA RUHRR 
SMMR CELLS 234 
HIRS/I1SU2 
CLIMATOLOGV 


TRIAD 

TRIAD 



(k) March 1982, HIRS 


RMS ERROR TABLE FOR DATA SET 3 


F I LE NAME : SSTSHOP : [SST lhcgs68003 . BLV 
DESCRIPTION: HIRC/MSU2 

LATITUDE BOUNDS= -55 TO 55 LONGITUDE BOUNDS 3 0 TO 360 
OUERALL RUERAGE RMS ERROR- 0.45 


12 3 4 

1 0.67 0.32 

2 2841 0.36 

3 — 

4 3989 2867 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RNS ERROR IN TRIAD 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIAD 

ROM AND COLUMN NUMBER OF TABLE DEFINE THE 2ND AMD 3RD DATA SETS IN TRIAD 


ROW/COLUMN FILE NAME 


DESCRIPTION 


1 SSTSHOP : [ SST J ftfl£s66003 . BLV 

2 SSTSHOP: [SST ]SDCs60603. BLV 

3 SSTSHOP : [ SST lhcg-r&0003 . BLV 

4 SSTSHOP : [SST 1 tgzs600 .bin 


NOAR AUHriR 
SMMR CELLS 234 
HIRS/M3U2 
CLIMATOLOGY 



(1) March 1982, Climatology 



C3 

F-O 

■O 


PMS ERROR T‘ LE FOR DATA SET 


4 


FILE NAME: SSTSHOP: ISST Jtgzs680 bin 

DESCRIPTION. CLIMATOLOGY 

LATITUDE BOUNDS 9 -55 TO 55 LONGITUDE BOUNOS* 0 TO 366 
OVERALL AVERAGE RMS ERROR* 0.4? 


<234 
< 0.69 0.45 

2 234 1 0 . 26 

3 3989 2867 


UPPER RIGHT TRIRNGLE OF TR8LE CONTAINS RMS ERROR IN TRIRO 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID F0 1 NTS IN TRIRO 

POU AND COLUMN NUMBER OF TABLE DEFINE THE 2J*J AID 3RD DATA SETS IN TRIAD 

ROM/COLUMN FILE name DESCRIPTION 


1 SSTSHOP: (SSTlRAEs6eeC3 BLV 

2 SSTSHOP: lS$Ti$0Cs6«ea3 6LV 

3 SSTSHi OP : (SST Jhcgj60O03 . BLV 

4 SSTSHOP: ISSTItgzs600.bin 


NORA AVHF.fi . 
SMMR CELLS 234 
HIRS/MSU2 
CLIMATOLOGY 


s 


/ 


1 
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(a) July 1982, AVHRR 


RMS ERROR TABLE FOR DRTft SET 1 


FILE NAME: SSTSHOP: [S3Tlfift£s6e804.BLV 

DESCRIPTION: MORA RUHrifi 

LATITUDE BOUNDS- -55 TO 55 LONGITUDE BOUICS- 0 TO 360 


OVERALL AVERAGE RMS ERRQR= 0.39 


1 2 


1 


2 


3 4 

0.20 8.44 


3 2661 0.53 

4 2664 3721 


UPPER RIGHT TRIRNGLc OF TR 8 LE CONTAINS RMS ERROR IN TRIAD 

LOWER LEFT TRIANGLE OF TA2LE CONTAINS NUMBER COftlQN GRID POINTS IN TRIAD 

ROW AMD COLUMN NUMBER OF TfiBLE DEFItfc THE 2ND AND 3RD DATA SETS IN TRIfiO 


P.CU /COLUMN FILE NAME 


DESCRIPTION 


1 Si .’SHOP : (SET ]AAE£fcS004 . BLV 

2 SSTSHOP : fSST JSOCs6e084 . BLV 

3 SSTSHOP : (SST lhecs0eS&4 . BLV 

4 SSTSHOP: :SSTltgzs600 bin 


HOAR 

Sit® CELLS 234 

HIR3/H3U2 

CLIMATOLOGV 
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(n) July 1982, SM-IR 

RMS ERROR TABLE FOR DRTR SET 2 

F I LE NAME : SSTSHOP : ISST JS0Cs60894 . ELY 

DESCRIPTION: SMMR CELLS 234 

LATITUDE BOUNDS 3 -55 TO 55 LONGITUDE BOUNDS- 0 TO 350 
OUERftLL fiUEfifiGE R11S ERROR- 0.55 

1 2 3 4 

0.67 0 .55 

2661 0.72 

2664 2653 

UPPER RIGHT TRIANGLE OF TABLE CONTAINS FJ13 ERROR IN TRIAO 

LOWER LEFT TRIfiNGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS IN TRIffi 

F;OU AND COLUMN HUMBER OF TABLE DEFINE THE 2ND AMO 3RD DATA SETS IN TRIAD 

ROU /COLUMN FILE NAME DESCRIPTION 


1 

2 
3 


1 SSTSHOP .ISST] AAEs60004 . BLV 

2 SSTSHOP: tSSTiSDCs50984. ELY 

3 SST SHOP : I S3T Jhcgs60O34 . BLY 

4 SSTSHOP : ISST 1 tgrs660.bin 


NCAA AUHAR 
SMMR CELLS 234 
HIRS/I1SU2 
CLIMATOLOGY 
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(p) July 1982, Climatology 
RMS ERROR TABLE FOR DATA SET 4 


FILE NAME: SSTSHOP: LSSTltgzs600.bin 

Description: Cliiaatology 

Latitude Bounds “ -55 to 55 
Longitude Bounds ■* 0 to 360 
Overall Average rms Error ■ 0.59 


1 0.64 0.69 

2 2664 0.46 

3 3721 2883 


UPPER RIGHT TRIANGLE OF TABLE CONTAINS RMS ERROR IN TRIRO 

LOWER LEFT TRIANGLE OF TABLE CONTAINS NUMBER COMMON GRID POINTS iN TRIAO 

ROW AND COLUMN NUMBER OF TABLE DEFINE THE 2ND RND 3RD DATA SETS IN TRIAD 

ROW /COLUMN FILE NAME DESCRIPTION 


1 SSTSHOP: [SSTJAAE560004.6LY 

2 SSTSHOP : I SST ] SDCs60004 . ELV 

3 SSTSHOP : I SST } hca'60004 . BLV 

4 SSTSHOP : t SST 1 tgzs600 .bin 


NORA AUHRR 
SMI1R CELLS 234 
H I RS/I1SU2 
CLIMATOLOGY 



Table B-3. Average Partitioned rms Error for Each Sensor: (a) No 

Smoothing; (b) 3x3 Spatial Filter Applied to all 
2* - Binned Data. 



November 

1979 

December 

1981 

March 

1982 

July 

1982 

(a) AVHRS 

0.53 

0.49 

0.60 

0.87 

SMHR 

0.87 

0.95 

0.98 

1.01 

HIRS 

0.56 

0.53 

0.68 

0.75 

Climatology 

0.63 

0.50 

0.43 

0.47 

(b) AVHRR 

0.82 

0.47 

0.42 

0.39 

SKMR 

0.87 

0.77 

0.75 

0.65 

HIRS 

0.37 

0.47 

0.45 

0.46 

Climatology 

0.42 

0.38 

0.47 

0.59 


APPENDIX C 
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In situ data distribution for November 1979 showing the 
number of SST measurements per month per 2* latitude- 
Longitude bin for (a) ships, (b) XBT's, (c) FCGE buoys. 
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(a) ships, (b) AVHRR 


Figure C-3 


SST anomalies for December 1981 
(c) MRS, (d) SMMR. 
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Figure C-U. SST anomalies for March 1982: (a) ships, (b) AVHRR, 

(c) HIRS, (d) SMMR. 



Figure C-5, Fcr July, 1982: (a) AVHRR SST data distribution 

(b) AVHRR SST anomalies, (c) ship SST anomalies. 







HIRS SST anomalies for: (a) March 1982 unweighted 
(supplied by J. Susskind), (b) March 1982 weighted 
(supplied by J. Susskind), (c) July 1982 weighted. 
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Figure C-7. SST anomalies for July 1982: (a) ships, (b) AVHRR 

(c) HIRS, (d) SMMR. 
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APPENDIX E 
GLOSSARY 


i 


AMTS 

AVHRR 

FGGE 

FNOC 

GARP 

GLA 

GMT 

GOES 

GOSSTCOMP 

GSFC 

HIRS 

JPL 

LST 

Me IDAS 

MCSST 

MSI 

MSU 

NASA 

NCAR 

NDBO 

NESDIS 

NOAA 

PODS 

SIO 

SMMR 

SR 

SST 

SSU 

STD 

TIROS 

TOVS 

VAS 


Advanced Moisture and Temperature Sounder 
Advanced Very High Resolution Radiometer 
First GARP Global Experiment 
Fleet Numerical Oceanography Center 
Global Atmospheric Research Project 
Goddard Laboratory for Atmospheres 
Greenwich Mean Time 

Geostationary Operational Environmental Satellite 

Global Operational SST Computation 

Goddard Space Flight Center 

Righ-Resolution Infrared Sounder 

Jet Propulsion Laboratory 

Local Standard Time 

Man-Computer Interactive Data Access System 
Multi-Channel SST 
Multi-Spectral Imaging 
Microwave Sounding Unit 

National Aeronautics and Space Administration 
National Center for Atmospheric Research 
National Data Buoy Office 

National Environmental Satellite, Data, and Information Service 

National Oceanic and Atmospheric Administration 

Pilot Ocean Data System 

Scripps Institution of Oceanography 

Scanning Multichannel Microwave Radiometer 

Scanning Radiometer 
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